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Co32aH B WeJLAx MH(opMupoBaHHua YMTAaTeIbCKOM ayTHTOpHn O HOBeMIMIMX AOCTHKCHHAX VU IepCieKTHBaXx B oOsacTu 
MexXaHHKH, MaliMHOCTpOecHHA, MH(OpMaTHKU WM BbIYMCUHTeIbHOM TeXHHKH. V3qaHue ABIIACTCA (bopyMom ASIA 
COTPYTHUYeCTBa poccniickux WM UHOCTpaHHbIxX YUCHBIX, cllocoOcTByeT cOUM2KeHHIO poccniickoro WM MHpoBoro 
Hay4HO-HH(opMallHOHHOTO TIpOCTpaHCTBa. 


2KypHasl BKUOUeH B HepeveHb pelleH3HpyeMbIX HAYYHbIX 131aHH, B KOTOPOM JOJDKHbI ObITh OyOJIMKOBaHbI 
OCHOBHbIe Hay4Hble pe3y.IbTaTbl AUcceprauHi Ha COMCKaHHe y4eHor cTeneHH KaHHAaTa HayK, Ha cONCKaAHHe 
yueHoH cTeneHu OKTOpa Hayk (Ilepeyens BAK) no caeaxyroutuM Hay4HbiM CihewMasIbHOCTAM: 


1.1.7 — Teoperwueckas MexaHnka, JMHaMHkKa MallHH (TexHu4ecKHe HayKH) 

1.1.8 — MexaHnka JedopmMupyemoro TBepsOro Tesla (TexHHYecKHe, (PH3MKO-MaTeMaTH4ecKHe HayKH) 

1.1.9 — Mexanuka 2KHJKOCTH, ra3a H Wa3MbI (TexHu4eCcKHe HayKH) 

1.2.2 — Marematw4eckoe MojesIMpoBaHHe, YMCJICHHbIe METOBI H KOMIVICKCbI IporpaMM (TexHHYecKHe HayKH) 

2.3.1 — CucTemuHbIii anau3, ynpaBsenve u oOpaboTKa HHPopMalluH, CTaTHcTHKa (TexHM4eCKHe HayKH) 

2.3.3 — ABTOMaTH3al[HA HM yIpaBleHHe TeXHOJIOTM4eCKHMH TIpoleccaMH MW NpOH3BOACTBaMH (TeXHM4eCcKHe HayKH) 

2.3.5 — Marematw4eckoe 1 IporpaMMHoe OGecieveHHe BLTHCIMTeJIBHBIX CHCTeM, KOMIUICKCOB H KOMMBIOTEPHBIX ceTeli (TeXHH4ecKHe HayKH) 
2.3.7 — KommbiotepHoe MoyesMpoBaHve VM aBTOMAaTH3alHA NpOeKTHpOBaHHA (TeEXHHMYeCKHe, (PH3HKO-MaTeMaTHYecKHe HayKH) 

2.3.8 — Uacdopmatuka 4 HHOpMallMOHHE!Ie Mpoleccsl (TexHu4eckHe HayKH) 


2.5.2 — Mammmuopeyenue (TeXHHYeCKHe HayKM) 

2.5.3 — Tpeuue 1 43HOC B MalliMHax (TeXHH4eCKHe HayKM) 

2.5.5 — Texnonorna 4 o6opyqoBaHve MexaHHueckoli U PU3MKO-TexHHYecKol OOpaboTKH (TexHHYecKHe HayKH) 

2.5.6 — TexHosIorMa MaluMHOCTpoeHHA (TeXHHYecKHe HayKH) 

2.5.8 — CBapka, poCTBeHHBIe Mpoleccbl U TeXHOMOruN (TexHHyeckHe HaykKH) 

2.5.9 — Merogpl 4 IpHOopbl KOHTpOsIA H WHarHOCTHKM MaTepHasioB, H3]eIMi, BELeCTB H NpHposHO cpebI (TexHM4eCKHe HayKH) 
2.5.10 — Pugpapmuueckne MauIMubl, BakKyyYMHat, KOMMpeccopHasd TeXHuKa, THApO- HW WHeEBMOCHCTeMBI (TEXHHYeCKHe HayKH) 
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ANNIVERSARY OF THE SCIENTIST 
FOBHJIENM YYUEHOTO 


Alexey Nikolaevich Beskopylny, Scientist, Teacher, Expert, Leader, is 65 


Alexey Nikolaevich Beskopylny — Doctor of Engineering Science, 
Professor, scientist and expert in the diagnostics and monitoring of 
building structures, Vice-Rector for Academic and International Affairs, 
Don State Technical University. 

In 1981, Aleksey Nikolayevich graduated from Rostov Institute of 
Construction Engineering majoring in Industrial Transport Engineering. 
He embarked upon a career as an engineer in the industrial transport 
research laboratory. After completing his postgraduate studies and 
defending his PhD thesis, from 1986 to 1994, he worked his way up from 
a Junior research fellow to the Head of the machine reliability management 
laboratory. 

In 1997, A.N. Beskopylny defended his doctoral thesis on the topic 
“Method for determining mechanical properties and quality control of 
structural steels by impact indentation” at Rostov State University of 
Civil Engineering (RSCU). From 1997 to 1998, A.N. Beskopylny 
worked as a Deputy director and, from 1998 to 2005 — as Director of 
the Road Transport Institute, RSCU. From 2005 to 2016, he successfully 
performed the duties of Vice-Rector for educational work, Vice-Rector 


for academic affairs, and Vice-Rector for the organization of educational 


activities at the university. 

As a result of the reorganization in early 2016, RSCU became part of the Don State Technical University (DSTU) as 
the Academy of Civil Engineering and Architecture, and Alexey Nikolaevich was immediately appointed its Director. 
From 2018 to 2021, he worked as Vice-Rector for Academic Affairs and Training of Highly Qualified Personnel, and 
from 2021 to the present — Vice-Rector for Academic Affairs and International Activities of DSTU. 

A.N. Beskopylny made a significant contribution to the modernization and development of the system of training 


highly qualified personnel and strengthening collaboration with foreign partners of the flagship university, including 
within the framework of the “Priority—2030” program. 

A.N. Beskopylny is a well-known scientist in the field of construction, non-destructive testing methods of dynamic 
structures. His teaching and scientific activities are reflected in 240 works, including 3 monographs, 5 textbooks and 
teaching aids, 9 patents. Over the past five years, he has published more than 220 scientific papers, 159 of which are 
indexed in the international scientometric databases Scopus and Web of Science. 

With the active support of Alexey Nikolaevich, the Academy of Publication Activity and a system of scientific, 
organizational and financial support for authors, including young scientists and specialists, were created at DSTU; 10 
scientific journals are published, 6 of which are included in the List of the State Commission for Academic Degrees and 
Titles. The Vice-Rector's initiatives enable the university to fulfill and exceed its publication plans, while developing and 
enhancing the professional competences of authors in the international scientific communication and academic writing. 
As Vice-Rector, Aleksey Nikolaevich made a great contribution to improving the educational process, introducing 
information and computer technologies, training highly qualified engineering and scientific-pedagogical personnel. Under 
his leadership, doctoral and candidate thesis papers were successfully prepared. 

A.N. Beskopylny has developed and is currently implementing scientific and educational programs for the engineering 
master's degree within the framework of 14 educational projects in cooperation with industrial partners of the Don Region: 
PJSC Rostvertol, OJSC PK NEVZ, PJSC Tagmet, LLC KZ Rostselmash, JSC Aluminum Metallurg Rus. 
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The results of Alexey Nikolaevich's scientific and innovative activities have been implemented by the Ministry of Transport 
of the Russian Federation, the Ministry of Transport of the Rostov Region at 53 enterprises in Russia, including Scientific 
Research Institute of Automobile Transport LLC (Moscow), Stroyproekt Engineering Group Association (St. Petersburg), 
the Department of Highways and Traffic Management of the City of Rostov-on-Don, and others. 

Alexey Nikolaevich Beskopylny is the recipient of the Certificate of Honor of the Ministry of Education and Science 
of the Russian Federation (2010), the badge “Honorary Worker of RF Higher Vocational Education” (2003), the badge 
“Honorary Builder of Russia” (2011), the medal “In Recognition of Development of the Construction Industry” (2011). 

The Editorial team sincerely congratulates Alexey Nikolaevich Beskopylny, Editor-in-chief of the journal “Advanced 
Engineering Research (Rostov-on-Don)”, on his 65th birthday and wishes him health, well-being, inspiration and energy 
for new scientific achievements for the benefit of his native university and construction education! 
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Abstract 

Introduction. One of the reasons for undesirable delamination of polymer composites with fabric reinforcement is low 
transverse shear properties. It is known that the reinforcement of polymer fabric composites in the Z direction reduces the 
sensitivity to delamination and increases the viscosity of interlayer fracture. Various methods of three-dimensional 
reinforcement of polymer fabric composites are proposed in the literature. However, they complicate the manufacturing process 
of the structure. The problem is solved by the method of three-dimensional reinforcement proposed in this article — felting. 


This is a local reinforcement of the composite in the Z direction with minimal production changes. The degree of 
Z-reinforcement is determined by the felting density, i.e., the number of needle punches per 1 cm? of the fabric package. The 
work is aimed at evaluating the effect of felting on the interlayer crack resistance of a composite material. 

Materials and Methods. The interlayer fracture toughness Gye was determined on a cross-woven fiberglass with 
felting of 10 cm?. The material was impregnated with Etal-370 resin and Etal-45 hardener. Experiments according to 
ASTM D7905M-14 and GOST 33685-2015 standards were carried out on an Instron 5900R test machine. The stress 
state at the crack tip was analyzed with regard to the nonlocal strength theory in the ANSYS Workbench program (option 
“static strength analysis”). The finite element method (FEM) was used. 

Results. The “load — displacement” curves were considered for the samples. Values Gy. were calculated. The results of 
EMF tests for felting density of 0 cm and 10 cm were summarized. Control samples and felting samples were compared. 
In the latter case, Gu. turned out to be ~33% higher. The stress state at the crack tip was calculated under DCB and ENF 
loading. The dependences of maximum normal and shear stresses, as well as displacements, were visualized in the form 
of graphs and color charts. To get the calculated “load — displacement” dependences using FEM, the reverse method of 
obtaining transverse shear constants was used. DCB loading showed that felting provided increasing the rupture strength 
in the Z direction to ~18%, by 39 to 46 MPa, and in the planes XZ— to ~16%, by 77 to 89 MPa. 

Discussion and Conclusion. Felting as a method of local three-dimensional reinforcement enhances the interlayer crack 
resistance of polymer fabric composites. It provides reducing the area of stratifications after local impacts during the 
operation of structures. Flexible felting technology makes it possible to create zones with an arbitrary impact density, 
increasing fracture toughness only in the required places of structures. The FEM analysis of the stress state at the crack 
tip within the framework of the nonlocal strength theory has shown that in strength calculations, the stratification crack 
can be considered as a stress concentrator. 


Keywords: reinforcement of polymer fabric composites, transverse shear strength, interlayer crack resistance, 
interlaminar fracture toughness, felting local three-dimensional reinforcement 
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AHHOTalna 


Beedenue. O7Ha 43 NpW4WH HeKeNAaTebHbIX PaccOCHHM NOJMMepHbIX KOMIMO3HTOB C TKAHEBbIM apPMHpOBaHueM — HH3- 
KWe TpaHCBepCasIbHO-CABHTOBbIe XapakTepHcTHKH. M3BecTHO, YTO aPMMpOBaHHe MOJIMMepHbIX TKAHEBbIX KOMIO3HTOB B 
HarpaBuleHHu Z yMeHBIULaeT YYBCTBUTEJIBHOCTh K pacCJIOCHHI0 H MOBbILMAeT BASKOCTbh MEXKCIIOHHOTO paspymieHua. B smTe- 
paType lpeyyiararorca pa3Hble ClocoObl TpexMepHOrO aPMHPOBaHHA MOIMMePHbIX TKAHEBbIX KOMIIO3HTOB. OZHaKO OHM 
YCJIOXKHAIOT IIpouecc H3rOTOBIeHHA KOHcTpyKuMH. IIpoOmemy peliaeT IpesIO%KeHHEIM B JaHHOM CTaTbe clOcoO TpexMep- 
HOrO apMupoBaHna — dbeuITHHT. ITO JOKabHOe apMHpOBaHHe KOMMO3HTa B HallpaBIeHHH Z IIpH MHHMMAaJIbHBIX MIpOH3- 
BOJICTBCHHBIX H3MeHeHHAX. CTerleHb Z-apMupoOBaHHA OMpeeAeTCA TWIOTHOCTHIO (pesITHHTA, T.e. KOJIMYCCTBOM y{apoB 
urbI Ha | cM? rKanesoro maxkerta. [ens paSoTbl — OLeHHTh BIIMAHHe (eITHHTAa Ha MEXKCIIOMHY!O TpelIMHOCTOMKOCTS KOM- 
MO3HTHOrO MaTepnalia. 

Mamepuajei u memoooi. MexcnoiHy0 BaA3KOCTb pa3pyleHHa Gy. OlIpeeMAIM Ha CTCKIIOTKAHM MOJIOTHAHOYO Mepe4e- 
TeHua c (pertunrom 10 cm’. Marepuam mponmtTsisaru cmMomoi Sran-370 u orBepaqurenem Sran-45. SkcrepuMeHrTs! 110 
cranfaptam ASTM D7905M-—14 u OCT 33685-2015 nposoyumm Ha HcnbItTaTesIbHOM Mamuue Instron 5900R. Hanpsxen- 
HOe COCTOAHHE Y BEPLIMHbI TPeLMHbI aHaIM3HPOBAIM C MO3HI[MM HeOKaJIbHOM TeOpHu NpOuHOCTH B IporpamMe Ansys 
Workbench (onmma «cTaTH4ecKHi MpOUHOCTHOM aHas3»>). 3a\eHCTBOBAaJIM MeTOJ, KOHCYHBIX 9IEMeHTOB (MK3). 
Pezynomamoi ucciedosanua. [Ina oOpa3yoB paccMOTpeM KpuBble «Harpy3Ka — TepeMeljeHHe». Baranciuiu 3Ha4e- 
Hua Gye. O60OujuIu MToru ENF-ncnpiranuii W14 WOTHOCTH (berTuHTa 0 cm 4 10 cm. CpaBHusIM KOHTpOUbHEIe OOpasiib 
WM OOpa3ubI c (berTHHToM. B mocueyqHeM ciy4ae Gre OKa3aacb BbILIe Ha ~33 %. PaccuuTaM HallpsKeHHOe COCTOAHHE y 
BePLIMHBI TpewHHb! mpx DCB- u ENF-varpyxennu. Busyarm3npoBasin B BAe rpa:bHKOB H I{BeTOBbIX {MarpaMM 3aBHCH- 
MOCTH MaKCHMA€JIbHBIX HOPMA@JIBHBIX M KaCaTeJIbHBIX HallpsKeHHH, a TakKoKe NepememjeHni. JI nowy4aeHuad pacueTHBIX 3a- 
BHCHMOCTeli «Harpy3Ka — TepemeljeHve» c TomMonbIo MKD ucnomb30BaiM OOpaTHbIM MeTOA MouyYeHuA TpaHcBep- 
CaJIbHO-CA/BHTOBBIX KOHCTaHT. Harpy2keHue 110 cxeme DCB noxa3aui0, 4TO (PeITHHT MO3BOJACT YBEIMYHT pees IpOuHo- 
CTH Ha pacTsKeHue B HallpaBsieHun Z Ha ~18 %, c 39 Wo 46 MIa, a B mIocKoctu XZ — Ha ~16 %, c 77 MI a no 89 MIa. 
O@cystcoenue u 3aknio4“enue. DeNTHHT Kak CIOCOO JOKaIbHOrO TPeXMePHOrO apMUPOBaHHA YCHIMBaeT Me2KCIIOMHYIO 
TPeLWIHHOCTOMHOCTh MOJIMMepHbIX TKAHEBbIX KOMIO3HTOB. OH NO3BOJIAeT COKPATHTb WIOWayb paccioeHu Mocse J0- 
KaJIbBHbIX Y1apoOB IIpH SKCIUIyaTaljuu KOHCTpyKuM. [uOKad TeXHONOrMA (PeNTHHTAa TaeT BOSMOXKHOCTL CO3LaBaTb 30HbI 
C MIpOH3BOJbHOM MWJIOTHOCTEIO YapoB, MOBbIMMAA TPeLIMHOCTOMKOCTb JIMIIb B HeOOXOAMMbIX MecTaxX KOHCTpyKUHH. 
MK®9-aHasIM3 HallpsoKeHHOTO COCTOAHHA Y BEPLIMHbI Tpel{HHbI B paMKaX HEOKAIbHOM TEOPHU NPOUHOCTH NOKa3all, YTO 
B IIPOUHOCTHBIX pacueTax TpelHHY PaccIOCHHA MOXKHO paCCMaTPUBATb KaK KOHICHTpaTOp HallpsKeHHH. 


Kunrouesble CJ10Ba: apMUpoBaHHe MOJIMMeCpHBIX TKAHCBEIX KOMHO3HTOB, TpaHCBepCaJIBHO-CABHTOBaxA TIPOUHOCTh, MexKcowHad 
TPeLIMHOCTOMKOCTS, MO?*KCJIOCBAA BAZSKOCTb pa3pyuieHHaA, (pesITHHTOBOe JIOKaJIBHOe TpexXMepHoe apMHpoBaHHe 
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Introduction. Fibrous polymer composites are widely used, in particular, in aviation and space engineering due to 
their significant rigidity and strength in the fiber orientation (plane XY) [1]. However, the transverse shear strength of 
these materials is quite low [2], as it is determined by the features of the polymer matrix [3]. Reinforcement of polymer 
fabric composites in the Z direction provides reducing the sensitivity to delamination, i.e., increasing the viscosity of 
interlayer fracture [4]. 

Various methods of three-dimensional reinforcement of polymer fabric composites are known [5]. However, they 
create additional difficulties in the manufacture of structures made of polymer fabric composites [6]. The method of three- 
dimensional reinforcement proposed in this work, felting [7], makes it possible to obtain a locally reinforced composite 
in the Z direction with minimal changes in the production process. The degree of Z-reinforcement is determined by the 
felting density, i.e., the number of needle punches per | cm? of the area of the fabric package [8]. 

The use of various methods for determining the fracture toughness of polymer composite materials [9] makes it 
possible to conduct studies on various samples [10] and with different loading methods [11]. One of the most common 
approaches is the three-point bending method. In this case, a beam-shaped delamination sample is used. We are talking 
about ENF tests (End-Notched Flexure — bending of a sample with edge delamination) [12], which involve transverse 
shear loading. This makes it possible to determine the interlayer fracture viscosity of Ge — mode II fracture. High shear 
stresses occur at the crack tip [13]. 

Another common way to determine transverse characteristics is the Double Cantilever Beam method (DCB tests). In 
DCB tests, the value of the interlaminar fracture toughness Gj). is determined under separation loading — fracture 
according to mode I [14]. The delamination crack spreads due to the action of normal stresses [15]. 

The presented work was aimed at the evaluation of the effect of felting on the interlayer fracture toughness of a 
composite material. To do this, ENF tests (bending of a sample with an edge separation) of a composite material with 
increased crack resistance due to the use of felting were carried out. Previously, the authors studied the effect of felting 
on the interlayer crack resistance of a composite material during DCB tests [16]. A computational model based on the 
nonlocal theory of strength has been developed. It provides for the calculation of the stresses that occur in ENF and DCB 
samples, for cracks of various lengths, using the finite element method (FEM). 

Materials and Methods 

Experimental determination of interlayer crack resistance by the ENF method. The samples were made of cross- 
woven fiberglass with a layer thickness of 0.2 mm. A package of dry two-layered fiberglass was punched on a felting 
machine with a felting density of 10 cm (10 punches with a felting machine needle per 1 cm? of dry glass fabric). The 
fiberglass package was punched in such a way that after impregnation and hardening, the initial crack did not fall on the 
felting area. To create an initial crack between two layers of fiberglass, an aluminum foil with a thickness of 11 um, 
coated with a Vs-M parting lubricant, was placed. Fiberglass was impregnated with resin Etal-370 and hardener Etal-45. 
For the production of reference samples, two layers of dry fiberglass were impregnated with Etal-370 resin and Etal-45 
hardener without punching on a felting machine. After impregnation, plates from glass textolite STEF (electrotechnical 
glass-cloth-base laminate) were glued to two layers of fiberglass (Fig. 1 a). Samples with a length of 150 mm and a width 
of 16 mm were obtained by cutting the hardened plates using a high-speed circular saw. 


Fiberglass STEF 2 mm 


Felting 
area 


a) 


Fig. 1. Configuration and parameters of ENF tests according to mode J/: 
a — loading scheme; b — photo of tests 
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The Instron 5900R test machine with a loading speed of 10 mm/min was used. The distance between the supports was 
2L = 100 mm. The initial crack length for all samples was ap = 25 mm. To perform compliance calibration over a wide 
range of crack lengths, one sample of each type was unloaded and reloaded. The resulting crack served as the initial crack 
in the next loading cycle. 

Calibration recommended by ASTM D7905M-14! and GOST 33685-2015? standards was used to process the test 
results. This approach provided determining parameters A and m for each felting sample and each non-felting control 
sample from the linear dependence of two quantities — compliance of sample C and cube of the crack length a?: 


C(8/P(8))=A+m-a’, (1) 
where P — load applied to the sample; 6 — displacement. 
After calibration and determination of parameters A and m, the crack length can be found from expression (1): 


o-(o=4)". re 


m 


The moment of the delamination onset is determined by the condition C(6) = C(Pmmax). Value of the interlayer fracture 
toughness at the separation onset (crack development): 


3m: Po a? 
Gite = —— (3) 
where Pmax — maximum load; a — crack length calculated by formula (2) at load Pina; b — sample width. 

Calculation of the stress state at the crack tip under loading according to the DCB and ENF schemes. The stress 
state of a crack-like concentrator is estimated from the perspective of approaches that use nonlocal stresses [17], i.e., 
averaged on some basis [18]. The calculation model also includes the assumption of linear-elastic behavior of the material 
up to destruction [16]. 

The main hypothesis is that the strength criterion of the composite, which includes all components of stress averaged 


on the basis A, is responsible for the development of a crack-like concentrator (Fig. 2): 


2 2 2 
max O Max O >» Max T MAX O ,, * Max oO 
2] { 2) +( =| zh XA <1, (4) 


Z, Xx, Ss Z,°X;, 
where Z, and _X;— rupture strength in the Z and_X directions; S — shear strength in the plane XZ. 


— 
{ Max On 4 MAX Ox, 
a) b) 


Fig. 2. Stresses averaged on the basis A at the crack tip: 
a— ENF tests; b — DCB tests 


' ASTM D7905/D7905M-14. Standard Test Method for Determination of the Mode II Interlaminar Fracture Toughness of Unidirectional Fiber-Rein- 
forced Polymer Matrix Composites. URL: https://cdn.standards.iteh.ai/samples/89096/03be6b5e53664f13a8703bb4342d98 la/ASTM-D7905- 
D7905M-14.pdf (accessed: 22.04.2024). 

2 GOST 33685-2015. Polymer Composites. Test Method for Determination of the Interlaminar Fracture Toughness under Shear. (In Russ.) URL: 
https://docs.cntd.ru/document/1200127774 (accessed: 22.04.2024). 
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Due to the presence of symmetry planes, a three-dimensional 1/2 crack sample model was constructed for ENF 
loading (Fig. 3), and 1/4 crack sample — for DCB tests (Fig. 4). Calculations were performed in the ANSYS Workbench 
program (option “static structural’’). 


Contact 
Frictionless 


Contact 
Bonded 


Fig. 3. Finite element model 1/2 of the sample and a fragment of the grid for ENF tests 


0.0 30.0 60.0 (mm) 


= SC 


15.0 45.0 


Fig. 4. Finite element model 1/4 of the sample and a fragment of the grid for DCB tests 


When creating the finite element grid, parameter 1, = 0.75 mm in width of the sample [19] and parameter 1, = 0.2 mm 
in thickness of the sample were set, which corresponded to the thickness of the modified layer [20]. One finite element in 
the thickness of the layer was set in accordance with the layer wise theory used in assessing the strength of layers within 
the framework of mesomechanics of composites [21]. In ENF tests, the total displacements in the sample were much 
greater than the local displacements from the loading roller (Fig. 1 5); therefore, the grid of finite elements was not 
condensed at the places of application of loads and supports (Fig. 3). Properties of fiberglass used in the calculation: 

- elastic modules E,. = Ey = 23 GPa, E, = 9 GPa; 

- shear modules Gy, = Gy: = Gy = 6000 GPa; 

- Poisson's coefficients [x = 0.15, Uy: = Hx: = 0.3 [22]. 

Since the volume fraction of transverse reinforcement is less than 1% [16], it is assumed in the calculations that the 
elastic properties of fiberglass do not change under felting. 

Dependence P(6) was calculated in accordance with the sequence described below. 

1. FEM-calculation of the maximum stresses maxoy, maxoy and maxtx, and displacements of point 6 of application 
of load P= 1H for cracks with the given lengths in the range a = 20...90 mm (DCB) and a = 25...40 mm (ENF) was 
performed. 

1. Approximation dependences o. = fla, P) = P:by-a; 6, = fla, P) = P-bra; %a.= fla, P) = P-b3:a; 5 = fla, P) = P-cy-a*-(DCB) 
and oa = fla, P) = P-by-a; 6n.= fla, P) = P-bra; twa = fla, P) = P-(bs-a + ds); 8 = fla, P) = P(cv-a? + cx.a® + c3:a + c4)(ENF). 

2. were constructed using the least squares method. 


3. Load P.,(ao) and displacement 6,,, at which the initial crack length ao would increase abruptly by Ax = 0.75 mm 
upon the violation of the strength criterion, were determined (4). 
4. With crack length agtnA, loads P(aotnA) and displacements 6(n) for n>0 were determined. 
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Research Results 

Results of the experimental determination of interlayer fracture toughness by the ENF method. Figure 5 shows 
the “load — displacement” curves for all tested samples. All “load — displacement” curves have an area with constant 
compliance (C;;,), corresponding to the linear “load — displacement” ratio. Values Cj, were used for calibration. 


P,kKN 


Fig. 5. “Load — displacement” diagrams of ENF tests 


without felting; = = = =— with felting 


Figure 6 shows the calibration curves. For felting and non-felting samples, compliance is proportional to the cube of 
the crack length. 


C, mm/N 
0.012 ---O 
y=8.44- 10°. + 6.69: 10° el _*2 
0.010 ae pt ee 
pot a 
0.008 =  -8- 
goer? y=6.98-10°x+6.29-10° 
0.006 
0.004 
0.002 
0.000 
0 10,000 20,000 30,000 40,000 50,000 60,000 a, mm? 


Fig. 6. Dependence of compliance of sample C on the cube of crack length a?: 
o — without felting; e — with felting 


To calculate values of the crack length a*, corresponding to the compliance at the delamination onset C(Pinax), the 
obtained calibration curves and equation (2) were used. For the found values of the crack length a*, values Gy. were 
calculated using equation (3). The results are shown in Table 1. 
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Table 1 
ENF Test Results 

ee ao, mm | Cin, mm/N ear a*,mm | Pax, N Gite, kI/m? ae ao 
0 25 8.032 8.845 29.4 504.6 1.723 
0 25 7.980 8.807 29.3 544.6 1.983 
0 25 8.299 9.392 31.7 489.4 1.908 pe One 
0 25 7.905 8.900 29.7 505.2 1.746 
10 25 7.587 8.337 30.8 627.3 2.432 
10 25 7.849 8.677 32.4 625.8 2.682 
10 25 7.937 8.811 33.0 578.4 2.376 2.441 + 0.154 (6.3%) 
10 25 7.824 8.594 32.1 581.4 2.261 
10 25 7.880 8.818 33.1 589.3 2.456 

* variation coefficient 


Felting samples showed a significant (by ~33%) increase in the interlayer fracture toughness Gy. compared to the 
control samples. After testing, felting samples were separated with a sharp knife and examined under a microscope. 
Micrographs of the zone without felting (area of the initial crack) and the zone with felting (area of crack development) 
are shown in Figure 7. When cracks develop, the fibers elongated during felting are destroyed, because their length is 
greater than the critical one [16]. 


Fig. 7. Micrographs of felting samples after ENF tests: 
a — zone without felting (area of initial crack); 
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b — zone with felting (area of crack development); 
c — zone with felting (enlarged scale) 
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Results of the calculation of the stress state at the crack tip under loading according to the DCB and ENF 
schemes. Figures 8—9 show the dependences of stress maxoa(a), maxow(a), maxtxx(a) and displacements 8(a). 
Conditions: P = 1 N, loading according to DCB and ENF schemes. 


max 6,, MPa max 6,, MPa 1 
6.0 y=2.07- 10x 
1.6 
5.0 R’=0.99 
4.0 1.2 
3.0 08 
2.0 
1.0 0.4 
0.0 0.0 
0 20 40 60 80 a,mm 0 20 40 60 80 a,mm 
a) b) 
MAX Ta, MPa 6,mm 
0.6 1.2 
0.4 0.8 
0.2 0.4 
0.0 0.0 
0 20 40 60 80 a,mm 0 20 40 60 80 a,mm 
c) d) 


Fig. 8. DCB-loading. Dependences of maximum stresses and displacements of the crack length at P= 1 N: 
a— dependence of normal stresses maxo.a(a); b — dependence of normal stresses maxoa(a); 
c — dependence of shear stresses maxtxa(a); d — dependence of displacements 5(a) 


max 6,, MPa >, —_ —___,_—__"——_ max 6;, MPa 
y=4.13-10°x y=4.13- 10x 
0.20 0.20 
R’ =0.999 R’ =0.999 
0.15 0.15 


0.00 + 0.00 + 
0 10 20 30 40 a,mm 0 10 20 30 40 a,mm 
a) b) 
Max T%,, MPa A 6, mm 1 ch, 5 : r 
y=-9.40-10 x + 1.63 -10°x — 4.63 - 10° x+ 
0.20 
+ 1.09- 10° 
0.15 0.12 ; } 
2 R= ae 
8 0.10 0.08 | —, 
=| 
3 0.05 0.04 | 
pI 
& 0.00 4 0.00 
g 0 10 20 30 40 a,mm 0 20 40 60 80 a,mm 
— 
Pa c) d) 
b= Fig. 9. ENF- loading. Dependences of maximum stresses and displacements of the crack length at P= 1 N: 


a — dependence of normal stresses maxo.(a); b — dependence of normal stresses maxoz(a); 
c — dependence of shear stresses maxtxa(a); d — dependence of displacements 8(a) 
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Examples of stress distribution at the crack tip are shown in Figures 10—11 with crack length ap = 30 mm. 
aS a 


2.108 1.325 0.361 0.512 —1.822 MPa 0.639 0.450 0.262 0.073 -0.115 MPa 
a) b) 


aS Eee 
0.364 0.234 0.104 -0.026 -0.156 MPa 


¢) 


Fig. 10. Stresses at the crack tip under DCB loading: a — normal stresses oy; 
b — normal stresses oz; c — shear stresses Txza 


0.649 0.361 0.072 —0.216 -0.505 MPa 0.053 0.292 0.006 —0.018 —0.041 MPa 


a) b) 
a | 


0.155 0.112 0.069 —0.026 —0.017 MPa 
°) 


Fig. 11. Stresses at the crack tip under ENF loading: a — normal stresses ox; 
b — normal stresses 6; c — shear stresses Tx 


To obtain the calculated “load — displacement” dependences using FEM, the strength characteristics of the composite 
in the main directions are taken into account, i.e., criterion (4). Direct obtaining of transverse shear constants involves a 
certain difficulty; therefore, the reverse method is considered below. With this approach, the constants vary, and their 
best combination is found. This means that the calculated and experimental loading diagrams are in good agreement 
(mean square deviation of displacements at specified loads is minimal). 

The results of the calculation under loading according to the DCB scheme were compared to the authors' experiment, 
which was considered in [16]. Samples were made in the same way. The tests were carried out in accordance with GOST 
R 56815-20153 and ASTMD5528-14¢ standards. 

The calculation was performed for loading according to the DCB scheme. When calculating dependence P(S) for 
samples without felting, the following stress limits were found and rounded to integer values: Z;= 39 MPa, X; = 360 MPa, 
S = 82 MPa. The obtained values 360 MPa and 39 MPa correspond to the data on the strength of fiberglass specified 
in [23]. For felting samples (density 10 cm”), the calculated values were X;* = 270 MPa, Z;* =46 MPa and S* = 97 MPa. 
Thus, the use of felting made it possible to increase the rupture strength in the Z direction from 39 to 46 MPa (by ~18%). 


3 GOST R 56815-2015. Polymer Composites. Method for Determination Specific Work of Exfoliation in Tearing Off Conditions. (In Russ.) URL: 


https://docs.cntd.ru/document/1200131393/titles (accessed: 22.04.2024). 
4 ASTMD5528M-21. Standard Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Compo- 


sites. https://doi.org/10.1520/D5528 D5528M-21 


Mechanics 


223 


https://vestnik-donstu.ru 


224 


Forental GA, et al. Increasing the Interlayer Fracture Toughness of Polymer Fabric Composites Using Local 3D-Reinforcement (Felting) 


When loading according to the DCB scheme, the shear strength limits in the plane XZ S and S* do not make a big 
contribution to criterion (4); therefore, values S = 82 MPa and S* =97 MPa obtained in calculations according to the 
DCB scheme need to be clarified according to the ENF loading scheme. Note that the effect of normal stresses in the X 
and Z directions is insignificant compared to shear stresses under loading according to the ENF scheme. Therefore, in the 
calculations, when searching for the values S and S*, values Z,, X;, Z;* and_X;* were taken from solving the inverse problem 
under loading according to the DCB scheme. 

Values S = 77 MPa (without felting) and S* = 89 MPa (with felting) were determined from the condition of the best 
consistency of the experimental and calculated curves P(5) (mean square deviation of displacements at given loads is 
minimal). Evidently, felting made it possible to increase the shear strength in the plane XZ by ~16%. 

Figures 12-13 show the experimental “load — displacement” diagrams, as well as calculated dependences P(6) for 
the found values of stress limits: 

- Z,;= 39 MPa, X; = 360 MPa, S= 77 MPa (for control samples without felting); 

— Z;* = 46 MPa, X;* = 270 MPa, S* = 89 MPa (for samples with felting). 


Fig. 12. Experimental diagrams “load — displacement” [16] and calculated dependences P(8): 

— DCB tests of samples without felting; — calculation without felting (Z;= 39 MPa, X;= 360 MPa, S = 77 MPa); 
- -- DCB tests of samples with felting; - - - calculation with felting (Z;* = 46 MPa, X;* = 270 MPa, S* = 89 MPa) 
PLN 

0.6 

0.5 

0.4 
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6, mm 


Fig. 13. Experimental diagrams “load — displacement” and calculated dependences P(5): == ENF tests of samples without felting; 
— calculation without felting (Z;= 39 MPa, X;= 360 MPa, S= 77 MPa); - - - ENF tests of samples with felting; 
--- calculation with felting (Z;* = 46 MPa, X;* = 270 MPa, S* = 89 MPa) 


Discussion and Conclusion. Studies of the fabric composite have shown that felting with a density of 10 cm~ 
increases the viscosity of the interlayer fracture Gye by ~33%. 

Using FEM, the stress state was analyzed in a quasi-static elastic formulation of the problem and with a nonlocal strength 
theory for the developed numerical models of a beam with cracks of the known length. The distinctive feature of the 
calculations was that they did not use contact algorithms, but only considered the destruction of the composite layer closest 
to the crack, and the corresponding change in the area of gluing the layers. That is, the crack was considered as a stress 
concentrator. The composite strength criterion, which contained three parameters and was recorded through averaged 
stresses, provided using the method of step-by-step crack advancement to predict the “load — displacement” curve. 

The use of felting with a density of 10 cm™ increases the stress limit of the composite in the Z direction by ~18%, and 
the shear strength in the plane XZ — by ~16%. This became known from solving the inverse problem, i.e., searching for 
the strength characteristics of the material according to criterion (4) and the “load — displacement” curve. 

The results of the presented research will find their practical application. They can be specifically used in problems 
of forecasting defects, such as delamination (e.g., in low-speed impacts on composites in aircraft skin). The research 
results will be useful for eliminating these defects with the help of felting. 
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Abstract 

Introduction. When studying composite materials for construction purposes, it is needed to consider the mechanisms of 
formation of the structure and properties of modern concretes in the process of strength development. In studies of modern 
composite materials based on cement binder, there is no information about the development of structural defects and 
destruction of the material at the initial stages of strength development. This information can be obtained using X-ray 
computed tomography, a promising method of nondestructive testing of the state of the material. Therefore, the objective 
of this work was to study the formation and propagation of cracks in samples of fine-grained concrete with different 
fractional composition of sand due to natural processes of cement shrinkage, as well as the mechanics of destruction of 
samples of modified fine-grained concrete when applying a compressive load at the early stages of strength development. 
Materials and Methods. The study used fine-grained concrete mixtures of three compositions with different sand 
gradation. The tomography samples were made by placing fresh mixtures in polymer cylindrical containers. Tomography 
of the samples immediately after manufacture, as well as after 8 and 51 days, was performed in a YXLON Cheetah 
microfocus X-ray machine. The composition with two-fraction sand was modified by mechanical activation of the 
components, 20x20x20 mm cube samples were made. Further, compression tests were performed at the Instron 
installation after 3 and 7.5 hours, and then — tomography of the destroyed samples. 

Results. It was established that the destruction of contact zones depended on the ratio of the size of the fractions. In the 
presence of a bulk of coarse sand grains in concrete, the destruction of contact zones was more pronounced and had a 
main mode. When using fine or polyfraction sand, contact zones were destroyed locally and had a visually smaller area. 
The images of the destroyed modified sample, tested 3 hours after manufacturing, showed clear cracks and indents on the 
edges, which indicated the elastic-plastic nature of the destruction. In 7.5 hours, the edges of the sample upon destruction 
were covered with a network of small cracks; inside the sample there were also numerous cracks and microcracks, which 
indicated brittle fracture. Based on the obtained images of the deformed structure of modified concrete, the mechanism 
of transition from elastic-plastic destruction of the material to brittle one was clearly visible. 

Discussion and Conclusion. The studied dependences of the influence of the size of fine aggregate on the mechanisms 
of formation and propagation of structural defects contribute to the theory of the processes of destruction of fine-grained 
concretes. The results obtained prove the prospects of using X-ray computed tomography as a method of nondestructive 
testing of the internal structure of fine-grained concrete, including at the early stages of strength development. 


Keywords: X-ray computed tomography, fine-grained concrete, deformation, elastic-plastic failure, brittle fracture 
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Opueunasbnoe aunupuyeckoe uccaedoganue 


UccueqoBanne nponeccoB 3BO.TIONHH Je@eKTOB CTIPyKTYPbI MeJIKO3ePHUCTHIX OeTOHOB 
MeTOJAaMH KOMIIbIOTepHOH TOMOTpadun 


A.B. Hysarona® > M.A. JImurpuera®, A.O. Tosnnnen®, B.H. Jleinun® 
Banrulickuit dbexqeparbubiit yHuBepcuteT uMeHu MMmanyuna Kauta, r. Kanmuuurpag, Pocculicxas Dexepanna 
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AHHOTalna 

Beeoenue. pu u3sysenuu KOMMO3HIMOHHBIX MaTepHasIOB CTPOHTeIbHOTO Ha3HadeHHA AKTYAJIbHbIMH ABJIAIOTCA HCCJIe- 
JOBaHHA MeXaHH3MOB (POPMHPOBaHHA CTPYKTYPbI MH CBOMCTB COBPeMeHHBIX OeTOHOB B TIpowecce Habopa mpouHoctTu. B 
MCCI€HOBaHHAX COBPCMCHHBIX KOMIO3HIMOHHBIX MaTepHasIOB Ha ICMCHTHOM BKYLIEM OTCYTCTBYIOT CBEeCHHA O pa3- 
BHTHH JedeKTOB CTPyKTypbl HM pa3pyWIeHHH MaTepHasia Ha HauyaJIbHBIX cpoKax HaOopa MpounocTu. Takve cBeqeHHAa 
MO2%KHO TIOJIYIHTb C MOMOLIbIO PeEHTTCHOBCKOM KOMMbIOTepHO TOMOrpaPun — MepcieKTHBHOrO MeTOsa Hepa3pyliaro- 
Wero KOHTpOsIA COCTOAHHA MaTepnana. IlosTomy WesIbIO JaHHOM paoorTsl ABHJIOCb H3y4eHHe OOpa3z0BaHHA HM paciipoctpa- 
HeHHA TpelHH B OOpa3ljax MeIKO3ePpHUCTOTO OeTOHA C pa3JIMYHbIM ()pakKUMOHHbIM COCTaBOM IlecKa BCJICACTBHe ecTe- 
CTBeHHBIX IIpoleccoB ycayKH W€MCHTHOTO KaMHA, a TaKXKe MEXaHHKH pa3spylieHHa OOpa3slOB MOAMPUUMpOBaHHOTO MelI- 
KO3e€PpHUCTOrO OeTOHa IpH NPHIOKeHHU CKUMAIOUIeH Harpy3KH Ha paHHUX cpoKax HaOopa MpoyHOcTH. 

Mamepuaavi u memoooi. B uccneqoBaHvi MCHOb30BasIMCb MCJIKO3CPHHCTHIe O€ETOHHBIe CMeCH Tpex KOMMO3HUHH Cc 
pa3IM4HBIM rpaHyJIOMeTpH4ecKHM COCTaBOM TiecKa. OOpa3ubl Ia TOMOrpaduu OBLIN U3rOTOBJICHbI IyTeM MOMeLe- 
HUA CBEX%KHX CMeCel B NOJIMMepHble WHIMHApH4eckune KOHTelMHepbl. Tomorpadusa oOpa3yoB cpa3y Wocue H3roTOBIJIe- 
HHA, a Take Yepe3 8 u 51 cyTKH NpoBOAMIAaCh B MUKpOdoKycHOl peHTTeHOBCKON ycTaHoBKe YXLON Cheetah. Co- 
cTaB C ByX(pakKUMOHHBIM IeCKOM ObluI MOAMPULMpOBaH MexaHHyecKON aKTHBaLMel KOMIOHECHTOB, H3rOTOBJICHbI 
oOpa3ubI-KyOuKu 20X20*20 mo. Jlanee Ha yctranoske Instron mpoBeyeHbI MCIbITAaHHA Ha cxKaTHE Yepe3 3 u 7,5 4acoB 
MW 3aTe€M — TOMOrpadua pa3spyWIeHHBIxX OOpa3si{oB. 

Pe3yibmamot uccnedosanua. YcTaHOBNeHO, ITO paspylIeHHe KOHTAKTHBIX 30H 3ABHCHT OT OTHOMICHHA pasMepos (ppak- 
qui. B mpucyTcTBHH OoNbUIOTO KOUM4eCTBAa KPyMHbIX YaCTHL MWecka B Tele OeTOHa pa3pylieHHe KOHTAKTHBIX 30H Ooslee 
BbIPaxKeHO HM HMCCT MarHCTpasibHbIt xapaktep. [pv ucnoub30BaHHM MeJIKOLO HIM MOUMPpakUMOHHOTO TecKa KOHTAKT- 
HbI€ 30HBI pa3pyWlaloTcA JIOKaJIbHO HM MMEIOT BH3YaIbHO MeHbITyr0 Tomjayb. Ha w300paxKeHHAX pa3pyleHHOro MOJH- 
(uUMpoBaHHoOro OOpa3la, HCIIbITAHHOrO Yepe3 TPH Yaca MOC H3TOTOBJICHHA, MpOCIeKUBAaIOTCA YCTKMe TpelHHbI 
BBIKOJIbI Ha TpPaHAX, TO TOBOPHT OO yupyro-acTH4ecKOM xapaktTepe pa3pyuienua. Uepe3 7,5 uacos rpaHu oOpa3ia pu 
pa3pyLWIeHHH MOKpbIBaIOTCA CCTKOM MeJIKHX TPeLIHH, BHYTpH OOpa3ila TakxKe OOpa3syeTCA MHOXKECTBO TpeLHH H MUKpo- 
TPeLIMH, TO CBUCTEMbCTBYeT O XPyHKOM pa3pyweHuH. Ilo nomyyeHHbIM H300paKeHHAM JePOpMUpoBaHHON CTpyK- 
TYpbl MOAM@UUMpOBaHHOrO OeTOHA HarATHO POCMexKUBAaeTCA MCXaHH3M TlepexOa OT yIIpyro-IWlacTH4ecKOro pa3py- 
WeHHA MaTepHasia K XPylKOMy. 

O6cystcoenue u 3axmo4yenue. Visy4enHble 3aBMCHMOCTH BIMAHHA Pa3MePOB MCJIKOFO 3allOJIHHTesIA Ha MCXAaHH3MBI 00- 
pa30BaHHa HW paciipocTpaHeHua JeteKTOB CTPyYKTYPbI BHOCAT BKIa B TEOPHio IPOMeccoB paspyWIeHHA MeJIKO3CPHHCTHIX 
Oberouos. Ilomy4eHHble pe3yIbTaTbl JOKa3bIBalOT IePCMeKTHBHOCTb MIpHMeHeHHA PCHTTCHOBCKOM KOMIIBIOTepHOH TOMO- 
rpaun Kak MeTOjla Hepa3pyllalolero KOHTPOJIA BHYTPeHHel CTpyKTypbI MeJIKO3epHUCcTOrO OeTOHAa, B TOM 4NCIIe Ha 
paHHux cpokax Hadopa MpouHocTH. 


Ks04uesbie C10Ba: peHTTeHOBCKad KOMIIbIOTepHat TOMOTpaHA, MeJIKO3epHUCTBIM OeTOH, JePopMallHA, yipyro-mlacTH4eckoe 
pa3pylleHue, xpylkoe pa3spylleHHe 


Baarojapuoctu. ABtTopbl BbIpaxaroT OsarogapHocts Lnasesoti Mapuu Bragumuposue, Korat Anune J[Imutpuesue 
3a TOMO B MIpHrOTOBJICHHH OOpa3slOB KOMMO3MIMOHHBIX MaTepHasIOB JIA MpOBeAeHHA HCCIeOBaHH CTpyKTypbI Ha 
PCHTTeHOBCKOM KOMIIbIOTepHOM Tomorpade. 


Asia waTapopanna. ITy3arosa A.B., Mutpvesa M.A., Tosrmuer A.O., JletiiH B.H. Uccneqospanne mpoleccos 9BOONMU 
TeeKTOB CIPyKTYPbI MeIKO3CPHUCTHIX OCTOHOB MeTO{aMU KOMIbIOTepHoH ToMorpadun. Advanced Engineering Research 
(Rostov-on-Don). 2024;24(3):227—237. https://doi.org/10.23947/2687-1653-2024-24-3-227-237 


Introduction. X-ray computed tomography is a promising technique of nondestructive testing of the state of the 
material. In the concrete industry, tomography is suitable for determining the structure of concrete samples, microcracks, 
internal fractures, and studying the distribution of pores and aggregate particles. X-ray computed tomography provides 
the construction of a model of the microstructure of cement paste, allows us to study the development of cement hydration 
processes [1, 2], to make forecasts of the formation of mechanical characteristics and fracture conditions [3, 4]. Computed 
tomography is actively used for studying the average density and porosity of high-strength lightweight concretes [5], the 
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morphology of structural heat-insulating concrete [6], building a mesoscale 3D model of foamed concrete [7], 
investigating the formation and distribution of pores in lightweight concrete [8], analyzing the microstructural 
characteristics of concrete samples with various aggregates [9, 10], developing three-dimensional mesoscale models for 
constructing a finite element grid in modeling the structure of concrete [11, 12]. In comparison to standard 2D X-ray 
methods, the construction of 3D models of multicomponent concrete samples is promising for the study of the 
fundamental mechanisms of formation of the structure and properties of modern concretes. 

The most vulnerable area of fine-grained concrete under loading is the contact zone — the contact areas of cement 
stone and aggregates. When exposed to external loads, it is from these areas that the formation of microdefects and 
microcracks starts, whose development causes the defect formation at the macro level, which can lead to loss of bearing 
capacity and structural failure. The destruction of the contact zones occurs due to the difference in characteristics of the 
bordering components (Young's modulus, Poisson's ratio, thermal linear expansion coefficient, sizes of the contacting 
phases, microdefects on the interface of the phases) [13]. The X-ray computed tomography method is promising for 
studying the evolution of contact zones, including at the early stages of hydration. It provides studying the structure 
without destroying the sample directly during the hardening process. Contact zones, as a rule, have higher porosity and 
low strength, as a result of which cracks are formed in these zones [4]. The strength of the contact zones is also affected 
by the size of the aggregate. It has been established that homogeneous and more durable contact zones are formed in 
concretes with combined aggregates (coarse fraction and crushed) [14]. There are numerous modern studies on the 
formation of contact zones of cement stone with reinforcement [15], cracks in popcorn concrete [16], defects at the contact 
boundary of sand-cement mortar with a coarse aggregate [17]. However, along with this, the formation of contact zones 
in fine-grained concretes with different fineness and packing density of sand grains is poorly studied. In the modern 
scientific literature, there is no description of the effect of the size of a fine aggregate on the formation of contact zone 
defects in fine-grained concretes. Thus, the study of the formation of defects in the structure of fine-grained concrete 
containing sand of aggregate fractional make-up, to obtain a visual picture of crack propagation by X-ray computed 
tomography is urgent. 

To reduce the stresses that occur in the contact zones, microfillers, similar in their properties to cement stone, are used. 
These fillers, having an increased specific surface area, create additional contact zones between which the stresses arising 
from the hardening of the binder are redistributed. Hardening of the contact zones can be achieved through introducing 
mechanically activated components into the concrete mixture [13]. Mechanical activation of separate components 
contributes to the formation of a denser structure, the preservation of uniformity in composition, the development of initial 
strength due to the acceleration of the hydration reaction and the growth of crystallohydrates of cement stone, as well as 
the reduction in the setting time [18, 19]. The development of structural defects in fine-grained concretes modified by 
mechanical activation of components has also been poorly studied. Among modern scientific papers, there are very few 
works devoted to the study of crack propagation processes in the modified structure of fine-grained concretes at the initial 
stages of strength development. Therefore, the use of the computed tomography method to study the mechanics of 
destruction of samples of modified fine-grained concrete at the initial stages of hardening is important today. 

This work was aimed at studying the formation and propagation of cracks in samples of fine-grained concrete with 
different fractional make-up of sand due to natural processes of shrinkage of cement stone, as well as the mechanics of 
destruction of samples of modified fine-grained concrete when applying a compressive load at the early stages of strength 
development. 

Materials and Methods. In the research, to determine the defects of the contact zones of fine-grained concretes with 
different sand gradation, three samples of fine-grained concrete mixture of the following compositions were used: 

— composition No. 1: Portland cement Eurocem 500 super; CEM I 42.5 N (“Petersburg Cement” LLC); monofraction 
sand (fraction 0.63—0.315 mm); water. The ratio of the mixture components was 1:2.56:0.67; 

— composition No. 2: Portland cement Eurocem 500 super; CEM I 42.5 N (“Petersburg Cement” LLC); polyfraction 
construction sand according to GOST 8736-2014 with ISO=1.85; water. The ratio of the mixture components 
was 1:2.56:0.67; 

— composition No. 3: Portland cement Eurocem 500 super; CEM I 42.5 N (“Petersburg Cement” LLC); two-fraction 
sand (fraction 2.5—1.25 mm — 80% of the total mass of sand, fraction 0.63—0.315 mm — 20% of the total mass of sand; 
no intermediate fraction 1.25—0.63 mm); highly active metakaolin (white) (manufactured by CG “Sinergo”, Chelyabinsk 
region); microsilica; water. The ratio of the mixture components was 1:1.75:0.43:0.23:0.15:0.67. 

Freshly-mixed compositions of fine-grained concrete mix were placed in a polymer cylindrical container with a 
diameter of 8 mm and a length of about 70 mm. The diameter of the container was determined in accordance with the 
size of the initial components of fine-grained concrete of the previously mentioned compositions based on the conditions 
of representativeness [20], and representativeness of the volume under study. To track changes in the structure under 
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cement hydration and cement stone shrinkage, tomography of the samples was performed immediately after mixing the 
components and after 8 days of hardening. Samples at the age of 51 days were taken as the final result. 

In the study on the evolution of the crack formation mechanism during the destruction of samples under the action of 
an external compressive force, composition No. 3 with two-fraction sand, modified by mechanical activation of the 
cement and sand composition, was used. Mechanical activation of the components (cement and sand) was carried out 
using a high energy ball mill Retsch EMax. The components were crushed at a rotation speed of the apparatus bowls of 
750 rpm for 5 minutes. Cubic samples of 20x20x20 mm in size were made from concrete mortars, to which an external 
compressive load was applied after 3 and 7.5 hours. After applying the load, tomography of the destroyed samples was 
performed to trace the evolution of the fracture mode of the material. 

The structure of the concrete samples was studied using a YXLON Cheetah microfocus X-ray computed tomograph 
with Y. Cheetah configuration. The device specifications are shown in Table 1. 


Table 1 
YXLON Cheetah Tomograph Specifications 

Characteristic Value 
X-ray tube Open-type 
Operating voltage range, kV 25-160 
Operating current range, mA 0.1-1.0 
Maximum tube power, W 64 
Maximum tube power on target, W 15 
Detector tilt angle +70° (segment 140°) 
Magnification (geometric/maximum) 2000x/17500x 
Maximum sample dimensions, mm 800x500 
Time from sample loading to first image, sec <10 
Time for full tomography of the sample, sec 7 
Time for laminography of the sample, sec 20 
Overall dimensions, mm 1650x1400x1850 
Weight, kg 2,200 


From the perspective of micromechanics of composite materials, the assessment of effective characteristics can be 
introduced by sets of properties of a representative volume of the body under study. The test sample for tomography 
should correspond to a representative volume of material, which makes sense of the elementary macrovolume of a 
microinhomogeneous medium. 

During the experiment, after X-rays pass through the sample, a set of flat X-ray images with an uneven distribution 
of grayness in the images is obtained. This is due to the uneven absorption of X-rays by the components of the studied 
material — the presence of pores, defects, dense inclusions, etc. After reconstruction of 3D images of the sample, the 
grayness gradient is inverted relative to conventional X-ray images: the materials that are most transparent to X-rays, e.g., 
pores, correspond to black, and the densest material is white. The grayness density in full-color representation is 
considered in the range 0-255, where 0 corresponds to black, and 255 corresponds to white. This algorithm is used to 
determine the minimum size of features that could be differentiated as structural components of the sample under study: 
porosity, cement mortar, and aggregate in volume. Further, using this grayscale, highlighting certain numerical ranges, it 
is possible to analyze a separate internal structure, distribution of components and porosity [21]. 

The shooting parameters for all the studied samples remained constant: voltage — 85 kW; current — 45 pA; 
approximation — 8.9; scanning angle — 360°. The survey results represent 1,024 consecutive images of the internal 
structure of the samples. Further processing of the resulting array of two-dimensional images took place in the “Volume 
Graphics Studio” program. To improve the visualization of inclusions, editing layered images by brightness and contrast 
levels was performed. The result of the tomography was a 3D model of the sample and its three projections with the 
possibility of studying the internal structure in any section. 
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Research Results. Images of the internal structure of samples of composition No. 1 with monofraction sand (fraction 
0.63—0.315 mm) at the age of 0, 8, and 51 days obtained by X-ray tomography are shown in Figure 1. The darkest areas 
in the images have the lowest density, in this case, the pores. The hardest particles correspond to the lightest areas. 


Sas y 


2mm 2mm a 2 mM 


Fig. 1. Internal structure of sample No. 1 (monofraction sand): a, d— immediately after preparation; 
b, f— at the age of 8 days; c, f— at the age of 51 days 


The images of the internal structure of the sample of composition No. 1 with monofraction sand showed no changes 
in the contact zones immediately after mixing. By the 8th and 51st days of hardening, dark stripes were visible around 
separate sand grains (indicated by red arrows), corresponding to voids that were formed due to shrinkage of the cement 
stone. Moreover, with the increasing age of concrete, such voids around the sand grains visually became more numerous. 
Enlarged images of separate fractured voids are shown in Figure 2. 
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Fig. 2. Enlarged fragments of Figure 1: a — fragment 1; b — fragment 2; c — fragment 3; d — fragment 4 
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In fragments 3 and 4 corresponding to sample No. | at the age of 51 days, the mode of the manifestation of defects in 
the contact zones around sand particles was most pronounced. This was confirmed by the fact that by 51 days, the process 
of shrinkage of cement stone was almost complete, whereas at the age of 8 days, the shrinkage was in the active phase. 

Images of the internal structure of samples of composition No. 2 with polyfraction sand at the age of 0, 8 and 51 days 
are shown in Figure 3. Separate enlarged fragments are shown in Figure 4. 


2mm 


Fig. 3. Internal structure of sample No. 2 (polyfraction sand): a, d — immediately after preparation; 
b, e —at the age of 8 days; c, f— at the age of 51 days 


Fig. 4. Enlarged fragments of Figure 3: 
a — fragment 1; b — fragment 2; c — fragment 3 
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The development of defects in the contact zones of samples of composition No. 2 with a polyfraction aggregate, as in 
composition No. | with a monofraction aggregate, manifested itself by the 8th day of hardening, the number of defective 
areas increased by the 51st day. It can be noted that the destruction of contact zones around individual large sand particles 
is not observed. 

Images of the internal structure of samples of composition No. 3 with two-fraction sand (fractions 2.5—1.25 mm and 
0.63—-0.315 mm, no intermediate fraction) at the age of 0, 8, and 51 days are shown in Figure 5. 


2mm 


2mm 


Fig. 5. Internal structure of sample No. 3 (two-fraction sand): a, d — immediately after preparation; 
b, e —at the age of 8 days; c, f— at the age of 51 days 


In the images of the internal structure of the samples containing two fractions of sand, at the age of 8 and 
51 days, a clear formation of cracks around large grains of sand was observed, and the main mode of the crack formation 
was traced when the cracks were connected to each other (Fig. 6). 
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Fig. 6. Enlarged fragments of Figure 5: 
a — fragment 1; b — fragment 2; c — fragment 3 


233 


https://vestnik-donstu.ru 


234 


Puzatova AV, et al. Study of Structural Defects Evolution in Fine-Grained Concrete Using Computed Tomography Methods 


The formation of main cracks near coarse sand grains indicated that the contact zones around coarse aggregate 
particles were most stressed and susceptible to destruction during the shrinkage of the cement stone. 

Based on the obtained images of the internal structure of samples with different sand sizes, it can be concluded that 
the development of contact zone defects due to shrinkage of the cement stone depends on the ratio of fraction sizes. In 
samples with mono- and polyfraction sand, the defective structure develops locally, the area of such defects is visually 
much smaller than in samples with two fractions of different sizes. In the presence of a coarse fraction of sand with a 
high-volume concentration, defects in contact zones develop near coarse grains and have the main mode. 

Cube samples of composition No. 3 with two-fraction sand, modified by mechanical activation of individual 
components, were subjected to the application of external compressive load at the age of 3 and 7.5 hours after preparation. 
Images of the deformed internal structure of the samples are shown in Figure 7. 


Fig. 7. Images of the internal structure of modified samples: a — frontal section at the age of 3 hours; 
b — frontal section at the age of 7.5 hours; c — horizontal section at the age of 3 hours; 
d — horizontal section at the age of 7.5 hours 


The images showed the evolution of hardening concrete samples from elastic-plastic to brittle fracture. The samples 
tested 3 hours after production had an elastic-plastic fracture pattern, with clear cracks and chips on the sample edges. At 
the age of 7.5 hours, the sample edges were covered with a network of small cracks during fracture, and there were also 
numerous cracks and microcracks inside the sample, indicating brittle fracture. 

Discussion and Conclusion. Thus, by means of mechanical tests and X-ray computed tomography, it is possible to 
track the processes of destruction in the structure of fine-grained concrete. The results obtained give rise to a new complex 
method for assessing the structural characteristics of modified fine-grained concrete at all stages of strength development. 
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It has been established that the destruction of contact zones depends on the ratio of fraction sizes. In the bulk of coarse 
sand particles in the body of concrete, the destruction of contact zones is more pronounced and has the main mode. When 
using fine or polyfraction sand, contact zones are destroyed locally and have a visually smaller area. This indicates that 
contact zones near coarse aggregate particles are most stressed and are primarily subject to destruction under the shrinkage 
of cement stone. 

The studied dependences of the influence of the fine aggregate sizes on the mechanisms of formation and propagation 
of structural defects contribute to the theory of destruction processes of fine-grained concrete. The results obtained prove 
the prospects of using X-ray computed tomography as a method of non-destructive testing of the internal structure of 
fine-grained concrete, including at the early stages of strength development. Computer tomography, along with traditional 
methods of investigation of the structure and properties of building materials, gives rise to a new complex method that 
provides studying modern multicomponent concretes at all stages of strength development, the mechanisms of formation 
and propagation of structural defects due to natural processes of changing the state of the material and under various 
external loading conditions. 
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Abstract 


Introduction Violation of mutual positioning and fixation of parts worsens the operation of the equipment. Traditional 
approaches to solving the problem under consideration have been sufficiently studied: interchangeability of parts and the 
use of special equipment. Both methods involve a significant number of additional elements and assembly operations. 
Fixation is often provided by means of force fitting and welding. Disadvantages of these methods include assembly, 
residual and other stresses, engineering constraints, etc. To solve these problems, alloys with thermoelastic phase 
transformations are used, which provide shape memory effects (SME) to manifest themselves. This article describes, for 
the first time, self-positioning and self-fixation using the example of parts specially made from an alloy with SME. 

Materials and Methods. The pin element under pressing mandrels the blind hole of the cup and enters the seat. The alloy 
with SME was Ti-55.7wt%Ni. The temperature of the onset of its austenitic transformation was A; = 95°C + 5°C. The 
elemental composition was determined by a Shimadzu EDX-8000 X-ray fluorescence spectrometer, the phase 


composition — by a Shimadzu XRD-7000 diffractometer. The temperature was specified through differential scanning 
calorimetry. The range was 20—300°C, the heating rate was 5 deg/min. A Guide T120 thermal imager and a RangeVision 
DIY 3D scanner with structured illumination were used. After pressing the pin into the cup at different angles, the 
alignment and deviations between the axes of the cup and the pin were examined. Then, the cup was heated to 
110—-120°C, cooled, and control measurements were taken. 

Results. Values of the deflection angle after pressing were 0.2—11°. With a rigid structure and an installation angle of 0°, 
the pin deflected in the mounting hole by 0.2—0.5°. The axes shifted and did not intersect. The pin was not always 
completely pressed in. This indicated uneven deformation of the metal and different stress values around the hole. Such 
a unit would soon fail. The pin took the required position after heating the cup to 110—120°C (this temperature was higher 
than at the end of the reverse martensitic transformation). The angular deviation of the axes was noted to be 0.03-0.1°. 
The maximum misalignment (0.04 mm) corresponded to high positioning accuracy. Heating during the reverse 
martensitic transformation created internal stresses that returned the initial geometry of the cup. They also formed the 
forces that positioned and fixed the pin in the hole. That is, it is the parts that provide positioning and fixation (this is self- 
positioning and self-fixation). 

Discussion and Conclusion. For self-positioning and self-fixation of parts due to the shape memory effect, it is necessary 
to avoid sharp transition lines between the surfaces of parts during design, select rounded corners or fillets, and get a clean 
surface without burrs. Self-fixation and self-positioning reduce defects and inaccuracies during assembly. The use of 
certain alloys increases the profitability of equipment production. 


Keywords: shape memory effect, thermoelastic phase transformation, self-positioning of a part, self-fixation of a part, 
mutual positioning of parts, shape restoration due to return stresses 
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AHHOTalna 

Beedenue. Hapyuienue B3aMMHOro MO3HIMOHMpOBaHHA HU PUKcaMN AeTael yxyquiaeT padoty oOopyzoBanua. Jlocta- 
TOUYHO H3YYeHbI TpayqHUMOHHBbIe MOAXOMbI K PeLWIeHHIO paccMaTpHBaeMOH TIpOOJIeMBI: B3AaHMO3aMeHACMOCTH JleTael 
MCIOJIb30BaHHe CHelHaIbHOH OCHAaCTKH. O6a MeTOa MIpeAMoMaraoT 3HaYMTeCIBHbIM OOBEM OMOIHUTEIIBHBIX IJIEMCH- 
TOB HM MOHT@KHBIX Ollepalui. Dukcal{MIo YacTO OOecMedHBAIOT C TOMOMI[bIO TOcaKH C HaTATOM UH cBapKH. HegoctaTKu 
39TUX METOJOB: MOHTAKHBbIC, OCTATOUHBIe U pyre HallpsKeHHA, TEXHHYeECKHe OrpaHvueHna U up. Jia pemienua yKa- 
3aHHBIX MpOOMeEM HCIONb3YIOT CIIaBbl C TEPMOYUPyrHMH (a30BbIMM TIpeBpalleHHAMH, KOTOPbIe MO3BOMAIOT TPOAB- 
JATBCA 9*(PeKTam Mamatu copper (QII®). B qanHo cTaTbe BiepBbIe OMMCaHbI CaMOTIO3HIMOHMpOBaHne HM CaModuKca- 
lua Ha IpuMepe JeTaei, ClewMasIbHO HV3FOTOBJICHHBIX 3 criaBa c OII@. 

Mamepuanoi u memoovi. Uccnenopanucb crakaHbl 43 criiaBa c DI1@ — Ti-55,7wt%Ni mp 3ampeccoske: LITbIpeBoli 
3IEMEHT JOPHYeT lyxoe OTBepcTHe H NoMaaeT B WOcayqouHoe MecTo. Temnepatypa Hayasia ayCTeHHTHOrO TipeBpalile- 
Hua — A,=95°C+5 °C. OnemMeHTHEIM cocTaB onpeyelaIH peHTTeHodsyopecieHTHbIM clieKTpoMetTpom Shimadzu 
EDX-8000, da3o0Bpii — qudpaxtometpom Shimadzu XRD-7000. Temnepatypy ompexesanu quddepenuHarbHon cKa- 
Hupyloulei KalopuMetpuen. J[uana30n 20-300 °C, ckopocts Harpepa — 5 °C/mun. 3ayeticTBOBasM TerIOBH30p Guide 
T120 u 3D-ckaHep co cTpykTypHpoBaHHbIM MoycBetoM RangeVision DIY. Ilocue 3anpeccoBKu ox pasHbIMH yriaMu 
IUTbIPA B CTaKaH MCCIeAOBAIM COOCHOCTh UM OTKJIOHCHHA MEY OCAMM CTaKaHa H WTbIpsA. 3aTeM CTaKaH HarpeBasIN oO 
110-120 °C, oxnaxyann M elas KOHTpOJIbHEIe 3aMepEl. 

Pe3yjemamol uccnedoeanua. 3Ha4eHUA yria OTKIIOHEHUA Tocsie 3armpeccoBKu — 0,2—11°. pu »xecTKOM KOHCTpyKuMu 
Wf yryie YCTaHOBKH 0° WITIPh OTKJIOHACTCA B MOCaQOUHOM OTBepcTHH Ha 0,2—0,5°. Oc cmMemjaroTca MW He HepeceKarorca. 
LU rpipb He BCerya MOJHOCTHIO 3alIpeCCOBbIBacTCA. ITO FOBOPHT O HepaBHOMepHOH AedopMallHu MeTaliIa H O pa3HbIX 10 
3Ha4eHHIO HallpsKeHHAX BOKpyT OTBepcTHA. Tako y3en ObicTpo BbliizeT u3 cTpox. LTEIpb 3aHuMaeT TpebyemMoe MoONI0- 
2%KeHHe Tocie HarpeBa cTakaHa oO 110-120 °C (aTa Temnepatypa BbILIe, 4eM B KOHI[e OOpaTHOoro MapTeHCHTHOrO Mpe- 
BpalweHua). OTMeTHIIM yroBOe OTKIOHEHHe Oceii — 0,03—0,1°. MakcuMasbHasa HecoocTHoctTs (0,04 mm) cooTBeTcTByeT 
BBICOKOM TOUHOCTH NO3HUMOHHpoBanHa. Harpesp mp oOpaTHOM MapTeHCHTHOM IIpeBpallleHHu co3jqaeT BHYTpeHHHe 
HallpsoKCHHA, BOSBpalalolHe MepBOHAauaJIbBHYy!O TEOMETPHIO CTaKaHa. OHM 2Ke POPMUPYIOCT YCUIMA, KOTOPble paciosa- 
TaroT H (PHKCHPYIOT WTbIPb B OTBEPCTHH. To ecTb MMCHHO AeTaIM OOeciHeuHBaIOT NOZSHUMOHUpOBaHHe UM (PUKcaLHIO (39TO 
CaMOMO3HUMOHHpoBaHHe UH CaModuKcalHa). 

O6cystcoenue u 3akuo4uenue. J[nA cCaMOMO3MIMOHMpOBaHHA HM CaMOUuKcalHH WeTael 3a cueT 9*pekta MaMsATH (POpMBI 
IIpH KOHCTpyHpoBaHHH CileyeT OTKa3aTbCA OT Pe3KHX JIMHUM MepexoOB MexKy MOBCPXHOCTAMH JleTaeH, BbIOUpaTb 
CKPYTJICHHBIe yIsibI WIM rare, JOOMBaTbCA YMCTOM MOBepXHOcTH 6e3 3aycenHues. CaMo*uKcallMA HW CaMOTMO3HIHOHH- 
poBaHHe coKpallaioT Opak HM HeTOUHOCTH mpu cOopKe. Ucnomb30BaHHe criaBoB c 3:pdpeKTOM MaMATH MO2KET MOBbILMATS 
PeHTaOeIbHOCTh MpOu3BOACTBA. 


Ksnouesbie copa: addexT namMaTu bopMbl, TepMoyripyroe cba30Boe MpeBpallenve, CaMOMo3HUMOHMpoBaHue jeTasu, 
caMO@UuKcallHA JeTaIH, B3AMMHOe MO3HI{MOHMpOBaHHe JeTasiei, BOCCTAaHOBIeHHe (PPOpMBI 3a CUT BOSBPaTHBIX HallpsKeHHH 


BuaarogapHoctu. Astop OnarogxapuT KaHaqMAaTa TexHHMyecKux HaykK, ouenta H.A. Moctaka u Hay4Horo coTpyqHHKa 
JladopaTopHu epcieKTHBHOrO MpoeKTHpoBaHHaA HedTera30Boro obopyyoBaHua KyOaHcKoro rocyyjapcTBeHHoro 
TexHOsOrM4eckoro yHuBepcuteta M.A. CamMapuua 3a NOMOLM|b B IWIaHHpOBaHHH VW TpoBeeHH 9KCIIepHMeHToB. ABTOp 
TakKe NPW3HATeeH peLaKUMOHHOM KOJWIErHU 2KypHasla U pelleH3eHTy 3a KOMIeTeEHTHYIO SKCMepTH3y HM peKoMeHalHu 
110 yJIYUWIeHHIO CTaTbH. 
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Introduction. Self-positioning should be understood as spatial orientation, installation and positioning of a part 
and elements mating with it in units and mechanisms. Self-fixation is the interface along the surfaces of structural 
elements with such geometric features that provide fixation of the mating parts and a given fixation force due to the 
shape memory effect. 

The accuracy of the mutual positioning of parts determines the correctness and reliability of the mechanisms. Valid 
positioning can reduce the inertia of the mechanical unit, it prevents the occurrence of system backlashes not provided by 
the design. Failure to comply with the requirements for the mutual positioning of parts can change the operation of the 
products, which is unacceptable in most cases [1]. 

The requirements for the accuracy of assembly of units are high in all industries, specifically, in machine-tool 
manufacture, aircraft, and shipbuilding. For correct, tight mating of parts, the method of group interchangeability is used. 
Custom tooling is often used when assembling units. As a rule, it is unique for this type of assembly operation. It provides 
positioning accuracy during installation and meets tolerance requirements [2]. 

In industrial practice (notably, in construction [3] and the oil and gas industry [4]), alloys with thermoelastic phase 
transformations, which provide for the shape memory effect (SME), are increasingly used. The uniqueness of such alloys 
is in the features of phase transformations, due to which shape memory effects are manifested, as well as superelasticity 
(pseudoelasticity) [5]. 

Shape memory in alloys with thermoelastic phase transformations is used, in particular, for power drives [6]. The effect 
is based on return stresses that restore the shape of the part. The phenomenon corresponds to the austenitic phase state of the 
elements of the power drive [7]. In the oil and gas industry, coupling joints of parts are used for pipe junction [8]. The method 
involves restoring the shape of elements that provide the envelopment and compression of mating parts [9]. Due to the shape 
memory effect during assembly, it is possible to provide such processes as spatial orientation, installation, and positioning. 
To do this, it is advisable to use two principles mentioned above at once: 

— embracing and squeezing the parts; 

— restoring the shape due to return stresses causing movement. 

This approach allows reducing the number of assembly stages and avoiding the use of dedicated equipment (tooling) 
when positioning and fixing a part made of an alloy with a thermoelastic phase transformation and for the elements of the 
unit mating with it. 

The proposed solution opens up the possibility of replacing such fixation methods as force fitting and welding. 
Consequently, it is possible to avoid the disadvantages of these methods: assembly, residual and other stresses. We also 
note the operating restrictions associated with the strength and reliability of the fixation. In addition, it is not always 
technically possible to perform force fitting and welding. 

It is necessary to point out another advantage of parts made of alloys with thermoelastic phase transformations. Their 
positioning and fixation in the assembly unit allows refusing group interchangeability of parts, from fitting parts and from 
additional use of control units and mechanisms in the design. This simplifies the installation and manufacture of parts 
with high positioning accuracy and reliable fixation of parts in the assembly unit. 

The presented research article, for the first time, proposes to use thermoelastic phase transformation for positioning 
and subsequent fixation of parts in an assembly unit. In this case, the shape memory effect of parts ensures their self- 
positioning and self-fixation. 

The author of this article observed the phenomenon under thermoelastic transformations, when it was required to 
provide fixation due to the shape memory effect of teeth. In this context, the following were considered: 

— drill bit roller cutter; 

— ball plug in the roller bit leg; 

— roller cutter on the bit leg in the design of a roller bit without a ball plug; 

— nipple and couplings of a tool joint [10]; 

— seats enclosed in the disc rotary valve; 

— disk on the stem of a disc rotary valve. 

The research is aimed at studying the possibilities of self-positioning (spatial orientation, installation, and positioning) 
and self-fixation during the manifestation of the shape memory effect of alloys with thermoelastic phase transformations. 
It is expected to find out what accuracy is provided in this case during the design and installation of assembly units, and 
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what conditions need to be maintained during assembly (design features, techniques, requirements for the implementation 
of self-positioning and self-fixation). 

Materials and Methods. For the experiments, two samples were made from an alloy with a thermoelastic phase 
transformation. In shape, these were cups with a special blind hole and a pin element made of heat-treated (hardened) 
40X steel. The cup had an opening — two conical steps with a right (lead-in) and inverted (fixing) cone. Transitions — 
by radius and by the line of intersection of cones. The pin element had a mandrel shape, so when pressed into the cup, a 
blind hole was mandrelized, and the pin element got into the seat in the cup. All surfaces were clean, without corners and 
burrs. This provided: 

— free sliding of the pressed-in pin element; 

— spatial orientation and installation in the seat under the pin element in the cup; 

— positioning of the pin element in the seat. 

Figure | shows the stages of pressing a pin element into a cup made of an alloy with SME. 


Fig. 1. Installation of pin element into seat of cup made of alloy with SME: 
a — 3D model of elements assembling after installing pin element 2 into cup 1; 


b — assembly drawing; c, d — installation of pin element 2 and cup | in hand press 3 before pressing the pin element into cup; 
e — angle between axes of symmetry o after pressing pin element 2 into cup | before heating (room temperature — martensitic 
phase of alloy); f— pin element and cup at room temperature after heating to 110°C 
(angle between axes of symmetry 0°, austenitic phase of the alloy) 
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Titanium nickelide Ti-55.7wt%Ni was selected as the alloy with SME. The temperature of the onset of its austenitic 
transformation was A; = 95°C + 5°C. Three properties listed below were determined. 


1. Elemental composition. For this purpose, X-ray fluorescence analysis was performed using a Shimadzu EDX-8000 
device (Japan). 

2. Phase composition. X-ray phase analysis was performed using a Shimadzu XRD-7000 X-ray diffractometer (Japan). 
Parameters: Cu Ka — 1.54 A°, 40 kV, 30 mA, angle range 20-90 degrees, shooting speed 1 deg/min. 

3. Temperature. It was determined by differential scanning calorimetry. The range of work was 20-300°C. Heating 
rate — 5 deg/min. 

This alloy allows for demonstration of the general principles of self-fixation and self-positioning of a pin element in 
the seat of the cup during the phase transition from martensite to austenite. At room temperature, the alloy is in the 
martensite phase. 

A Nordberg manual press was used for pressing. The cup made of an alloy with SME was heated with a GHG 23-66 Bosch 
technical hot air gun. Temperature control provided uneven heating of the cup and maximum imitation of the assembling process 
at enterprises. Specifically, the worst installation conditions and operation with violations of the process were recreated. 

The temperature and hot zones, as well as heating over the entire surface to the specified temperature, were monitored 
using a Guide T120 thermal imager. The positioning accuracy was determined by a stationary 3D scanner with structured 
illumination RangeVision DIY. Daheng cameras provided a measurement accuracy of 0.02 mm. 

The sequence of operations for installing the pin element into the cup is described below. 

1. The pin element is pressed into the blind hole of the cup. As a result, the hole is mandrelized. It deepens until the 
pin element enters the seating surface of the hole. The process takes place at room temperature, corresponding to the 
martensitic phase state of Ti-55.7wt%Ni alloy with shape memory effect. 

2. The relative position of the parts is measured. 

3. The cup is controllingly heated until the temperature is higher than it was at the end of the austenitic transformation. 
At this temperature, the initial (before deformation) shape of the cup is restored, and the force from the resulting return 
stresses is sufficient to shift and orient the pin element inside the blind hole of the cup. After this, the cup is cooled to 
room temperature. 

4. A control measurement of the relative position of the parts is performed. 

The pin element was pressed into the cup at different angles between the axes of symmetry of the cup and the pin 
element. The angle sizes were from 0 to 12 degrees. The step was 3 degrees. The alignment and the angle of deviation 
between the axes were determined. Then the cup was heated to 110—120°C and cooled in various ways — from free 
cooling in air to forced cooling (by lowering into water). Then, a control measurement of the alignment and the angle of 
deviation between the axes was made. 

Research Results. After pressing, the angle of deviation between the axes of the cup and the pin element had values 
in the range from 0.2° to 11° (angle o in Fig. 2 a). This indicates that with the rigidity of the structure and the installation 
angle of 0°, some displacement of the pin element in the mounting hole still occurs (from 0.2° to 0.5°). 
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¢) d) 


Fig. 2. Deviations arising during installation of pin element in cup mounting hole: 
a — angle between the axes of symmetry; b — displacement as a result of underpressing; 
c, d— displacement of the axes of symmetry in the section plane perpendicular 
to the axes of symmetry of the pin element and the cup 


The axes were offset and did not have an intersection point, as shown in Figure 2 c and 2 d. If we take the axis of 
the cup (reference point) as the center, the offset of the axis of the pin element will be p = 0.2—0.8 mm (Fig. 2 d), 1.e., 
the pin element in the cup will have an offset mounting position. In a number of cases, the pin element was not 
completely pressed into the cup (Fig. 2 6). Firstly, this caused uneven deformation of the metal around the hole. 
Secondly, it indicated that as a result of pressing, different stress values arose along the circumference of the hole. 
Such a unit will quickly fail when used. 

The pin element occupies the required position in the cup hole and relative to the hole seat after heating the cup to 
110—120°C, i.e., at a temperature higher than at the end of the reverse martensitic transformation. 

The positioning accuracy of the pin element in the cup was studied using a stationary 3D scanner with structured 
illumination RangeVision DIY and Daheng cameras. The angular deviation of the axes was recorded in the range 
of 0.03—0.1°. The maximum value of the misalignment was 0.04 mm, which corresponded to high positioning accuracy 
for assembly units. 

It should be noted that there is a similarity between the two mechanisms that provide: 

— occupation of the structurally required position of the pin element in the cup hole; 

— action of the working element made of an alloy with SME of the power engine. 

The specified shape of the part made of an alloy with shape memory effect of the cup in the austenitic state corresponds 
to the shape in which the pin element exactly occupies the position required by the design of the unit. At a temperature 
corresponding to the martensitic phase state, the mechanical characteristics of the alloy are lower than at a temperature 
corresponding to austenite. In this case, it is possible to specify a shape into which the pin element can be easily mounted. 
With such installation, the positioning of the pin element does not correspond to the required. 

Subsequent heating of the glass as a result of the reverse martensitic transformation creates internal stresses that allow 
the original shape to be restored. They can be conventionally called return stresses. In this case: 

— the cup takes on a shape that corresponds to the operational geometric characteristics; 

— the forces arising as a result of the action of return stresses are sufficient to move, orient, install and position the pin 
element in the seat of the cup (it is located and fixed in the hole). 

Evidently, it is the parts that provide self-positioning and self-fixation. No special equipment or special assembly 
techniques are needed. 

There is a process similar to the one described above. To operate the satellite antenna deployment mechanism, the 
power drive uses the force developed by the return stresses during the reverse martensitic transformation of the working 
element made of an alloy with SME. 
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Discussion and Conclusion. The displacement of the pin element inside the cup mounting hole is provided by 
design features that can be applied to all mating parts under self-positioning and self-fixation due to the shape memory 
effect of alloys with thermoelastic phase transformations. Techniques that allow maintaining the required design 
features are listed below: 

— rejection of sharp transition lines between surfaces of parts in favor of smoother ones; 

— selection of rounded or variable radius chamfers and external corners; 

— replacement of internal corners with fillets or curves with variable radius; 

— high-quality surface cleaning, elimination of burrs. 

These requirements should be taken into account when designing. This will provide the desired mutual 
displacements of parts with significantly less force. In this case, the force is compared to that given by the return 
stresses during the reverse martensitic transformation of the alloy with thermoelastic phase transformations (the cup is 
made of such material). 

Self-fixation and self-positioning can maintain assembly accuracy, eliminate mounting and engineering inaccuracies 
during installation work, and reduce defects. The use of some alloys has economic substance, i.e., it works for the 
profitability of the final product [11]. These and other advantages of the approach described in the article can be 
implemented in machine tool manufacturing. The method seems useful for developing oil and gas, precision and other 
equipment with high requirements for the geometric accuracy of parts. 
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Abstract 

Introduction. When drilling oil or gas wells, rock from the bottomhole is brought to the surface by drilling mud, which 
is cleaned of sludge by shale shakers at the first stage. Shale shakers are equipped with such a screen and create such 
trajectories of vibration of the frame in order to solve a dual problem: on the one hand, to provide the cleaning of the 
drilling mud coming from the well mouth, and on the other hand — to maintain the proper degree of cleaning. To correctly 
select the nomenclature of screens, it is necessary to reliably determine the throughput capacity of a shale shaker in the 
real-time setting. This will allow you to pre-order the required number of screens with the required cell size. Previously, 
studies were conducted by sieving a mixture of calibrated granules on a standard shaker or by straining mineral oil through 
a fixed screen. However, this does not fully correspond to the actual conditions of drilling mud screening. The objective 
of this article was to simulate the throughput of a drilling shale shaker under real conditions by calculating the movement 
of drilling mud along a vibrating screen using its previously obtained specific throughput capacity. 

Materials and Methods. When creating a mathematical model of the full throughput capacity of a shale shaker, a previ- 
ously known mathematical model of fluid flow in an open channel and a finite-dimensional approximation in the form of 
a sequence of sections of concentrated containers connected by hydraulic conductivities (by analogy with the approxima- 
tion of long electrical lines and extended gas pipelines by the finite element method) were used. The fluid flow rate over 
the screen was determined according to Chézy law. In this case, the cross-sectional area of the flow above the screen for 
a specific section was specified as the arithmetic mean between the initial and final values of the drilling mud height in 
the corresponding section. 

Results. A chain calculation scheme for sieving drilling mud was compiled. Based on the material balance, mathematical 
models were developed for determining the throughput and height of drilling mud above the screen: 

— for the first section of the screen and the first nodal point of the screen; 

— for subsequent screen sections and subsequent screen nodal points; 

— for the last section of the screen and the last nodal point of the screen. 

A mathematical model of the drilling mud movement along the shale shaker screen has been obtained, which, due to the 
finite-dimensional representation of this flow along n sections of the screen length, is a system of n integral equations 
describing concentrated tanks, and n — | algebraic equations describing the flow of drilling mud between tanks. 
Discussion and Conclusion. The movement of drilling mud along the length of the shale shaker screen corresponds to 
the movement of a steady uneven flow in an open channel. However, in the first case, the volume flow along the length 
is variable (decreasing to zero), and in the second case, it is constant along the length. Therefore, the use of only one 
differential equation for steady uneven flow in an open channel is not sufficient in this case. The resulting mathematical 


model of the drilling mud movement along the shale shaker screen is a mathematical basis for the development of a 
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modeling program to determine the throughput of the shale shaker for real operating conditions, 1.e., for a specific drilling 


mud, a given screen or a given drive system of the shale shaker. 


Keywords: shale shaker, vibrating screen, height of the mortar layer, material balance, throughput capacity, edge condi- 


tion, nodal points 
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AHHOTalna 

Beedenue. [pu Oypenun HepTaubix WIM Ta30BbIX CKBAKHH MOpoOa U3 300A BEIHOCHTCA Ha HOBEPXHOCTb OypoBbIM pac- 
TBOPOM, KOTOpbIM OUMMaeTCA OT WiIaMa Ha TepBol cTayquH BHOpocutTamu. BuOpocnta ocHaljaroTca TaKoli CeTKOM H 
CO3al0T TakHe TpaeKTOpHu BHOpoKoNeOaHHli paMBbl, UTOObI C OAHOM CTOPOHbI — OOecHeuHBaTb OYMCTKY MOCTyMarouyero 
W3 YCTbA CKBaKHHbI OypoBoro pacTBopa, a c Apyrol — oOecieunBaTb JOMMKHYIO CTeleHb OUNCTKH. J[It KOppeKTHOTO 
ToOopa HOMeHKAaTYPbI CeETOK HYKHO HaJeKHO OMpeeATh MPOMYCKHY!O CIOCOOHOCTb BHOPOCHTa B peaJIbHBIX YCIIO- 
BMAX. ITO MO3BOMMT 3apaHee 3aKa3aTb TpeOyeMOe KOJIMYECTBO CETOK C HY2KHbIM pa3MepoM sYeeK. PaHee HccIeqOBaHHA 
IIPOBOAWIMCb IyTeM IIpoceHBaHHA CMeCH KasIMOpOBaHHbIX rpaHyJI Ha CTaHAapTHOM BHOpocTeHAe WIM MpowexuBaHHeM 
MHHepasIbHOrO Macyla Yepe3 HeMOABYWKHY!IO CeTKY. OWHaKO 9TO HE B MOJHOM Mepe COOTBETCTBYECT PeasIbHbIM YCIOBHAM 
lipocenBanna OypoBoro pacTBopa. Lesbo WaHHO paooTEI ObLIO MOZeMpOBaHHe MponycKHOH ciocoobHocTH OypoBoro 
BHOPOCHTa B PeaIbHBbIX YCIOBHAX IyTeM pacueTa JBYWKeHHA OypoBoro pacTBOpa 10 BHOpUpyroleli CeTKe C HMCIIOJIb30- 
BaHHeM eé paHee MOJyYeHHOH yAebHOM MpOmycKHOH crocoOHocTH. 

Mamepuanot u memoooi. IIpu co3qaHun MaTeMaTH4eckKOl MOJeIM NOHO MpolycKHoH clocoOHocTu OypoRoro BHOpo- 
CHTa MCHOIL30BaIacb H3BeECTHad paHee MaTeMaTHUeCKad MOJeJIb TEXCHHA X9KHAKOCTH B OTKPbITOM KaHasle H KOHC4HO- 
MepHas arllipOKCHMallHA B BUe MOCICOBaTeIbHOCTH Y4aCTKOB M3 COCpeOTOUCHHBIX EMKOCTeH, COC AMHEHHBIX THApaB- 
JIMYCCKUMH TPOBOAMMOCTAMH (10 aHaIOIH C alIIIPOKCHMallHeli IICKTPHYeCKHX JVIMHHBIX JIMHUM UM WpOTs.KeCHHBIX Ta30- 
IIPOBOOB MCTOJOM KOHEYHBIX 3JIEMCHTOB). Pacxoy, TedeHHA 29KMAKOCTH Hay CeTKOM onpeyensaica 10 3aKony Ile3n. pu 
9TOM TIOWayb CeCYeHHA MOTOKA Hal CeTKOM JIA KOHKpeTHOTO y4yacTKa OMpeyeraacb KaK cpeqHee apudMeTuyecKoe 
MOXY HavasIbHbIM MW KOHCUHBIM 3HAYCHHEM BBICOTHI OypoBOro pacTBOpa Ha COOTBETCTBYIONeM Y4acTKe. 

Pe3yiemamot ucciedoeanua. CoctaBieHa WelHad pacueTHad CxeMa IIpoceuBaHua OypoBoro pacTBopa. Ha ocHoBaHun 
MaTepHaJIbHOro OaslaHca paspaOoTaHbl MaTeMaTHYeCKHe MOI OMpeeeHUA IpOIMyCKHOM CIOCOOHOCTH H BBICOTEI Oy- 
poBoro pactTBopa Haj ceTKOH: 

— [JIA HWepBOro y4acTKa CeTKH H MepBOH y3IOBOM TOUKH CeTKH; 

— [IA MOCIICLYIOUIMX YAACTKOB CeTKH HM MOCICAYIOWIMX Y3IIOBbIX TOUCK CETKH; 

— [IA MOCHIeMHerO y4acTKa CeTKH HU Kpalineli y3IOBOM TOUKH CeTKH. 


Tosyuena MaTeMaTuuecKad MOJCJIb TBYOKCHHA OypoBoro pacTBOpa II0 CeTKe BHOpocuta, KOTOpad BCIICACTBHe KOHCUHO- 
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MepHOPo IipeACTaBJICHHA STOO TWOTOKa 0 1 YAaCTKAaM JVIMHbI CCTKH NpeCTaBJIaAeT cobok CHCTeMY U3 1 UAT pasIbHbIx 
ypaBHennii, OTIMCBHIBAIOWIHX COCHpeOTOUCHHbIC CMKOCTH, Un — 1 anreOpawueckux ypaBHennii, OMMCHIBAaIOWINe WepeTOKH 
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O6cystcoenue u 3akur04enue. J\ByoKeHue OypoBOro pacTBOpa IO JWIMHe CeTKH BHOpOCHTa COOTBeTCTByeT JBHKCHHIO 
YCTaHOBUBINerocd HepaBHOMepHOro IOTOKa B OTKPbITOM pycsie. OWHakO B TepBoM ciry4ae — OObeMHBIM pacxoy, 110 
JJIMHe TlepeMeHHbIi (yORIBaIOINIHH 0 HyJIA), a BO BTOPOM cilyuae — NOCTOAHHBIM Mo WIHHe. Ilootomy Ucnomb30BaHHe 
TONbKO OAHOTO ZupdepeHuMasbHOTO ypaBHeHHA yCTaHOBMBIIeroca HepaBHOMepHOLO IMOTOKAa B OTKPBITOM pycile B jaH- 
HOM cylyyae HeJOCTaTOUHO. 

IlomyueHHad MaTeMaTH4UecKad MOJesIb JBMXKCHHA OypOBOTO pacTBOpa IO ceTKe BHOPOCHTa ABIIAeCTCA MaTeMaTH4YeCKHM 
oOOocHOBaHHeM JIA pa3pa0oTKH MOJeMpylollei MporpaMMBbI 10 olpeeseHurO IPOMycKHoH clocoOHocTu BHOpocuta 
JIA PeasIbHBIX YCIOBHM ero paOoTHI, TO eCTb Id KOHKpeTHOrO OypoBoro pacTBopa, JaHHOM CeTKH WIM 3ayjaHHol cH- 


CTEMBI IIPHBOJIOB BHOpocnta. 


Ksrouesbie c10Ba: OypoBoe BHOpocHTo, BHOpHpyloljat CeTKa, BbICOTa CJIOA pacTBOpa, MaTepHaJIbHbI OamaHc, 


TIponyckHaa ciHocoOHocTs, TpaHu4Hsle YCJIOBHA, Y3JIOBbIe TOUKH 
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Introduction. The capacity of the drilling mud cleaning unit of the rig circulation system is limited by the flow rate 
of the circulating drilling mud, which must be cleaned properly [1]. The capacity is the volume of drilling mud sieved 
through the screen, and the flow rate is the volume of drilling mud that flows without sieving. The shale shaker is the first 
processing vessel in this unit, and it performs primary purification. Therefore, its throughput mainly determines the 
capacity of the entire cleaning system. 

The capacity of the cleaning unit of the circulation system (CS) of a standard drilling rig should exceed the throughput 
of the drilling pump and be about 30-80 1/s. Insufficient capacity of the CS limits the rate of penetration [2]. 

The circulating drilling mud is fed to the shale shaker directly from the well mouth via a chute. There are no 
intermediate tanks; therefore, even if its consumption briefly exceeds the shale shaker's capacity, the drilling mud will 
flow through the screen into the tank waste, which will result in its loss. If the driller's assistant opens the gate valve, the 
raw mud will flow past the vibrating screen into the tank. This is an extremely undesirable case, since the sludge settles 
in the tank, preventing it from being pumped into the desander. This dramatically reduces the degree of drilling mud 
purification, which leads to increased drilling costs and even complications in the well [2]. 

The throughput of a particular shale shaker is largely determined by the size of the screen cell with which it is 
equipped. With a decrease in the cell size, the throughput of the shale shaker decreases significantly, but the degree of 
cleaning of the drilling mud increases. In practice, to prevent the solution from passing by the shale shaker, a screen with 
a larger cell with a large but unknown margin of throughput is installed on it [1]. 

The actual measurement of the throughput of the shale shaker is difficult, because after sieving through the screen, 
the mud flows into the tank under the shale shaker by gravity, without passing through the pipeline. Therefore, the 
throughput of the shale shaker was determined by indirect methods. 

In [2], based on production experience, the parameters that change the throughput are indicated: an increase in the 
length of the screen, a decrease in the flow rate, a decrease in the angle of inclination of the screen, a change in the 
direction of particle movement, a decrease in the amplitude of screen vibrations, the simultaneous use of two consecutive 
or parallel screens. However, the authors do not provide quantitative indicators of these parameters. In [3], values of the 
shale shaker throughput when working with various drilling fluids are given, but the general formula is not given. In!, 


only the formula for calculating the change in capacity is indicated. The problem of selecting a screen for a specific case 


' American Association of Drilling Engineers Shale Shakers and Drilling Fluid System. Houston: Gulf Publishing Company; 1999. 335 p. 
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is solved in a review by a team of American engineers’. It consists of selecting a screen with the smallest cell, at which 
the shale shaker would have a throughput greater than the specified one. Several methods have been developed for testing 
screens, but all of them are far from the operating conditions of shale shakers on drilling rigs. Thus, in the API 
Recommended practice 13C standard, screens are tested either by sieving a mixture of calibrated granules on a standard 
shaking table or by filtering mineral oil through a stationary screen. In [4], an attempt was made to analytically describe 
the movement of the drilling mud flow through the screen and the process of sieving it through the screen based on the 
Navier-Stokes equation. The simulation results do not take into account the dynamic parameters of the screening machine. 
These parameters are taken into account in [5], but as a parameter in question not for the drilling fluid, but for coal sludge. 
The amount of all the sludge removed from the total mass is relatively small [6]. Such a concentration of sludge does not 
have a noticeable effect on the rheology of the drilling mud [3]. Accordingly, the separation of the dewatering process 
given in [7] is inappropriate. A high concentration of sludge is observed in the drilling mud leaving the desanders and 
desilters, which is cleaned on the second shale shaker?. 

From all has been said it follows that it has not yet been possible to determine the value of the drilling shale shaker 
capacity using theoretical studies. Therefore, the objective of this work was to model the full throughput capacity of the 
shale shaker based on previously conducted laboratory experiments to determine the specific capacity. 

Specific throughput capacity 

Previously, the team with the participation of the author developed a methodology [7], technical and software support 
for the experimental determination of a mathematical model of the specific throughput of a specified screen of the first 
shale shaker for a real drilling mud for a special shale shaker, and experimental studies were performed [8]. If all the 
points of the support frame located along the screen line have the same vibration trajectories, or, as they say in vibration 
engineering, “the orbits are homogeneous”, then it is enough to conduct one experiment. If the trajectories are 
nonhomogeneous, then it is enough to conduct three experiments with trajectories at the edges and in the center of the 
screen. In any case, the vibration trajectory parameters of the frame points of the production shale shaker must be known 
and implemented in the experiments. 

Experimental benches of other scientists are known, e.g., in [9], a bench that creates polyharmonic trajectories of 
vibrations is described. This is achieved by the corresponding arrangement of the multi-speed drive system. During the 
experiments, even Lissajous figures were obtained [10]. The resulting best design may work less efficiently due to external 
actions. At the same time, the purpose of automatic control is to compensate for the negative impact of external actions. 
In [11], a control method is given, in [12] — the dynamic parameters of the shale shaker under control, and in [13], an 
example of step-by-step optimization using the Adams and Nastran modeling units is given. 

The experiments performed on sieving a bentonite-water solution with a density of 1.19 g/cm? through a screen sample 
with square cells measuring 0.40.4 mm allowed us to obtain the following mathematical model of the specific capacity 


of the screen [9]: 


q(h, Ay) = 0.0088-h'® - Ay'*?, (1) 


where g(h, Ay) — specific throughput capacity of the screen, m?/m’s;  — height of the rubble mortar layer on the screen, m; 
Ay — amplitude of vertical vibration acceleration, m/s’. 

Boundary conditions of the sieving process 

A mathematical model of specific throughput is not enough to determine the throughput capacity of a drilling shale 
shaker. It is also necessary to know the distribution of the height of the drilling mud layer along the length of the screen 
h(x) at a known flow rate of the solution entering the screen Q(0). It is not possible to take a linear distribution of height 
along the length of the screen as a first approximation, since the boundary condition at the end of the sieving process (0) 
is unknown. The boundary condition for the beginning of the sieving process, as described earlier, is known and has the 
form (0.8 - L) = 0, where L — length of the screen, m. 


? Tbid. P. 3-335. 
3 Tbid. P. 3-335. 
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In classical hydraulics, a mathematical model is known for the distribution of the height of the liquid layer along the 
length of the flow with a non-pressure steady uneven flow of liquid in an open channel [14]. However, the flow rate of 
the liquid in this case is constant along the length of the flow, and in our case, the flow rate decreases due to the sieving 
of part of the drilling mud. Therefore, the steady motion of the drilling mud flow along the shale shaker screen is an object 
with variables h(x) and Q(x) distributed along the length of the screen. 

Analysis techniques for distributed objects 

The search for engineering solutions to such problems has shown that in electrical engineering, long electric lines are 
replaced by a chain of concentrated resistances, capacitances and inductances [15]. Extended gas pipelines are replaced 
by a chain of alternating tanks and hydraulic resistances to simulate their operation [16]. This approach is especially often 
used in machine dynamics, where chain calculation schemes are built according to general drawings, consisting of 
alternating concentrated masses or moments of inertia and concentrated pliability for vibrating machines [17] and control 
systems [18]. This is essentially a finite element approximation of distributed objects. 

Construction of a mathematical model of the full throughput capacity of a shale shaker 

Let us construct a calculation scheme of the drilling mud flow from concentrated capacities and concentrated hydraulic 
conductivities, shown in Figure 1. To do this, we divide the length of the shale shaker screen by nodal points into equal 
sections of length Ax = L/m, where m — number of sections. Nodal points must be at the beginning and end of the screen. 
We place concentrated capacities in the nodal points. The initial and final capacities correspond to a section of the flow 
of length 0.5 - Ax, and the remaining capacities correspond to Av. 

L 


— = 


Fig. 1. Chain calculation scheme of the hydraulic system of drilling mud flow along the screen 


The material balance equation for the first container in Figure | has the form: 
Ax 


(QO: (t)-Q2 (t)- On (t))-dt = : -b-dhi(t), (2) 


where Q(t) — flow rate of the drilling mud entering the shale shaker, m3/s; Q.(t) — flow rate of the drilling mud 
flowing out of the screen from the tank of the 1st nodal point to the tank of the 2nd nodal point, m?/s; Q,1(¢) — flow 
rate of the drilling mud sieved in the first section, m?/s; 41(¢) — height of the mud layer at the first nodal point, m; 
b — screen width, m. 


The sieved flow rate at the first nodal point, taking into account the specific throughput, is equal to: 
Ax 
Om (t)=q(n (¢), yi) 5-5, (3) 


where Ay, — vibration acceleration amplitude at the first nodal point. 
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Substituting expression (3) into (2) and transforming the equation, we obtain the expression for /(¢): 


n()= || “0, (t)-02(0) a(hi(t).Av) (4) 


Ax-b 


The material balance equation for the capacities of all intermediate nodal points in Figure | has the form: 
(Oi (t)— Ois1 (t)- Oni (t))- dt = Ax -b- dh; (t), (5) 


where Q,(t) — flow rate of drilling mud flowing into the tank of the i-th nodal point, m/s; Q;+:(¢) — flow rate of drilling 
mud flowing along the screen from the tank of the i-th nodal point to the tank of i+ 1st nodal point, m°/s; 
Ont) = q(hit), Ay) - Ax - b — flow rate of drilling mud flowing along the screen from the tank of the i-th nodal point to 
the tank of i+1st nodal point, m*/s; A(t) — height of the solution layer at the i-th nodal point, m. 


By transforming equation (5), we obtain the expression for all intermediate nodal points: 


(0) ff Lp -2.a(0))-a( h(t.) fa 6) 


Ax-b 


At the last nodal point, based on physical considerations, we take the height of the solution layer to be zero (C), for 
not to support the flow in the event of its flow into the tank waste. This is essentially the right boundary condition for the 
height of the solution layer. The flow rate of the drilling mud entering the shale shaker is assumed to be constant in this 
mathematical model and equal to the specified value Q)(4) = Q1 say. This is an indirect assignment of the left boundary 
condition for the height of the solution layer. 

All intermediate costs in Figure 1 are determined by the difference in heights of the solution layer h;,(4 — h(#) and the 
concentrated hydraulic resistance of the screen section. Drilling mud is a non-Newtonian fluid, but with the turbulent flow 
of clay solutions, it is possible to proceed from the usual formulas of hydraulics [19]. 

To determine the flow rate QO), we do not specify the value of the hydraulic resistance of a section of the screen with 
a length of Ax, but we use the well-known differential equation for a steady-state non-uniform flow of constant flow 


rate [20], but for a raised screen: 


QO? 
ee 
dh(x)__"" F2.C7-R (7) 
d. atl 
7 ee 
gF 


where i, — geometric slope of the flow channel, in Figure | i, >0i, >0; F — cross-sectional area of the flow, m’; 
C— Chézy coefficient; R — hydraulic radius of the flow cross-section, m; « — Coriolis coefficient, « = 1.10 — 1.15 [21]; 
g — gravitational acceleration. 

Let us express the derivative in differential equation (7) by finite differences and use the variables from the calculation 


scheme in Figure | in it: 


2) 
Ina(t)-mi(t)_ "FP? -C?-R, s 
Ax oO Q;’ (t) 
J-—. =. 
ag sc 


where F; — cross-sectional area of the flow in the section from the 7 — 1st nodal point to the i-th nodal point. 
Here, it should be clarified that Ot) — flow rate in the section from the i — 1st nodal point to the i-th nodal point. The 

cross-sectional area of the flow in this section varies from /;-:(f) b to h(t) b. Therefore, we define value F; as the average 

value of these areas: 

hy (t)+ hi (t) 


FF; =b- 
2 


(9) 
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Value of the hydraulic radius of the flow is equal to the ratio [15]: 


Rot! _ py fal thi(t) 1 (10) 


Aj 2 pe 


where A; — length of the wet perimeter of the flow in the section from the i — 1st nodal point to the i-th nodal point, m. 
The width of the industrial shale shaker screen is from 1 meter, the thickness of the solution layer on the screen is 
several centimeters [4]. Therefore, from expression (10), it follows that R; ~ (Aii() + A{d)/2. Substituting this value and 


value F; from expression (9), we obtain: 


ip + Qi" (*) 
a(t) (05-(hea(e)+h()) -b*-C* “ 
= (eo Qi" (t) 
& (0.5-(hia(t)+hi(t))) <b? 
We solve equation (11) with respect to Of) for i> 1: 
hy-1 (t)-hi (t) ee 
{i= Ax TGS he ' * pe? 
Q;(t) moe 4 [.0.5-(di-a(t)+hi(t)) | +B. (12) 
Ax ae Oe 


Discussion and Conclusion. The mathematical model of the flow of drilling mud through the shale shaker screen is 
a system of n equations (6) and n — lequations (12). 

The initial conditions for variables {0) cannot be taken as zero, since this will lead to division by zero in equations (12). 
Therefore, it is recommended to take these conditions as non-zero, but sufficiently small, e.g., 4(0) = 0.0001. 

Modeling of the system of equations (6) and (12) can be performed in simulating environments, such as MATLAB 
with the Simulink or SimInTech application. In any case, it is required to determine the number of nodal points of the 
chain calculation scheme in Figure 1. The throughput capacity of a shale shaker equipped with a particular screen is equal 
to the flow rate Q) 3a,, at which the last 20-25% of the screen length would be free of drilling mud and used to dewater 
the sludge’. Therefore, the penultimate nodal point with a minimum number of these points should be at the same distance 
from the edge of the screen. This implies Ax ~ 0.15 - L, what the seven nodal points correspond to. 

In this work, the task was to determine the steady-state flow of the drilling mud, or more precisely, its profile. The author 
was forced to use the transient process of this flow, since it was impossible to establish the left boundary condition /)(:°), 
instead of which an indirect assignment of this boundary condition Q1(f) = Qi sa, was used. The final number of nodal points 
is determined only in the process of modeling based on the accuracy of the drilling mud profile. For this purpose, the 
simulation of mathematical models (6) and (12) with n and n + | nodal points is performed. If the greatest difference in the 
flow profiles exceeds, e.g., 5% of the value of the mud layer height at this point, then the number of nodal points is increased 
by one to m + 2 [22]. Then, a simulation is performed with n + 2-nd node points and the flow profiles of this simulation are 
compared to the result of the previous simulation with n + 1 points. Upon reaching the specified error, the mathematical 
model with the latest value of the number of node points is used in the future. 

The obtained mathematical model of the drilling mud flow along the shaker screen, together with the previously 
obtained mathematical model of the specific throughput capacity of the screen for specific drilling mud [8], makes it 
possible to reasonably recommend the installation of a screen with certain cell sizes for various costs of drilling mud 


flowing from the well mouth. 


4 American Association of Drilling Engineers Shale Shakers and Drilling Fluid System. Houston: Gulf Publishing Company; 1999. 335 p. 
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Introduction. The scientific literature examines various approaches to quality management in information 
technology (IT). The issues of identifying and correcting defects are worked out, and the possibilities for minimizing 
them are shown. There are materials on quality management in complex engineering processes. At the same time, there 
is no detailed description of quality management at each stage of the IT product life cycle, including testing. It should be 
noted that the coordination of software releases is closely related to quality management, but this process is rarely or 
fragmentarily considered in the literature. Additionally, the interprocess communication is not taken into account; 
therefore, there is no comprehensive understanding of quality management in the creation, testing and refinement of 
software. This study is designed to fill these gaps. The research is aimed at presenting a comprehensive approach that 
links the theory, practice and methods of software quality management. 

Materials and Methods. Theoretical and applied literature on the subject were studied, analyzed, and reviewed. The 
author's professional experience in managing the quality of IT products was used. The practices of global suppliers of 
digital goods and services were taken into account. The author has used these materials and methods to study in detail the 
issues of software testing and code deployment. 

Results. A comprehensive model of quality management in software development is elaborated, described and presented in 
the form of a diagram. Its interconnections with the project management model and the product life cycle, namely: analysis, 
design, development, testing, deployment, and support, are identified. Principles of quality management at each of these 
stages are specified. The processes and checks during code deployment are systematized and presented in the form of a 
diagram. Their features are shown in three environments: during development, testing, and production. The algorithm allows 
quality experts to build the sequence of actions to eliminate detected defects in the future, understand the situation when it is 
possible (or impossible) to deploy code and determine the moment when the software should be transferred to the user. In 
addition, the proposed scheme can be the basis for automating code deployment. The solution will reduce development time. 
As a result, the product will enter the market faster, which will speed up the payback of costs. 

Discussion and Conclusion. The implementation of the model created within the framework of this scientific work into 
the production practice of IT companies presupposes strategic changes. Their implementation requires significant time 
and other resources; therefore, the overall transformation process should be divided into parts. The proposed approach is 
adaptable to the needs of different organizations and products. You can work with individual components to create an 
optimal plan for achieving quality management goals. 


Keywords: code deployment, quality management in IT, IT product life cycle, IT product testing, software release 
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AHHOTalna 

Beedenue. B naywHol WuTepatype paccMaTpuBaloTcA pasHble MOAXObI K MCHePKMeHTY KayecTBa B cibepe HHPOpMal{HOHHBIX 
TexHosiorni (MT). popadotaHnl Bompocnl BbIABJICHHA H HCiIpaBIeHHaA Je:eKTOB, 1OKa3aHbI BO3MO2KHOCTH HX MHHMMV3allHH. 
Ectb MaTepnasibI 00 yipaBJIeHHH KadeCTBOM B CJIOXKHBIX TEXHOJIOTM4eCKHX IIporeccax. JloKa3aHo, YTO paboTa c KayeCTBOM 
IM@poBblx MpOyKTOB TpebyeT B 4YHcIe Tpouero NpoACHeHMA BOTIpPOCoB KayecTBa Koga. IIpu sTOM HeT AeTasIbHOrO OMCAaHHA 
ylipaBJIeHHa KaYeCTBOM Ha K@KJOM 9Talle 7KU3HEHHOrO Waka MT-mposyKta, BKIOUad TecTHpoBaHie. OTMeTHM, 4TO KOOpAH- 
Hal[HA PeJIM30B (BBIILYCKOB) lIporpaMMHOoro ObecrieyveHHA TECHO CBA3AaHA C YIIPaBJICHHeM KaYeCTBOM, OJHAKO JaHHBbIi TpoLecc 
peyzKo wi (bparMeHtTapHo paccMaTpuBaeTca B JIMTepatype. K Tomy 2Ke He yYHTbIBAeTCA B3aHMOAeHCTBHE TIPOWeCcoB, HOITOMY 
HeT KOMIVICKCHOTO TIpey{cTaBIeHHA OO yipaBJIeHHH KayeCTBOM IPH cO31aHH, TeCTHpOBaHHH U WOpadoTKe mporpaMMHoro 
odecneyenua (I1O). JJanHoe uccieqoBaHie IpH3BaHO BOCTIOJIHHTb yKa3aHHble mpoOespl. Ero esb — mpeCTaBUTb KOMIUICKC- 
HBIM MODXO], CBASbIBAIOLHM TCOpHIo, IIpakKTHKy H MeTOAbI ypaBsIeHHuaA KayecTBOM IO. 

Mamepuaaoi u memooot. Uccnexosana, ipoanasM3upoBana H OTpedbepupoBaHa NpopusbHad TeopeTHueckad UH MIpHKay- 
Has IMTepatypa. 3afelicTBOBaH Mpo*eccHOHabHbI OBIT aBTOpa B yIIpaBsIeHHH KayecTBOM MT-npogzyKTos. YuTeusi 
TpakTHKH PIOOaIbHBIX MOCTABLIMKOB WHpoBbIx TOBapoB UM ycuyr. ABTOP HCHOMb30BaIl ITH MaTepHasIbl HM MCTOBI JIA 
WeTaIbHOu mpopadoTKH Borpocos TecTupoBaHua IIO u pa3sBepTEIBaHHA Kosa. 

Pe3yibmamot uccredosanua. CpopMuposana, onucana u MIpeycTaBJIeHa B BUC CXeMBI KOMIVIeKCHad MOJeuIb yIpaBe- 
HHA KaYeCTBOM Tipu co3qaHHH IO. BaraBseHbl ee B3AMMOCBA3H C MOJCJIbIO MCHEKMEHTA MIpOCKTOB H 2KH3HEHHBIM LMK- 
JIOM IIpOAyKTa, a MMCHHO: aHasI30M, JM3aiHOM, paspaOoTKoH, TecTHpOBaHHeM, pa3BepTbIBAaHHeM H MOAWepxKKo. Yka- 
3aHbI UPHHIHMbI yupaBIeHHA KayecTBOM Ha KaxkKOM W3 9THX cTaqMi. CucTeMaTH3HpOBaHbI H IpesCTaBJIeHbI B BUC 
CX€MBI MPOLeccbI MH MpoOBepKH IPH pasBepTEIBaHHH Kosa. [loKa3aHbI HX OCOOCHHOCTH B Tpex cpeax: pH paspaooTKe, 
TECTHPOBAaHHH H MpOu3BOCTBE. 

O6cyatcoenue u 3ak0uenue. ANTOPUTM NO3BOsIAeT CHeWMasIMCTaM 10 KaYeCTBY BbICTPOHMTb MOCHeAOBATeIbHOCTh Jeli- 
CTBHH JI UCKIIOUCHHA B OYYUIeM BbIABJICHHBIX je(eKTOB, MOHHMaHHA CHTyalMH, KOra MOXKHO (WIM HeIb3A) pasBep- 
TbIBATh KOA H ONpeseseHHA MOMEHTa, Korya cuezyeT nepegzaTs ITO nomp30Batemo. Kpome Toro, ipeqioKeHHast Cxema 
MOxeT ObITb 6a301i JIA ABTOMATH3aLHH pa3sBepTbIBAaHHA KOZa. PewteHHe MO3BOJIMT COKpaTHTb BpeMa Ha pa3paboTKy. Kak 
CJIeACTBUe, MIpOAYKT ObicTpee BbIiieT Ha PbIHOK, YTO YCKOPHT OKyIaeMOCTsb 3aTpaT. BHeszpeHue B MpOH3BOACTBeHHY!O 
mpaxtTuxy MT-kommaHuit Moje, CO3{aHHOH B paMKax JaHHOl Hay4Hon paoorsl, UpesuoaraeT crpaTerMueckue H3Me- 
Henna. Ux pearm3auua TpeOyeT 3Ha4NTeJIbHBIX 3aTpaT BPCMeHH UM Apyrux pecypcos, nosromy oOOMIMi Mpoecc TpaHc- 
cbopMaluit crefzyeT pa3OuTs Ha yacTH. IIpenoxKeHHbI MOAXO aqantTupyetca NOX WYK AbI pasIH4HbIx OpraHv3zalHi u 
TmpozyKToB. MoxkHo paobotatTp C OTACIbHBIMM KOMIOHCHTaMH, YTOOKI CO3LaTb ONTHMAJIbHbIM MW1aH JIA OCTHKeHHA Le- 
lei 10 yrIpaBJICHHIO KayecTBOM. 


Karouesble c10Ba: pa3BepTbIBaHHe Koja, yipaBlieHue KayecTBoM B MT-cdbepe, xu3HeHHEIM “HK MT-npogyxta, 
TecTupoBanue MT-npogykta, pesm3 mporpaMMHoro obecrieueHHA 


Baarojapnoctu. ABtop BbIpaxkaeT OaroyapHOCTb BeAYIIMM opraHH3alluiM B Ciepe pa3spadoTKH MporpaMMHoro 
oOecneueHua — EPAM Systems, Amazon u Google 3a ux WpakTHKH, KOTOpbIe aHaIM3HpOBaJINCb IPH MOATOTOBKe 
cTaTbu. OcoOaa OnaroyqapHocTb KoMnaHHu Project Management Institute 3a oOHOBNeHHe MOAXOAOB K yipaBJIeHH1O 
TIpoeKTaMH. 
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Research (Rostov-on-Don). 2024;24(3):255—263. https://doi.org/10.23947/2687-1653-2024-24-3-255-263 


Introduction. Within the conditions of high competition on the software market, users focus not only on the marketing 
but also on the product utility, pay attention to the user-friendly interface, efficiency, and operational stability. All this 
should be taken into account by manufacturers and divisions of software companies that work with quality. There are 
enough theoretical and applied publications devoted to quality management and automation of these processes in open 
access. The authors consider innovations in this sphere, compare them with traditional practices. 

Defect handling (specifically, their early detection and analysis) allows preventing errors at the software development 
stage. This approach reduces time, money and other resources spent on fixing logical, syntactic, compilation and other 
software errors. The quality of the product is improved, and its value in terms of reliability, ease of maintenance, and 
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cost-effectiveness is increased [1]. In [2], quality management in complex processes is described. The author shows how 
to identify the relationship between technology and quality. According to his assumption, there is an information 
correspondence between probabilistic models of technology and quality. It is known from [3] that in the information 
technology (IT) sector, a 49% increase in budget leads to the creation of new quality management models. As an example, 
we can mention the ISO/IEC standards [4]. In [5], it is shown how quality management affects manufacturing processes. 

An important and complex task is to automate the control of processes on which quality depends. The key form of 
control is statistical. It works with indicators that are critical to quality of the endproduct. These indicators are monitored 
and compared to target values. As a result, a conclusion about the statistical controllability or uncontrollability of the 
process is obtained [6]. 

Software development lifecycle analysis is covered in numerous sources, but there is no detailed description of 
managing each stage of the lifecycle, including testing [7]. 

A separate line of research is the relationship between the basic quality management tools and their impact on the 
operational, financial and market performance of the manufacturing company. At the same time, integrated use and 
management of such tools is an urgent task [8]. 

In [9], it is shown how different quality management practices affect innovation in ISO 9001 certified companies. 
Some studies focus on the incorporation of practice guidelines in product life cycle development. It is known from [10] 
that 60-80% of potential errors in innovative product development are related to misunderstood requirements. The authors 
of [11] argue that when discussing product quality, the issues of code quality, cost of achieving a given quality, and time 
to enter the retail market should be clarified. Continuous integration (or continuous delivery) of program code plays a 
special role in quality automation [12]. This is the process that provides continuous updating of the code during the 
development of a software product. In [13], recommendations are given on how to take into account the level of 
competition in the industry when releasing new versions of a product sequentially. 

Each of these papers considers different elements of quality management. At the same time, the system and structural 
interprocess communication is overlooked. As a result, there is no comprehensive view of quality management in software 
development, testing and revision. This study is aimed at filling this gap. It describes an integrated approach linking 
theory, practice and methods of quality management in the IT project. It should also be noted that due to the flexibility of 
the proposed model it can be adapted to the needs of a particular manufacturer. 

Role of Quality Management in Software Development. Quality is a complex category that can be considered from 
different standpoints: philosophical, social, engineering, economic, legal [1]. In this article, quality is discussed as a set 
of properties of a product or service that can, to one degree or another, satisfy the needs of the target audience [3]. The 
product is considered as the total result of different types of activities, and each has its own input data (parameters), which, 
upon the production process, turn into output (value). The product is a set of values. It is designed to meet needs that are 
determined in advance, during the marketing development of a product or service. The quality of the product is judged 
by how much the actual consumer satisfaction coincides with the planned one. 

In the information technology, two groups of product requirements (user needs) are considered the basic ones: 
functional and non-functional. As on any market, in the retail sale of IT solutions, precisely structured work with value 
provides loyalty of the target audience, and therefore, increases the profitability of the release of software, applications 
and similar products. 

Defective goods sold spoil the image of the manufacturer [8]. Negative feedback is spread online and in messengers. 
The loss of users leads to a deterioration in financial indicators. A quality crisis can result in the bankruptcy of the 
manufacturer. 

In [6], it is shown that quality management through process management is associated with all five types of innovation. 
Let us list the factors that form them: 

— new equipment and technologies; 

— products with new properties; 

— new raw materials; 

— new production organization; 

— new market channels. 

Materials and Methods. The scientific research, whose results are summarized in the presented article, was based on 
the global practice of production and sale of software products and services. In addition, the literature devoted to the 
creation and promotion of IT solutions was studied. The theory has been correlated with the practice of the largest software 
developers. The theoretical part gives an idea of the specifics of different quality management methods. Examples of their 
practical use in product life cycle management are given. From this point of view, some approaches to project management 
are considered — flexible, cascade (second name — waterfall), and hybrid. 


Information Technology, Computer Science and Management 


257 


https://vestnik-donstu.ru 


258 


Birulia MD. Quality Management in Software Development 


The organization and coordination of software releases is studied. It is closely related to quality management, but is 
rarely or fragmentarily considered in the literature. 

The practice of implementing large projects is related to a well-studied theory of quality management. This allows us 
to determine precisely how large organizations use the theoretical basis to provide quality of the information products 
(and ultimately, their commercial success). 

It should be noted that the author of this article is experienced in software products management, and his professional 
expertise was also used as research materials. 

Types of software testing. As shown in [7], on the market (i.e., outside the manufacturing company), the perception 
of quality (as an element of competitive advantage) is determined in the process of product design, statistical control and 
feedback. The perception of quality within the company depends on what percentage of products pass all tests and 
ultimately do not require revision. This indicator is associated primarily with process management, and not with statistical 
control and feedback. 

It is important that the elements of the system that affect quality are supported by senior management and aligned with 
human resource management. In addition, quality practices and indicators should be taken into account when setting up 
internal corporate relations, interaction with suppliers and other contractors. 

Testing is performed to determine the quality of future and end products; therefore, it is important to understand what 
exactly is the object of tests. In this case, it is not enough to say: “product”. This is too general a concept, it should be 
specified. Often, certain functional and non-functional requirements are imposed on the product. In the first case, we are 
talking about a set of functions that the system must perform in a certain way. Let us say, when typing an address, the 
corresponding page should appear in the browser bar. Selecting certain menu items on this page should lead to known 
results in advance. At the same time, there may be several website interaction scenarios — for example, for registered 
and unregistered users. Thus, most online resources, when registering, open up the opportunity to comment on posts, 
receive collections of materials, etc. 

As the name implies, non-functional requirements are not related to the functions available to the user. However, they 
often provide the consumer with more significant benefits than functional ones. This could be a certain level of data 
protection, integrated analytics collection, compliance with legislation (e.g., on personal data protection). 

Product compliance with functional and non-functional requirements is verified differently. 

Figure | shows the functional requirements testing pyramid. Test types are ranked by importance, speed, cost, and success rate. 


Time Cost 


Acceptance 
testing 


System 
testing 


Integration 
testing 


Unit testing 


n> 
Number of tests 


Fig. 1. Types of testing [14] 


Unit testing, or modular testing, is a check of the correctness of individual modules of the program source code, 
individual functions of the application or service. This establishes whether the code is operational, whether it gives correct 
results with different inputs, and whether the logic matches the expected one. An example of such a test: when adding 
two numbers, their sum should be obtained — invariably. The input values required for testing are entered into the 
function. The unit test is passed if the expected result of the function is equal to the one actually obtained. The metric is 
calculated by comparing the number of lines of code that passed the test to the total number of lines of code in the 
repository. The best indicator is considered to be 80%. It shows that the logic of any function matches the expected one 
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by 80% even before the next stage of testing starts. This means that problems are identified quite early and solved much 
cheaper and faster. Therefore, a unit test is the first step in the functional testing pyramid. Its absence (like any other type 
of testing) creates gaps, whose accumulation causes an unjustified increase in the cost of the project [15]. 

Integration testing. The next stage is testing the interaction of the components of the future application. The 
integration of modules is checked. These tests are more complex than unit testing, since the ability to perform sequential 
actions in various components is checked. Example: the user logs in — is redirected to the main page — places an ad. 
That is, several modules are tested at once. The basic options for integration testing are listed below. 

“Big Bang”. All components are assembled and tested together. This allows testing to be performed once after 
development. However, with a large number of modules, some may be overlooked, and this is the weak point of the 
method. In addition, the feedback loop increases, since ready-made solutions are tested when the development is 
complete. We should also note the difficulties with error localization. They are explained by the fact that it is necessary 
to analyze the entire development process (scenario) to find out the cause of the problem. This means that testing time 
will increase. The method is definitely convenient for small systems. 

Incremental approach. Unlike the “big bang”, this testing can start even if only two modules are ready for it, and 
add the rest as they are developed. In this case, it is easier to localize the initial errors, the feedback loop is shortened. 
There is no need to wait until all services are ready. This means that the testing will start and finish earlier. Removing 
defects will be cheaper. Incremental testing can be done in two ways: bottom-up and top-down. 

Bottom-up integration. The least critical modules are tested first, then — the more important ones, and finally — the basic 
ones. Therefore, the upper components cannot be tested until there is data on the quality of the lower-level elements. This 
increases the time between the development and the start of testing. 

Top-down integration. Testing starts with the most significant application layers. In this case, the testing of lower- 
level modules may not be completely adequate. However, even if this turns into problems later, they will be uncritical. 

Hybrid integration. Modules of different levels are tested together. The success of such testing is determined by 
whether the developer understands the application architecture and can determine the optimal approach for a specific 
product [12]. 

System testing. All components of the program are tested as a single application. The specialist must ensure that the 
product correctly handles various scenarios and situations. This option may include testing non-functional requirements. 
This approach is based on requirements or on use cases. Whatever choices are made, the result will be test cases. Their 
successful completion confirms that the system works exactly as expected. 

System testing requires considerable time, because as the product evolves, the number of possible use cases grows. 
Therefore, the software needs to be updated in a timely manner. This is done by the support service. With each product 
update, it is necessary to check whether the previously tested functions work. This process is called regression testing. If 
the tests worked before, but stopped working after the introduction of a new function, then there is a problem. The 
procedure of such checking sometimes takes up a huge amount of operation-use time for testing engineers, and they do 
not have time to write new tests. This explains the importance of automating tests that are performed manually. With this 
approach, testing is optimized to a simple analysis of the results of automatic checks. The best option is 100% automation. 
It is not always possible to achieve this figure, since some tests cannot be automated due to their complexity [16]. 

Acceptance testing. At the last stage, the product is tested by key specialists and customers. They check in real time 
to what extent their expectations and requirements have been met. It is determined if the functions work as intended, the 
interface is convenient, how likely errors and incorrect operations are, what problems the future user will face. As a result, 
the project is accepted or not accepted. Acceptance is the final stage of testing. 

Fixing any bugs or defects found at this stage will be especially costly, since solving the problem will require going 
through the entire process from the very beginning, from development and unit testing to full system verification. Therefore, 
every effort should be made to avoid this situation at earlier stages. Non-functional testing checks the application for: 

— performance; 

— reliability; 

— compatibility; 

— security; 

— utility; 

— scalability. 

More than 15 types of testing are used to check these properties. The most common ones are listed below. 

Performance testing. The speed and efficiency of an application or system is assessed. The focus is on the speed of 
response to a user request. As an example, testing may include evaluating the loading time of the main page of a website. 
Performance targets are defined, and actual results are compared to them. 
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Load testing. This determines the maximum load a system can handle without failure or significant degradation in 
performance. Specifically, for a website, this means determining how many users can interact with the platform 
simultaneously before crashes or failures occur. 

Fault tolerance testing. The ability of a system or application to remain operational during failures or unexpected situations is 
tested. The purpose of such testing is to detect vulnerabilities and develop measures to maintain continuous operation. 

Compatibility testing. Compatibility of an application or system with other software, with different operating 
systems, browsers and devices is checked. 

Security testing. The level of protection of a system or application from attacks and potential security threats is assessed. It 
is determined whether there are vulnerabilities. Protective mechanisms and the effectiveness of security measures are analyzed. 

Disaster recovery testing. The ability of a system or application to recover from failures, crashes, or other 
emergencies is determined. The goal is to verify the effectiveness of recovery procedures and minimize downtime in the 
event of a disaster. 

Research Results. Within the framework of the presented research, a comprehensive model of quality management 
in software development was formed. Its interrelations with the development life cycle, releases, costs, and the project 
management model (agile, waterfall or hybrid) were described. 

Quality management process in the software development life cycle. Standard digital product development 
consists of six production processes, or stages. Each of them should contribute to quality assurance (Fig. 2). 


Fig. 2. Software Development Life Cycle [17] 


Analysis and determination of product requirements. Problems are investigated, solutions are proposed taking into 
account the stated requirements. This process sets the key indicators for checking the endproduct. Depending on the 
development model, the strategy, approaches and testing tools are described and agreed upon. One of the documents is 
called “Product Vision and Capabilities”. It records the identified problems, notes the links with the business plan, 
describes the market situation, lists the main functional and non-functional requirements. The document is approved by 
the partner parties (i.e., representatives of the contractor and the customer). If necessary, adjustments are made to the text. 
At this stage, it is important to involve security and application architecture specialists who will help to accurately define 
non-functional requirements that are important for future testing. 

Design of architectural and visual solution. Taking into account the requirements identified at the previous stage, 
an architectural solution is developed to provide their implementation, and the design of the user interface. All this is a 
guideline for the development team, who will test and, if required, adjust the product. 

Development. Code is created to solve the tasks and problems of future users. Unit tests are written to check the logic 
of individual functions. The main goal is to achieve the required code coverage rate (at least 80%). However, other 
practices are also widely used, such as code reviews. In this case, the code is reviewed by other team members, and errors, 
shortcomings, omissions, and vulnerabilities are identified early in the development process. Reviewers can leave 
comments on specific lines of code, indicating the need for refactoring. In terms of product quality, the ideal scenario is 
dual code review, where at least two specialists unrelated to the author review and comment on the code. 

In addition to unit testing and code review, it is necessary to use linters. These are automated tools that analyze code 
in real time, detect shortcomings, errors, and suggest ways to fix them. Other types of testing can be used depending on 
the strategy and requirements for the product. 

Testing. Test engineers develop test scripts that are adequate to the tasks. As a rule, they are written manually, which 
slows down the implementation of this stage. It is advisable to automate the process, but this is not always possible. At 
this stage, specialists check and, if necessary, re-check the results of manual and automated tests. Particular attention is 
paid to cases where automated tests detect errors. The detected defects are ranked according to their importance to the 
system or users, and then are passed on to the developers for correction. After the defects are fixed, it is useful to perform 
a root-cause analysis (RCA) of the problem using the RCA method. This provides developing measures to prevent 
problems in the future, add new unit or integration tests. To fully understand the causes of the problem, it is recommended 
to use the “5 Why’s” technique. The method scheme looks like this: the problem is formulated and the question is asked: 
“Why did this happen?”. It is followed by an answer describing a certain fact or situation. The question is asked again: 
“Why did this happen?” And this is repeated five times. Five answers allow you to get closer to understanding the causes 
of the problem. 
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Testing is a large-scale process. The major stages are listed below. 
. Setting up and configuring the testing environment. 

. Writing and executing manual test scripts. 

. Automating previously written manual test scripts. 

. Checking the results of automated tests. 

. Maintaining up-to-date automated test scripts and fixing errors. 
. Recording defects. 

. Documenting defects. 

. Analyzing the causes of defects. 

. Developing steps to prevent future defects. 

Deploying the application code. Software requires various environments for deployment. It is important to have a 
good understanding of how continuous integration and continuous deployment (CI/CD) maintain stable software quality 
with minimal effort. This is facilitated primarily by the automation of code deployment — the transfer of its changes from 
one environment to another. The process is divided into stages that depend on the product, organization, and accepted 
standards. 

The degree of automation, the number and rigor of checks determine the expected quality of the product before the 
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code is released to end users. Figure 3 shows a possible code deployment scheme. With a high level of automation, labor 
costs will be minimal, and only a few manual checks will have to be performed. 


Development Testing Production 
environment environment environment 
Development Performing 
automatic checks 
Fixing Unit tests (80%) Conducting Unit tests (80%) Defects No defects 
defects automated tests found found 


100% pass Analysis 100% pass 
Analysis of causes of integration tests of results of integration tests 
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latest version 
of code 


of defect 
No defects 
Vithewbilivehee || Vulnerability 
i i check Fixing defects 
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the process 


Fig. 3. Code Deployment Process and Checklist Diagram 
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Following the presented scheme, you can control the quality and automate the deployment of the code. This will 
reduce the development time, which means the product will enter the market faster. 

Product support. At this stage, it is necessary to build feedback with users: conduct surveys, collect and analyze 
metrics. All this will be the basis for tasks and ideas for improving the software. Assumptions should be confirmed by 
A/B tests, i.e., comparison of product options. It is also important to provide feedback to the support team, which 
directly contacts users. These are the specialists who share valuable observations, convey wishes of the target audience, 
on the basis of which analysis is performed and new strategies are formed. Evidently, the process does not end even if 
the product is completely ready. New data appears for analysis, collection of requirements, and development of a 
digital product or service. 

Discussion and Conclusion. The quality management model presented in this article can be adapted to the needs of 
various organizations and products. Its components can be further studied to create an optimal plan for achieving quality 
management goals. However, it is important to keep in mind the key factors for commercial software success: developers 
and testers must think critically, establish collaboration, and continuously improve processes. 
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The proposed model involves strategic changes, whose implementation requires a significant amount of time; 
therefore, the overall transformation process should be broken down into parts. 

The presented material will be useful for management that administrates quality issues related to the company's overall 
strategy and optimization of individual processes. 

It should be noted that achieving significant results requires a high degree of testing automation, a developed 
engineering culture, and significant costs for creating and maintaining infrastructure. 
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Abstract 

Introduction. Numerous researchers of the reliability of building structures pay attention to hardness, an important 
characteristic of the structural material. It is determined by indentation — pressing the tip of the tool into the surface. The 
advantages of dynamic indentation methods and the distribution of stress intensity on the surface and inside the sample 
are investigated. However, the condition of layered materials on impact has been poorly studied. The objective of the 
presented work is to consider indentation for a two-layer sample and determine the sensitivity of the top layer to the 
strength of the substrate. This will allow us to identify significant characteristics of the strength properties of 
homogeneous and heterogeneous structures. 

Materials and Methods. An elastoplastic model of material behavior and a shock indentation scheme were used, which 
took into account the masses of the indenter and the striker coupled by linear springs. The surface of the indenter was 
conical, the opening angle was 120°. The impact was simulated in the MATLAB system. Finite element model in Ansys 
APDL was used to verify the data and analyze the results of the experiment. Traditional models of elasticity theory were 
used for calculations. The behavior of the material in the zone of plastic deformation was described using the options of 
multilinear isotropic hardening and the von Mises plasticity criterion. 

Results. The results of comparing three versions of varying the level of yield strength in the bottom layer are presented: 
when the yield strength in the bottom layer is half as high as the top one, equal to it, and twice as high. Displacements at 
different observation points for samples with a top layer of 2 mm and 1 mm were analyzed. In the first case, under horizontal 
shear, the displacement indices inside the sample did not change if the yield strength level was twice lower or higher than in 
the top one. If these indicators were equal, the difference became noticeable. In the second case (layer 1 mm), the difference 
in displacement was visible at all observation points. Thus, it can be reasonably concluded that a structure with a smaller top 
layer is more sensitive to impact. In the course of the research, it became known that vibrations associated with the transition 
to the plasticity zone occurred in the 2 mm zone, and elastic damping vibrations occurred below this zone. We solved the 
classification problem for the top layer of the material with changing characteristics of the base. The indicator for comparison 
was the Brinell hardness (HB) in the range of 200-600. The results were processed using a neural network and visualized in 
the form of graphs. The accuracy of its calculations was 98%. 

Discussion and Conclusion. To determine the strength properties of homogeneous structures, it is sufficient to 
characterize the speed of displacement inside the sample. For an inhomogeneous structure, additional parameters should 
be introduced — displacements on the surface and inside the sample at fixed observation points. An integrated approach 
to determining the strength properties of an inhomogeneous structure improves the accuracy of calculations, and the use 
of neural networks increases their speed. 


Keywords: multilayer structure, layered material on impact, yield strength level, Brinell hardness, strength of 
heterogeneous structure 
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AHHOTalna 

Beedenue. Muorne uccieqopatesM Hae%KHOCTH CTPONTeIbHbIX KOHCTPyKUMH yeAIOT BHHMaHHe TBEPAOCTH — BarK- 
HOW XapaKTepHCTHKe KOHCTpyKUMOHHOrTO MatTepuaa. Ee ompeyesiai0T HHJCHTHpOBaHHeM — BJ{aBJIMBaHHeM HaKOHeY- 
HUKa MHCTPyMeHTAa B HoBepxHocts. UccnezyroTca IpeuMylecTBa MeTOOB JHHAaMHYeCKOTO HHCHTHpOBaHHA, paciipe- 
JeIeHWe MHTCHCHBHOCTH Halips.KeHH Ha WOBCPXHOCTH UM BHYTpH OOpa3ita. OWHAKO MaJIO H3y4eHO COCTOAHME CHIOHMCTBIX 
MaTepHasoB Ip yqape. Llenu npeyctaBseHHol padboTbl — paccMoTpeTb HHCHTUpOBaHHe JIA AByXCOMHOTO OOpa3ila 
Wf OMpeeIMTb YYBCTBUTCJIBHOCTb BEPXHETO COA K MPOUHOCTH NOWIOKKU. ITO MO3BOIMT BbLABHTb 3HAYMMbIe XapakTe- 
PHCTHKM NpOUHOCTHBIX CBOHCTB OJHOPOAHBIX H HEOJHOPOAHBIX KOHCTpyKUHH. 

Mamepuanoi u memoooi. Ucnions30Baiu yupyro-liacTH4ecky!0 MOJesIb NOBEACHHA MaTepHasia U CXeMY yJapHOro HH- 
JeCHTMpOBaHHA, KOTOPad YAMTbIBAeT MACCbI HHJeCHTOpa MW yapHuka, CLeIJICHHbIX JIMHeMHEIMM Wpy2KvHamu. IloBepx- 
HOCTb HHCHTOpa — KOHMYeCKaA, yro packpbiTua — 120°. Yoap MoyesMpoBau B cucteme Matlab. KoneyHoomementT- 
Hy!0 Moye B Ansys APDL npuMenvin AIA BepHUKAaLMH WaHHbIX VM aHasI3a pe3yIbTaTOB 9KcIepuMenta. JIA pacue- 
TOB B3AJIM TpayqMUMOHHble Moen Teopuu yipyroctu. Ilopeyqenue MaTepuana B OOacTH MacTH4ecKoro AetopMupoBa- 
HMA ONMCaIM C MOMOM{bIO ONWUMM MYJIbTHIIMHEMHOrO H30TPONHOLO YIpOuHeHHA HW KpuTepua WiacTH4HocTH Mus3eca. 
Pezynbomamoi uccnedosanua. I{puBoqaTca UTOrM COMOCTaBJICHHA TpexX BaPHaHTOB BapbUpOBaHHA YpOBHA Tpeyesa Te- 
KY4eCTH B HWKHEM CyI0e: Kora peel TeKYYeCTH B HWDKHEM C10 BABOe MEHbILe BEpXHErO, paBeH eMy H BABOe OobUe. 
I[poanamm3upoBaHbl NepeMelleHHA B pa3HbIX TOUKAX HaOJOAeHHA JJId OOpa3slOB C BepXHHM cyI0eM 2 MM H | MM. B nep- 
BOM CJlydae Ip TOpH30HTAIbHOM CBHTre He MCHAIOTCA MOKa3aTesIM MepeMeleHHi BHYTPU OOpa3iia, CCIM YpOBeHb Mpe- 
Wella TeKYYeCTH BABOe HWKe VI BbILe, YM B BepxHem. IIpu paBeHcTBe 9TUX NOKa3aTeseH pa3Hula CTAHOBUTCA 3aMeT- 
HOH. Bo BTopom ciryyae (co | MM) pa3HHia MepemMeleHuii BUAHa BO BCex TOUKaX HaOmOAeHUA. Tak MOX%KHO OOO0CHO- 
BaHHO 3aKJIIOUHTb, YTO KOHCTPyKUHA C MCHBIIHM BepXHHM CJI0eM Ooslee YyBCTBUTeJIbHa K yapHoMy Bo3yelicTBHIO. B 
XO H3bICKaHHHM CTaJIO H3BECTHO, YTO B 30He 2 MM COBEpIalOTCA KOeOaHHA, CBAZAHHBIC C TlepeXOJOM B 30Hy IylacTH4- 
HOCTH, HWKe STON 30HbI — ypyrve 3aTyxalomue KoueOaHna. Penn 3ayauy KaccHpukalluH Id BepXHero CIOA Ma- 
TepHasia C MCHAIOLWMMUCA XapakTepHcTHKaMn OcHoBaHHA. Iloka3aTelb JI cpaBHeHuaA — TBepzocTb 0 bpuHemto (HB) 
B Auara30He 200-600. PesymbTaTbl BU3yYaIM3HpoOBasIM B BUAe rpaduKoB u OOpabotaH C MOMOLLbIO Helipocetu. TOYHOCTE 
ee BLIMHCIeHHH CocTaBusia 98 %. 

O6cyorcoenue u 3akmo4uenue, Ja orpeqeueHvA NPOUHOCTHBIX CBOMCTB OJHOPOJHbIX KOHCTpyKUMM JOCTATOUHO xapak- 
TePHCTHKH CKOPOCTH TepeMeleHHA BHYTpH OOpa3zya. Ja HeOMHOpOAHOM CTpyKTypbI HEOOXOAHMO BBOAMTb OMOJHH- 
TCJIbHbIe TapaMeTpbl — TlepeMeljeHHA Ha MOBepXHOCTH MU BHYTpH OOpa3ia B (PMKCHPOBAaHHBIX TOUKAaX HaOMOAeHuH. 
KommieKcuHbiii NOAXO K Ope eeHHIO MIPOUHOCTHBIX CBOMCTB HEOTHOPOAHOM KOHCTPYKI[MH MOBbIIMAeT TOUHOCTS pac- 
YeTOB, a HCIOJIB30BaHMe HelipoceTeli — HX CKOPOCTb. 


Kurouesbie CI0OBa: MHOTOCIOMHad KOHCTpyKUHA, CIOMCTHIM MaTepHal pH yape, ypoBeHb peyeya TexyyecTu, 
TBepAOCTb 10 bpvHesw10, IpO4HOCTb HEOAHOPOAHOM KOHCTpyKWMU 


BaarogxapHocru. ABTODBI BbIpaxKaroT OaroqapHOcTB peak U peweH3eHTaM 3a BHHMAaTeCJIbBHOe OTHOIMICHHE K CTaTbe 
HW 3aMedaHhA, KOTOPbIe NO3SBOJINJIN MOBbICHTh Ce KadeCcTBO. 


Aisa untTupopanua. badyumuna H.E., Jianmn A.A. Onpeyenenue quHamvueckux HallpmKeHHM HW MepemMeljeHui upu 
TelcTBMM yJapHon Harpy3kKH Ha JByXCNIOMHYy!O KOHCTpyKUMIO B Tipoljecce HHAeHTHpoBaHHa. Advanced Engineering Research 


(Rostov-on-Don). 2024;24(3):264—273. https://doi.org/10.23947/2687-1653-2024-24-3-264-273 


Introduction. When designing and operating building structures, strict requirements are specified to their strength 
and reliability [1]. The literature discusses the issues of reliability of the structure [2], as well as the preservation of its 
operational properties throughout its entire service life [3]. 
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From an engineering point of view, hardness is one of the important material specifications. It is associated with 
mechanical properties, such as yield strength, tensile strength, endurance, etc. Hardness of the material is determined by 
pressing the tip of the tool into the surface. There is a distinction between static and dynamic indentation. The static 
indentation test involves the indentation of a solid indenter into a flat and smooth surface of a soft material (target), whose 
mechanical properties are determined by the measured dependence of the applied load on the depth of penetration. Static 
methods do not allow us to evaluate the physicomathematical characteristics of the material under dynamic loading 
conditions [4]. Dynamic indentation methods are used for this purpose [5]. 

GOST R 56474-2015! presents dynamic indentation as the introduction of an indenter into material under the action of a 
single shock pulse created by a special acceleration device or gravity. This approach applies to non-destructive testing methods. 
Their basic principle is safe inspection, determination of the integrity, and basic working properties of the object [6]. The primary 
advantage is the ease of use. 

Indentation testing is a complex process that includes contact mechanics, material nonlinearity, and fracture 
mechanics. For general cases, it is very difficult to obtain analytical solutions. Therefore, the understanding of the 
processes under consideration is based mainly on experiments and finite element modeling. 

Elastic contact was first studied by G. Hertz [7] at the end of the XIX century. He tried to find an accurate definition 
of hardness using an elastic process. Later, J.V. Boussinesq developed a method based on potential theory for calculating 
stresses and displacements in an elastic body loaded with any rigid axisymmetric indenter. As for the depressions 
associated with plastic deformation, early research focused on yield strength and tensile strength. 

Publications pay attention to the methods of dynamic indentation to determine the strength properties of structures. 
Thus, N.N. Avtonomov and A.V. Tololo considered the problem of pressing a spherical indenter into an elastoplastic 
material [8]. The authors analyzed the distribution of stress intensity on the surface and inside the sample. Thus, it became 
known that the zones of maximum intensity are located at a short distance from the contact zone of the indenter and the 
sample, and expand with increasing load. 

Foreign authors [9] have developed a method of dynamic indentation, which consists in measuring the depth and 
reaction to the load of the sample during the indentation process. To determine the depth, a displacement measurement 
method based on moire interferometry was used. The load was measured with a quartz sensor. The results of numerical 
simulation by the finite element method indicated that this approach was in good agreement with the values obtained 
using traditional methods for determining the strain rate. 

Scientists considered strength properties of materials in the indentation process. Many papers are devoted to this, for 
example [10]. The study of the behavior of the “target” under the mechanical loads plays a key role in the development 
of new materials, structures and products [11]. Comprehensive research makes it possible to solve the problems of 
optimizing operating procedures, providing the reliability of structures, and preventing failures in their operation. 

The solutions obtained open up the possibility of creating more efficient and sustainable structures [12], which 
determines the relevance of research in this direction. At the same time, the issue of the state of layered structures under 
impact is still insufficiently studied. 

The work is aimed at analyzing stresses and displacements in a two-layer structure and identifying the sensitivity level 
of the top layer of the sample to the strength of the substrate. The study will allow us to establish significant characteristics 
for determining the strength properties of homogeneous and heterogeneous structures. 

Materials and Methods. The diagram of the impact indentation device is shown in Figure 1. 


mi 


F(t) 


AIA AEA TLS S AS 


Fig. 1. Diagram of the impact indentation device 


' GOST R 56474-2015. Space Systems. Non-destructive Testing of Physical and Mechanical Properties of Space Technique's Materials and Coatings 
by Dynamic Indentation. General Requirements. Electronic Fund of Legal and Regulatory Documents. URL: https://docs.cntd.ru/docu- 
ment/1200122009 (accessed: 14.05.2024). 
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The system of differential equations corresponding to the dynamics of the mechanical elements of the shock 
indentation device has the form: 

mx, =k, (x, —x,)-F(t), 

—_ (1) 

MX. = —ky xX, —k, (x5 —x, i 

This is how the oscillatory process of a dynamic system with two masses is described: m; — mass of the indenter; 

mz — mass of the striker. The masses are coupled by linear springs with stiffness k; and ko, respectively [13]. The values 

of vertical displacements x(t), x2(f) as degrees of freedom of massive elements are unknown in time. The striker was 
cocked to height 4. After descending, at the moment of contact with the indenter, he reached the speed 


h2 
Vo9 = 2gh+= : 
2 


Thus, the initial conditions should be added to the system (1): 
x, (0) =0, x, (0) =0, x, (0) =0, x, (0) = va9. (2) 

Value F(¢) of the resistance force from the injected material is unknown. 

We assume that the surface of the indenter is conical with opening angle a = 120 degrees. Therefore, it is required to 
use an elastic-plastic model of the behavior of the test material. At the apex of the cone, the onset of impact and plastic 
deformation coincides. 

The impact process was modeled in the MATLAB system. This made it possible to select the mechanical parameters 
of the installation. A finite element model in the Ansys APDL environment was used to verify the data and analyze the 
results of the experiment. A two-layer construction was considered. Its displacements and stresses arising in a two-layer 
sample under the action of a dynamic load were analyzed. Traditional models of elasticity theory were used in the 
calculations. To describe the behavior of the material in the zone of plastic deformation, the options of multilinear 
isotropic hardening and the von Mises plasticity criterion were used: 


foe, \=6.-4;.=0. (3) 


: ; 1 ieee 
where o, — von Mises equivalent stress, o, = | (a o- Ful (o)’ } o, — uniaxial yield strength. 


To control the convergence accuracy, a finite element grid was pre-selected. When studying the parameters of a two-layer 
structure, we considered three options for varying the level of yield strength in the bottom layer in three cases — when the yield 
strength in the bottom layer is half as high as the top one, equal to it, and twice as high (Table 1). 


Table 1 
Structure and Parameters of the Studied Structures 
No. ; ; : Level, mm 
Level of yield strength in the bottom layer relative to the top one ? 
(strength coefficient K) top bottan 

1 1/2 

2 1 1 9 

3 2 

4 1/2 

5 1 2 8 

6 2 


A detailed analysis of the displacements at different observation points clearly shows the difference in the nature of 
displacements in the zone of plastic deformations. To study the displacements, the following scheme was selected (Fig. 2). 
The top layer of the structure was highlighted in a darker gray color for clarity. 
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I-33: +5 Sample surface 1 3°55 Sample surface 
-1.5mm @ -1.5mm 
2.5mm @ 2.5mm @ 
a) b) 


Fig. 2. Determination of points for measuring vertical and horizontal displacements: a — construction with a top layer of 1 mm; 
b — construction with a top layer of 2 mm. Numbers 1, 3 and 5 indicate the level of yield strength in the bottom layer relative 
to the top one in accordance with Table 1 


Research Results. Figure 3 shows a comparative characteristic of the displacement level on the sample surface with 


a horizontal shear of the observation point from the impact point. 
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Fig. 3. Displacements on the sample surface at various contact points: a, b — at the 1st contact point; 
c, d — at the 3rd point on the surface; e, {, — at the Sth point on the surface 
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Evidently, for a two-layer structure with a top layer of 2 mm, a change in the yield level at the Ist and 3rd contact points 
during horizontal shear does not affect the displacement parameters inside the sample. At the 5th point, the difference 
becomes noticeable. For a structure with a top layer of 1 mm, the difference in the level of displacement is noticeable at all 
observation points. It can be concluded that a structure with a smaller top layer is more sensitive to impact. 

Figure 4 shows the distribution of stress intensity at different levels of the top coating. These data suggest that 
fluctuations associated with the transition to the plasticity zone occur in the 2 mm zone, and elastic vibrations that fade 
over time occur below this zone. 
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Fig. 4. Stress intensity at different values of strength coefficient K, MPa: a— K= 0.5; b— K=1; c— K=2 
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Figure 5 shows graphs of stresses on the surface at various observation points. 
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Fig. 5. Stresses on the sample surface at various contact points: a, b — at the contact point; 
c, d —at the 3rd point on the surface; e, f, — at the 5th point on the surface 


To analyze the influence of the base on the strength properties of the sample, the problem of classification [13] by the 
value of Brinell hardness (HB) of the top layer of the material with varying characteristics of the base was solved. Within 
the framework of the study, five classification groups were identified (Table 2). 


Table 2 
Classification Groups 
No. Group Value of HB 
1 I 200 
2 II 300 
3 I 400 
4 IV 500 
5 Vv 600 
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The calculated values obtained using a mathematical model were fed to the input of a neural network. Figure 6 
shows the values of speeds and displacements for the selected hardness groups. 
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Fig 6. Calculated values for classification groups: 
a — rate of indenter introduction; b — displacements 


For further data processing, a created and trained neural network with a fully-connected structure was used [14]. It is 
worth noting that neural networks have extensive capabilities in solving this type of tasks [15], in particular, in the sphere 
of construction [16]. Neural networks are often used for qualitative assessment, forecasting and monitoring of the 
condition of building structures [17], as well as for parametric identification of objects [18]. The increase in computing 
power of neural networks directly affects the development of their functionality [19]. 

The accuracy of neural network calculations was 98%. Therefore, the neural network correctly analyzes the 
experimental data and is able to adequately take into account the behavior of the material under shock loading. 

Discussion and Conclusion. Thus, the analysis of the distribution of displacements and stresses in the sample allows 
us to draw a number of conclusions. 

1. Level of vertical displacements in the wave process is much higher than in the horizontal one. 

2. Level of stiffness of the substrate affects the distribution of displacements and stresses. 

3. In the more solid bottom layer, the displacement is noticeably less than in the soft one. 

4. Top layer is sensitive to the strength of the substrate if its height does not exceed 1 mm. 

5. With distance from the indentation site, the deformations become elastic, the stress does not go out into the zone of 
plastic deformations. 

Based on the research results, it can be argued that to determine the strength properties of homogeneous structures, 
one value characterizing the speed of displacement inside the sample is sufficient. For an inhomogeneous structure, it is 
required to introduce additional parameters, such as displacements on the surface and inside the sample at fixed 
observation points. An integrated approach to determining the strength properties of an inhomogeneous structure will 
improve the accuracy of calculations. And the use of neural network technologies to solve this type of problem opens up 
the possibility of faster calculations. 
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Introduction. Digital solutions make the operation of the company clearer, reduce staff costs, and provide data security. 
Various aspects of automation, digitalization and cloud technologies are described in the literature. The question is raised 
about the quality of the methodological basis for such transformations. Global and narrow technical approaches are 
presented. As a rule, materials are presented from the perspective of experts in the implementation of digital technologies. 
In this paper, for the first time, the author's scheme is proposed that can allow not only providers, but also their customers 
to navigate the upcoming transition to the cloud. The latter will receive systematic information on how to select a 
contractor and the most economically feasible option of cooperation. 

Materials and Methods. The information model was based on the semantic network as a system of nodes, their characteristics 
and connections. Management of the cloud migration project and the migration itself were visualized. The contraction practice 
between providers and their customers was summarized. The specifics of the tasks of such projects were taken into account. 
The part of the subject area related to the implementation of a cloud service is algorithmic — a step-by-step transition to the 
cloud, a generalized scheme of the process taking into account the hierarchy of elements are presented. 

Results. For the first time, a method of self-preparation of a company for the implementation of cloud solutions is 
proposed. The algorithm systematizes the cloud migration processes. The activities related to goal setting, audit, selection 
of cloud environment and services, calculation of the economic efficiency of the project, planning and implementation of 
migration, technical support and scaling of processes are described. The possibilities of determining the economic 
feasibility of measures for the transition to the cloud are shown. The costs of equipment, data storage and processing, 
software licenses, salaries, information security, etc., are taken into account. The amount received is compared to the 
providers' offers. For the final decision, the costs of infrastructure support are taken into account — by the customer or 
the outsourcer. The best option is selected. As a result, the customer gets the opportunity to work with better profitability 
and scale the project. Feedback is provided, and processes are adjusted, starting with IT reaudit. 

Discussion and Conclusion. The proposed solution will give the customer's management a system view of the execution 
sequence when migrating to the cloud, the issues and tasks to discuss with a potential outsourcer. Providers can use the 
algorithm to typify and unify projects, which can eventually simplify the coordination of the list of services and the 
migration procedure with customers. In this way, the parties can free up significant resources in terms of time, labor and 
other costs. In addition, customers and providers can partially use the described semantic network to develop not only the 
organizational, but also the technical aspect of the project. 


Keywords: digital transformation methodology, step-by-step transition to the cloud, cloud migration, profitability of 
cloud migration 
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AHHOTalna 


Beedenue. Undpossie petiicHia 1O3BONAIOT CeaTb paOoTy KOMMaHHH 4YeTUe, COKLPaTHTb pacxoJ[bI Ha MepcoHall, o0ec- 
TICUHTb Oe30NacHOCTb JaHHBIX. B uTepaType ONMCcaHbI pa3sIM4HBIle aciieKTbI aBTOMATH3AaIMN, WMdpoBu3allMu u oOmay- 
HbIX TEXHOJOrHM. CraBHTCa BOIIpOc O KayecTBe MeTOOOrHYeCcKOH OCHOBbI TakHXx MpeoOpa30BaHul. IIpeyqcTaBeHEl 
TIO00asIbHbI UW y3KHH TexXHuyecku NoOAXo HI. Kak npaBMsi0, MaTepHasIbI MOMAIOTCA C MO3MIMM cleluaMcTOB No BHeL- 
PeHHIO WHdpossix TexHouorui. B qaHHo padoTe BIepBble MpesaraeTca aBTOPCKad CXeMa, KOTOPad MO3BOJMT OpHeH- 
THpOBAaTbCA B IIPeCTOAIMeM Mepexoye B OOaKO He TOJbKO WpoBaliyqepaM, HO MW Ux 3aka3uuKamM. IlocneqHue NoryyatT 
CHCTeMaTH3HpOBaHHble CBe]{eHHA O TOM, Kak BbIOpaTh NOApATINKa UW BapHaHT COTPyAHWYeCcTBa, HawOosee Werecoobpa3- 
HbIM C IKOHOMHYCCKOM TOUKH 3peHHA. 

Mamepuaaci u memooot. UndopmatonHas MOjelIb CTpowsach Ha 0a3e CeMaHTHYeCCKOM CeTH KaK CHCTeMa Y3JIOB, HX Xa- 
pakTepuicTHK Vu cBa3ei. Bu3yalM3upoBaHbl yipaBlieHHe IpoeKTOM MUrpallMv B OOWAYHbIe CepBHCbI MW caMa MUTpallna. 
OdoOmnieHa IpakTHKa 3aKIHOUeHHA JOFOBOPOB Mex Ty MpoBaliepamMn UW Ux 3aKa34uHKaMu. YATeHbI OCOOCHHOCTH Tex3aaHuit 
Takux poekTos. UacTb IpesMeTHOM OOMacTH, Kacarolllaica peaM3alluu OONadHOrO CepBuica, alropHTMMpoBaHa — Tpeyi- 
CTaBJICH MOmMaroBbli Wepexoy B oOaKo, OOOOMIeHHad CXeMa TIpollecca c yUeTOM HepapXHH 3JIEMeHTOB. 

Pe3yivmamot uccaedosanua. Buepsble ipeqI0KeH MeETOA CaMOCTOATEIBHOM MOATOTOBKM KOMMaHHH K BHe{peHHlo 00- 
JIa4HBIX pellieHui. ANTOpHTM CHcTeMaTH3UpyeT Mpoleccbl Murpalluv B OOako. OnMcaHbI MeporpHaATHA, CBA3AHHBIe C 
leslenonaranHuemM, UT-ayqutom, BEIOopoM obma4dHON CpezbI HU CepBucoB, pacueTOM 3KOHOMHYeCKON IdeKTHBHOCTU 
TIpoekTa, WIaHHpOBaHHeM HU peasM3allueli MUTpalluu, TexMosAepxKON HM MaciITaOupoBaHHeM mporeccos. IloKa3aHBI 
BO3MOXKHOCTH OMpedeseHHA IKOHOMMYCCKON WerecooOpa3HOCcTH MepolpHaATHH Mo Mepexoyy B OONaKo. Y4MTHIBAIOTCA 
3aTpaTbI Ha OOopyOBaHHe, xpaHeHue MU OOpaooTKy aHHBIX, JIMIeCH3HH Ha COOT, 3apmaThl, oOecrleyeHue HH(opMaln- 
OHHOM Oe30racHocTH u up. lomyaeHHylo CyMMY CpaBHMBaloT C MpesOxKeHHAMU MpoBaligepos. JI OKOHYATeIbHOTO 
pelleHud IIpHHUuMaroTCA BO BHUMAHHe pacxojbl Ha MOA epxKy UHpacTpyKTypbI — 3aKa3dHKOM MIM ayTcopcepom. BpI- 
OupaeTca ONTHMAIbHbIM BapuaHT. B uTore 3aka34HK MOydaeT BO3MOXHOCTh paootats Cc Jy4Iel peHTaOeIbHOCTHIO U 
MacilTaOupoBatTs mpoext. IIpeqycmotpena oTpadoTkKa oOpaTHOH CBA3H MW KOPpeKTUPOBKa Mpoieccos, HadMHad C MOBTOp- 
Horo MT-aygnta. 

Odécyotcdenue u 3akmouenue. IIpeqnoxKeHHoe pellieHve act MeHeJDKMeHTy 3aKa34HKa CHCTeMHOe IipeycTaBliIeHve O 
TOM, B KaKOW NOcweOBaTeIbHOCTH J[eCHCTBOBATb IIpH MUMrpalluu B OOMaKO, KaKHe BOMpOChI HU 3aqa4H OOcyKAaTb Cc 
MOTCHIMMAaIbHBIM ayTcopcepom. IIposaiiqepbl MOryT IpHMeHUTh asIrOpUTM JIA THIM3alMu, YHU@UKAMN WpoeKTOB, YTO 
B WTOre yIPOCTHT cormacoBaHve Cc KIHeHTaMH MepeyHA ycuyr WM MopayqKa Murparuu. TakuM o0pa30M CTOpOHBI 
BbICBOOOAT 3HAYHTEJIBHBIC PCCYPCHI TO BPCMeHH, TPyOBbIM HM MHbIM 3aTpaTaM. K Tomy 2Ke 3aKa34HKH MW MpoBalijepbl 
MOTYT OTYACTH 3a}eHCTBOBAaTb OMMCaHHYIO CCMaHTHYCCKy10 CeTb, YTOOI OTPAaOOTAaTh He TOKO OpraHH3allMOHHy!0, HO 
MW TEXHMYCCKy!O CTOPOHY Tpoekta. 


Koro4ueBble C10Ba: MeCTOJIOJIOTHA I{M(pOBBIXx TIpeoOpa30BaHnit, TMOWarOBblii mepexoy, B o0mako, MHI pala B oOsako, 
peHTaOeJIBHOCTB MH palit B o0nako 


BaarogqapHocrn. ABTop BbIpaxkaeT OaroapHOcTE pewakWHu UW peleH3eHTaM 3a BHHMAaTeCJIbBHOe OTHOIMCHHE K CTaTbe HU 
3aMedaHHA, KOTOPbIeC MO3SBOJIMJIM MOBbICHTb Ce KAaYeCTBO. 


Aina WaTHpoBanua. Bopouosa A... Tumu3aiua mupoektos mepexoya Ha oOsauHbIe cepBucnl. Advanced Engineering 
Research (Rostov-on-Don). 2024;24(3):274—282. https://doi.org/10.23947/2687-1653-2024-24-3-274-282 


Introduction. The current socio-economic environment creates conditions for automation and digitalization of the 
operation of enterprises and organizations. This is due, in particular, to a shortage of personnel, new requirements for the 
efficiency and sustainability of production and management processes. Various aspects of automation, digitalization and 
cloud technologies are being investigated. In [1], a global approach to the digitalization of entrepreneurial activity is 
considered. Paper [2] shows the need to move to new digital business models. Stepwise digital transformation of an 
enterprise involves the consistent development of different levels, which are described in [3]. In [4], the reconfiguration 
of management processes required for digitalization is investigated. 
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Main stages of digitalization are as follows: 

— preliminary research — goal definition and modeling of business processes; 

— selection and implementation of hardware, software and hardware-software solutions. 

Organizations can carry out these procedures on their own, but in practice, they often turn to IT companies. 

One of the priorities of digital transformation is the development of cloud technologies and services [5]. RF 
President Vladimir Putin set the task to provide the operation of cloud infrastructure in the country and develop 
domestic cloud technologies!. 

Adequate implementation of cloud services requires a high-quality methodological framework. In [6], a methodology 
for decision support in the selection of cloud IT services is proposed. However, the issue of algorithmization of the cloud 
transition has not been sufficiently studied, specifically, for customers of cloud services who do not work in tech sector. 
At the same time, they can plan and implement some migration activities themselves. The presented work fills this gap. 
Its objective is to form a standard structure for the project of cloud transition services in the context of digital 
transformation. The proposed solution will be the basis for the creation and implementation of such projects. In addition, 
it will allow you to collect, process and use analytics on the processes being implemented. 

Materials and Methods. A semantic network is selected to build an information model and display knowledge in the 
subject area. It contains objects (nodes and connections) and relations between them. This visual structure simplifies the 
understanding and interpretation of knowledge. Figure 1 shows the semantic network of project management for the 


Provider company 


part 


transition to cloud services. 


Customer concludes concludes 


part 


coordinates 
Technical 
requirements 


for example 


part 


initializes 


IT project 
has 
has 


for example 


Hardware and has_ (Tasks of the work 
software 
requirements 
| q has as 


Project 
manager 


Characteristics 


for example 


regulates implements 


Cloud 


Service 


for example Duration 


Cost 


for example 
Human resources 


hi 


implements 


Formulation 


Project 
management 
software 
product 


for example 


Ld Level 


: HalipuMep 


implements 
includes 


‘+ Current project 


“LeaderTask” 


includes 
-______ Realized project 


Fig. 1. Semantic network of project management for cloud transition: 


o- — nodes; —z| — characteristics; 


— connections between nodes 


' Putin set the task to support the development of domestic cloud technologies. TASS. (In Russ.) URL: https://tass.ru/ekonomika/16418975 (accessed: 
16.01.2024). 
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The provider of cloud solutions enters into an agreement with the customer. Technical requirements (TR) are drawn 
up, and an IT project is launched, which is managed by an expert from the provider company. TR regulate the content, 
composition and requirements for the service (cloud service). The cloud service has the following requirements: 

— to the hardware and software (infrastructure placement, virtual resources, cloud services); 

— to safety assurance; 

— to technical support. 

The customer's representative and the project manager from the provider determine the part that is outsourced and 
will be implemented within the framework of the IT project. 

The part of the subject area related to the implementation of a cloud service requires more formalization. In this case, 
algorithmization makes it possible to describe a step-by-step cloud transition and show a generalized scheme of this 
process. In the algorithm: 

— main processes and subprocesses of the implementation of cloud services are identified; 

— sequence of transition between them is determined; 

— connections that provide returning to previous stages to make other decisions are indicated. 

To describe the technical specifications of the cloud transition, we have analyzed the best practices known from open 
sources, including scientific literature. The IT project management part is formalized using a semantic network. 

The basic actions in the development of the project are taken into account: 

— identification and analysis of requirements; 

— clarification of the customer 's wishes; 

— preparation of technical specifications. 

Based on the project method, the basic tasks (jobs) that are required for the cloud transition are identified, and their 
hierarchy is built. The WBS structure of the project includes a list of tasks-jobs. Their duration, connections are specified, 
and resources are assigned for execution. When resources are loaded, it is possible to recalculate the duration of jobs. 
Then critical tasks are determined, whose total duration corresponds to the duration of the project as a whole. 

Thus, the major features of the IT project are the following: duration, cost, labor resources. There are special software 
products for managing an IT project. As an example, these are “Jira” (created by the Atlassian company, Australia) [7] 
or the domestic analogue of “LeaderTask” (developer: Organizer LeaderTask, LLC) [8]. 

Note that the WBS structure is being built in the project management software product along with the tasks-jobs. The 
provider needs to unify the tasks-jobs. In a unified system, the provider, among other things, maintains owned project 
documentation. Analytics are also generated here to track the implementation of the project and its operation after going. 

On the subnet “IT project has”, the requirement to unify the formulation and execution of tasks-jobs is identified. A 
possible option is proposed in this article. Cloud services are considered from the point of view of business development, 
increasing its competitiveness [9]. 

Research Results. Thus, professional providers offer ready-made cloud services, configure them and accompany the 
operation. However, the company can independently create, implement a solution and supervise its performance. 

The project management methodology? [10] assumes the following stages: initialization, development, execution, and 
completion [11]. Note that the same tasks can be formulated variously, especially if the projects are realized by different 
managers. Typification of solutions reduces labor costs for the creation of project documentation, enhances its quality, 
and improves analysis. 

Main steps for migrating to the cloud are listed below. 

1. Defining goals. 

2. Selecting a provider. 

3. IT-audit. At this stage, the initial IT infrastructure is evaluated (efficiency, security, etc.). 

4. Selecting the type of cloud environment and cloud services. 

5. Cost-effectiveness analysis of the cloud transition. 

6. Migration planning (a roadmap is drawn up, priorities, stages, deadlines, resources, budget are determined). 

7. Migration (application dependency scheme is drawn up, cloud infrastructure is designed, test run is conducted, data 
is transferred to the cloud, and service is put into operation). 

8. Operation monitoring and technical support. 

9. Scaling. 

The transition to cloud services can be represented as an algorithm (Fig. 2). 


? A Guide to the Project Management Body of Knowledge. PMBOK Guide. Newtown Square, PA: Project Management Institute, Inc.; 2017. 579 p. 
URL: https://prothoughts.co.in/wp-content/uploads/2022/06/a-guide-to-the-project-management-body-of-knowledge-6e.pdf (accessed: 12.04.2024). 
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Fig. 2. Cloud transition algorithm: 


[ — process; {Lf} — subprocess, —————»— transition to stage 


278 


Voronova AG. Typification of Projects for the Transition to Cloud Services 


We comment on the presented scheme. 

1. The goal of cloud transition should be consistent with the company's business strategy. As noted above, the use of 
cloud services makes business more flexible and efficient. 

Tasks in the project: 

— negotiations with stakeholders and customers; 

— formulation of the planned business result. 

2. It is necessary to find out in advance whether a foreign company can be a partner in the project. In some cases, 
cooperation is possible only with domestic providers. Restrictions are related to the scope of the company's activities, 
top-secret information, processing of personal data, etc. In these cases, the provider must have a set of documents: 

— license of the Federal Service for Technical and Export Control (FSTEC) for the technical protection of confidential 
information; 

— FSTEC certificate; 

— license of the Federal Security Service (FSB) for cryptography; 

— FSB license to work with the national security information; 

— FSTEC license for the protection of national security information; 

— certificate of compliance when working with personal data. 

Reference offers from cloud providers are as follows: implementation and maintenance of virtual machines, data 
centers, etc. However, the services of various companies differ in the nuances of settings, payment models, quality of 
service and user support. 

If the partner's goodwill is important to the customer, he contacts an organization that is well known in the market. 

At this stage, first of all, the certificate and license of the provider should be checked. Secondly, it is necessary to 
agree on the general characteristics of the service: 

— hypervisor; 

— reliability category of data centers according to the Uptime Institute standard; 

— technical support terms; 

— service level agreement; 

— test access; 

— payment model. 

The above are the tasks in the project. 

Hypervisor is software that is used by a cloud provider for virtualization. It allows for creating multiuser logically 
independent cloud environments. 

The reliability of the data center infrastructure can meet the Uptime Institute standard [12] and, in some cases, it 
should be at least “Tier IIT’ [13]. 

The service level agreement specifies the areas of responsibility of the customer and the client, and fixes guarantees 
for services. The terms and parameters of the test access are approved. 

Various payment models for cloud services are practiced. The most common is hourly. 

3. The customer's IT infrastructure, applications and their connections with the IT environment are analyzed. IT audit 
allows for determining which functions can be transferred to a cloud provider and which ones can be left in owned 
infrastructure. It is not uncommon for a customer to decide on a full transition to the cloud. 

Tasks in the project include: 

— conclusion of an IT audit contract; 

— acquisition of information about the company's software; 

— gathering information about the technical support of the enterprise; 

— acquisition of information about the topology of the enterprise; 

— building a business process model “as is”; 

— building a business process model “as it will be”; 

— compiling a list of functions that are transferred to the cloud provider; 

— submission of the IT audit report. 

4. A cloud solution is selected. First, you need to decide on the type of cloud environment [14]: 

— JaaS (infrastructure as a service); 

— PaaS (platform as a service); 

— SaaS (software as a service). 
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Types of cloud environment differ in the elements that are transferred to the provider for maintenance (data, network, 
servers, operating systems, software, etc.) [15]. Elements of the environment may remain under the control of the 
customer's services. 

Tasks in the project: 

— select the type of cloud environment; 

— define a set of cloud services. 

5. Under the cloud transition, an enterprise needs to estimate the costs on operating a standard data center. It is about 
the costs of a data storage and processing system, communication equipment, software and hardware licenses, servers, 
salaries of IT specialists of the enterprise, information security, etc. The amount received is compared to the offers of 
cloud providers. They can do all these calculations. 

The following costs are determined: 

— for the support of the infrastructure by the client; 

— for the selected type of cloud environment. 

The above are the tasks of this stage of the project. 

6. The selection of migration method depends on the scale, structure and IT infrastructure of the organization. Full 
migration is often recommended for small and medium-sized businesses, partial migration is practiced by large 
companies. Project deadlines range from a few weeks to a year. When switching to the cloud model, issues of supporting 
the access to data sources, using archiving and recovery tools, are being resolved. 

There are tasks of this part of the project. You should create an application dependency map and determine 
the following: 

— migration type; 

— list of critical services; 

— transfer time. 

7. At the migration stage, much is determined by the individual characteristics of the customer, therefore detailed 
solutions are formed for specific conditions. Thus, for example, specifics of migrating physical infrastructure to a virtual 
environment and transferring existing virtual infrastructure depend on this. 

Tasks in the project: 

— transfer (installation from scratch; “hot” or “cold” backup); 

— test run; 

— application recycling (when required). 

8. Problems are possible when using any resources. Cloud services are no exception. One of the central issues of their 
use is safety. As a rule, it is determined by the quality of interaction between the provider and the customer's employees. 
The basic task is continuous improvement and correction of shortcomings. 

9. The diagram shows the scaling process separately — the development of the solution, its distribution (e.g., to other 
divisions of the customer). It is recommended to systematically work out the feedback. In this way, you can get 
information that will indicate the further direction of the project. In some cases, it is advisable to return to the previous 
steps of the algorithm — up to the third, i.e., IT audit. 

Discussion and Conclusion. The cloud transition allows organizations to defend digital processes more flexibly 
and economically. The introduction of clouds can significantly reduce the cost of hardware, software, licenses, salaries, 
etc. The Russian cloud services market is developing. Numerous providers deliver services. The literature devoted to 
the study of these processes is published. Within the framework of this work, a scheme is visualized and explained. It 
can be the basis for starting a digital transformation of an enterprise. It gives management a systematic idea of what 
sequence to act in, what to pay attention to, which issues to discuss with the provider, and which issues to work out 
with their own specialists. 

Providers can use the proposed algorithm to typify projects, coordinate the list of services with customers and the 
stages of transition to cloud technologies on-the-fly. 

The unification of task complexes for each step of the described algorithm should provide an improvement in the 
quality of project documentation due to the elaboration of organizational solutions for the project. It should be noted that 
both customers and providers can partially take into account the presented material to design and coordinate the technical 
aspect of migration. Thus, organizational and management processes should be taken into account when forming 
approaches to data coordination and control in cloud solutions. In addition, information management technologies in 
cloud solutions that will enable migration between providers should be developed. 
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Abstract 

Introduction. Computer vision is widely used for semantic segmentation of Earth remote sensing (ERS) data. The method 
allows monitoring ecosystems, including aquatic ones. Algorithms that maintain the quality of semantic segmentation of 
ERS images are in demand, specifically, to identify areas with phytoplankton, where water blooms— the cause of 
suffocation — are possible. The objective of the study is to create an algorithm that processes satellite data as input 
information for the formation and checking of mathematical models of hydrodynamics, which are used to monitor the 
state of water bodies. Various algorithms for semantic segmentation are described in the literature. New research focuses 
on enhancing the reliability of recognition — often using neural networks. This approach is modified in the presented 
work. To develop the direction, a new set of information from open sources and synthetic data are proposed. They are 
aimed at improving the generalization ability of the model. For the first time, the contour area of the phytoplankton 
population is compared to the database — and thus the boundary conditions are formed for the implementation of 
mathematical models and the construction of boundary-adaptive grids. 

Materials and Methods. The set of remote sensing images was supplemented with the author's augmentation algorithm 
in Python. Computer vision segmented areas of phytoplankton populations in the images. The U-Net convolutional neural 
network (CNN) was trained on the basis of NVIDIA Tesla T4 computing accelerators. 

Results. To automate the detection of phytoplankton distribution areas, a computer vision algorithm based on the U-Net 
CNN was developed. The model was evaluated by the calculated values of the main quality metrics related to 
segmentation tasks. The following metric values were obtained: Precision = 0.89, Recall = 0.88, F1 = 0.87, Dice = 0.87, 
and IoU = 0.79. Graphical visualization of the results of CNN learning on the training and validation sets showed good 
quality of model learning. This is evidenced by small changes in the loss function at the end of training. The segmentation 
performed by the model turned out to be close to manual marking, which indicated the high quality of the proposed 
solution. The area of the segmented region of the phytoplankton population was calculated by the area of one pixel. The 
result obtained for the original image was 51202.5 (based on information about the number of pixels related to the bloom 
of blue-green algae). The corresponding result of the modeling was 51312. 

Discussion and Conclusion. The study expands theoretical and practical knowledge on the use of convolutional neural 
networks for semantic segmentation of space imagery data. Given the results of the work, it is possible to assess the 
potential for automating the process of semantic segmentation of remote sensing data to determine the boundaries of 
phytoplankton populations using artificial intelligence. The use of the proposed computer vision model to obtain contours 
of water bloom due to phytoplankton will provide for the creation of databases — the basis for environmental monitoring 
of water resources and predictive modeling of hydrobiological processes. 


Keywords: environmental monitoring of water resources, phytoplankton boundaries, water bloom contour, blue-green 
algae bloom, space image data segmentation 
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AHHOTalna 

Beedenue. Komibtotephoe 3peHve WIMpOKO HCHONL3yeTCA JIA CCMaHTHYeCKON CerMeHTalH JAaHHbIX JHCTAaHWMOHHOrO 
30H_upoBanua 3emmu (7133). MeToa Mo3B0MeT KOHTPOJMPOBaTb IKOCHCTeMBI, B TOM 4CcIIe BOAHbIe. BoctpedoBaHbl 
ayITOPHTMBI, OOecrieuMBarolliMe KaYeCTBO CCMaHTHYeCKOM CerMeHTallHu CHUMKOB J{33, B YaCTHOCTH, JIA BBIABIICHHA 00- 
JlacTeli c (PUTOMIAaHKTOHOM, re BO3MO2XHO WBeTeHHe BOJbI — UpH4aHHa 3amopos. [en uccieqoBaHud — co3faHue all- 
ropuTMa, OOpaGaTEIBarollero CIYTHHKOBbIe aHHble KaK BXOJHYIO MH(OpMallM1o Id (POpMUpoBaHuA U BepuduKalNu 
MaTeMaTHYeCKHX MOJeseH TUIPOAMHAaMUKH, IO KOTOPbIM OTCIeKUBACTCA COCTOAHHE BOAHBIX OObEKTOB. B WuTepatType 
OMMCaHbI pa3IM4Hble alITOpHTMbI CeMaHTHYecKol CerMeHTalHu. Hosble UccreqOBaHuA COCpeAOTOUeHEI Ha MOBBIINeCHHH 
HaJ@xKHOCTH pacio3HaBaHua — Yallje C TMOMOMbIO Helipocetei. ITOT MOAXOA COBepMIeHCTByeTCA B NpecTaBseHHoH 
padote. Ja pa3BuTua HallpaBlleHHa MpeyararoTca HOBbIM Habop CBeJeHHM U3 OTKPBITHIX HCTOUHUKOB MU CHHTeTM4eCcKHe 
TaHHBle Wd yryaenuad oOobmarollel clocoOHocTu Moyzenru. BiepBbie oOsacTb KOHTypa (PUTOMIAHKTOHHOM Monys4- 
IMM CpaBHMBaeTca C 6a30i TaHHbIX — HW Tak POPMUpyrOTCA rpaHW4Hble yCOBHA IA peau3alluHW MaTeMaTHYeCKHX MO- 
Teel W MOCTpoeHHA TpaHH4HO-aallTHBHBIX CeTOK. 

Mamepuaaoi u memoovi. Hadop cHumxos J{33 JOMONHMIIM C MOMOLIbIO ABTOPCKOTO ayrMeHTAaL|MOHHOYO asIrOpHTMa Ha 
a3bike Python. KommbroTepHoe 3peHve CerMeHTHpOBalO OOMAaCTH (PUTOMMAHKTOHHBIX MOMyIAWM Ha CHUMKax. CBepTou- 
HY!O HelipoHHyto ceTb (CHC) U-Net o6yaumu Ha 6a3e yckopuTerel BEraucineHuii NVIDIA Tesla T4. 

Pe3yravmamot uccnedoeanua. Jia aBtoMaTH3aluu OOHapyxKeHHA OOacTel pacnpoctpaHeHua (PUTOMIAHKTOHA pa3pa- 
O0TAaH aJITOpHTM KOMIbIOTePpHOLo 3peHHA, OCHOBaHHBIM Ha CHC U-Net. Mogenb oleHuH 10 BEIYHCIICHHBIM 3HadeHHAM 
OCHOBHBIX MeTPHK KauyecTBa, OTHOCAINMXcA K 3aqayamM cermMeHtaruu. IlomyyeHbi cieqyroluue 3HadeHHaA MeTPHK: 
Precision = 0,89, Recall = 0,88, Fl = 0,87, Dice = 0,87 u IoU = 0,79. [paduyeckas Bu3yamn3alua pesyibTaToB oOy4eHHA 
CHC ua o6yuarollieM H BaIMaltwOHHOM HaOopax Moka3asla Xopoliee KayecTBO OOy4eHHA Moe. OO 3TOM CBU eTEIb- 
CTBYIOT MaJIble W3MeHeHHA (YHKIWU MoTepb B KOHIe OOyaeHHA. BrmowHeHHad MOJeJIbIO CerMeHTAalIMA OKa3aslaCcb 
O1M3Ka K py4HOH pa3MeTKe, YTO TOBOPUT O BBICOKOM KayecTBe IpesO*KeHHOrO pemeHna. Ilo nIomWaqM ONHOTO WuKcelsA 
paccuHuTalH WI0lmlaqb CCOrMeHTHpOBaHHOH OOMacTu PUTOMMAHKTOHHON NonyranMn. omryyeHHE pesyibTaT Wa UCxo]- 
Horo w300paxeHua — 51202,5 (10 HHopMallMu O KOMYECTBe MMKCeeH, OTHOCAIIMXCA K IIBeETCHHIO CHHe-3eJICHBIX 
BoNopocien). CooTBeTCTBYIOINMM UTOr MoyemMpoBanua — 51312. 

Ooécyocoenue u 3akniovenue. UccneqoBaHne pacilinpseT TeopeTH4eckHe H MpakTH4eckHe 3HaHHA O IIpHMeHeHHH CBep- 
TOUHBIX HeEMPOHHBIX CeTel JIA CCOMaHTH4eCKOM CerMeHTallHH JaHHbIX KOCMHYCCKHX CHUMKOB. YYTHIBad HTOTH padorst, 
MO2XKHO OICHHTb MOTeHIMas ABTOMATH3allMM Tpollecca CeMaHTH4eCKOM CerMeHTalMN WaHHBIx J[33 Wa onpeneneHua 
TpaHull (PUTOMMAHKTOHHEIX MOMMYyIAWM C MOMOMIbIO HCKYCCTBeCHHOrO MHTeseKTa. IpumMeHeHuve IpeoxKeHHOM MOTeIu 
KOMIbIOTEPHOLO 3peHHA Jd MOJYYCHHA KOHTYPOB I[BCTCHHA BOJEI H3-3a PUTOMMAHKTOHA WO3BOJMT CO3aTb Oa3bl aH- 
HbIX — OCHOBY JIA IKONOIM4YECKOrO MOHMTOPHHTa BOJHBIX PecypcoB H MpOrHOCTHYeCKOrO MOJeIHpOBaHHA THApoOOHo- 
JIOTHYECKHX TIpOLWeccos. 
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Introduction. Automated algorithms for processing information received from satellites are needed in various 
fields of activity. Solving fundamental and applied problems of ecology requires segmentation of regions in accordance 
with the focus of attention of researchers. This optimizes the process of studying and modeling hydrobiological 
processes. An example of such local interest is the bloom of water due to the spread of phytoplankton. The phenomenon 
is important for current and complex monitoring of water resources. It is clearly visible from satellites during remote 
sensing of the Earth (ERS). 

Water bloom affects significantly its quality in surface sources used for domestic water supply systems [1]. The 
reaction of phytoplankton populations in the hydrological environment can reliably assess the general state of the aquatic 
ecosystem [2]. The negative consequences of uncontrolled algae growth are mass death of fish (suffocation), increased 
load on water purification plants [3], and pollution of shores and beaches [4]. 

Systematic measurements at automatic water quality monitors, as well as obtaining data from research expeditions, 
are labor-intensive and expensive activities. An additional source of information on the state of the phytoplankton 
community is modern satellite systems equipped with survey instruments. They allow remote recording of the state of the 
algae biomass, tracking its dynamics in a given time period. 

A significant advantage of satellite data as a tool for monitoring water resources is the possibility of full-scale and 
operational control at any point on Earth. A wide view of the water area, as a rule, gives researchers a significant amount 
of useful information. But, despite the active development of systems based on computer vision algorithms, the problem 
of identifying the contours of regions of interest in remote sensing data has not yet been fully solved. 

Good results are obtained by various algorithms of semantic segmentation on images. With their help, it is possible to 
identify and clarify the boundaries and structure of natural objects. In [5], the efficiency of the LBP method (local binary 
patterns) for recognizing objects consisting of curvilinear contours is shown. LBP provides high edge sharpness and detail 
of Earth satellite sensing data. In [6], it is noted that to increase the reliability of recognition, it is required to combine 
artificial intelligence algorithms and such classical methods of image edge detection as Sobel, Kirsch and Laplace 
operators. In [7], a comprehensive approach is proposed for semantic processing of satellite images of unlimited size 
using U-Net neural network models, which showed an F1-score value from 0.78 to 0.91 when detecting objects. 

Paper [8] provides an overview of intelligent methods for solving the problem of semantic segmentation of data on 
satellite images. The authors conclude that in this case, neural network algorithms are the most effective and productive. 
As an example, a convolutional neural network (CNN) is given, trained on several thousand satellite images of 
Massachusetts (USA). The accuracy of the model was 85.31%. In [9], semantic segmentation, instance segmentation, and 
panoptic segmentation are considered. The advantages of using deep learning methods implemented in the architectures 
of such CNN as SegNet, U-Net, and DeepLab are specified. In [10], automated processing of satellite images is based on 
a combination of the SpaceNet dataset and progress in computer vision which are made possible by deep learning. This 
paper presents five approaches based on improvements to the U-Net and Mask R-Convolutional Neural Networks models. 
The metric values for the best models are as follows: average precision (AP) and average recall (AR) are 0.937 and 0.959, 
respectively. An effective application of CNN for detecting contrails on satellite images is described in [11]. It is proven 
that in large-scale monitoring of contrails with measurement of their impact on climate, the approach based on CNN of 
U-Net architecture demonstrates Fl-score equal to 0.52, with an overall average detection probability of 0.51. 
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The 2023 models Segment Anything (SAM), Language-Segment-Anything (Lang-SAM), and HQ-SAM are of 
particular interest. These are dynamic deep learning tools that can predict object masks from images using input hints. 
Several researchers have already applied this approach to the analysis of aerial photographs and ERS data. The accuracy 


of identifying areas of interest has proven to be high [12]. In [13], the F1- score value reaches 86.5%+4.1%. In the future, 
models of similar architecture with various modifications (Polyp-SAM, Grounding DINO, etc.) will provide for both 
interactive (requiring user intervention) and automatic segmentation. 

Intelligent technologies are increasingly being introduced to process remote sensing data. High accuracy of models is 
noted. Particular attention is paid to methods based on such CNN as SegNet, U-Net, DeepLab in combination with 
classical methods of image preprocessing. A generalized approach to segmentation is actively developing. 

This paper considers the solution to a problem in the ERS data assimilation using computer vision. The application of 
U-Net CNN for segmentation of areas containing phytoplankton populations is shown. The algorithm created by the 
authors provides for the segmentation of regions of interest and calculation of their areas, which is required for further 
analysis when solving problems of hydrodynamics and hydrobiology. 

The following four points describe the scientific novelty of the presented study. 

1. A data set was formed from open sources. 

2. Synthetic data were generated to improve the generalizing ability of the model. For this purpose, the authors’ own 
augmentation algorithm was used to make the model more resistant to noise in practical use [14]. 

3. An intelligent model based on the U-Net CNN architecture was implemented in the high-level Python language. Its 
key hyperparameters were optimized using the Optuna library and checked on a test dataset. 

4. The areas of the found contour containing phytoplankton populations were compared to the existing database. In 
this way, boundary conditions have been formed for the subsequent implementation of mathematical models and the 
construction of boundary-adaptive grids. 

To achieve the set goal, it is required to solve a number of problems: 

— to prepare an ERS database containing regions of interest segments of water bloom; 

— to validate and describe the topology of the U-Ne SNS; 

— to perform data augmentation to create an extended representative set; 

— to implement, optimize, debug and test the CNN of U-Net architecture; 

— to determine values of key metrics of the model quality for segmentation; 

— to calculate the areas of the segmented contour given the scale of the original image. 

The theoretical significance of the study is due to the expansion of ideas about the possibilities of using computer vision 
technology in the field of water resources monitoring. The practical significance consists in the development of an applied cross- 
platform and scalable tool for analyzing remote sensing images to record regions of interest in aquatic ecosystems. 

Materials and Methods. For geospatial analysis, we use open-source software that is often applied to solve 
environmental problems [15]. 

The study is based on current satellite data. The authors focus on the state of water bodies during the bloom of blue- 
green algae. Analysis of this information allows: 

— predicting the volume and distribution of phytoplankton in the water area [16]; 

— checking physical and biological processes that determine the rate of phytoplankton growth and biomass 
accumulation [17]; 

— analyzing climate change based on the forecast of the dynamics of the bloom process [18]; 

— studying in detail the process of CO exchange between a water body and air [19]. 

To automate the process of detecting regions of phytoplankton populations and calculating their areas, it is proposed 
to develop a computer vision algorithm based on the U-Net CNN architecture. 

As a training sample for the deep learning algorithm, 20 space images of water bodies such as the Black, Caspian, 
Azov Seas, etc., were taken. The photos were obtained at different points on the earth's surface. 

The first step was to label the images to transform the information into a format that could be understood by the computer vision 
algorithm for performing the segmentation. There were two common approaches to providing annotations: 

— creating a pixel-level mask; 


— selecting polygon boundaries for the region of interest. 
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We have used the first option, where the pixel-level mask files represent regions of interest for the algorithm. The 
marked masks are files with the extension jpeg or png. The proportions correspond to the image they annotate. Figure | 
shows an example of the original image and its mask, where green indicates land, blue indicates water surface, and red 
indicates the region of the phytoplankton population. 


a) b) 
Fig. 1. Image mapping: a — original image; b — image mask 


To increase the number of images in the data set, we used the authors' augmentation code, supplemented with noise 
effects. When creating the extended data set, we used the following modifications of the original images: 

— rotation by an arbitrary angle; 

— display along the OX and OY axes; 

— cropping; 

— scaling; 

— color correction. 

All changes were made taking into account the noise that may appear on real images obtained through remote sensing, 
and were segmented using the developed algorithm. 

Let us note the advantage of the authors' algorithm for creating additional source data. Under conditions of a limited 
set of real images, the use of artificially created images for training will allow for a more fine-tuning of the developed 
model, optimizing its parameters and making it more resistant to distortion in practical application. 

The U-Net CNN architecture is designed to solve the problem of biomedical data segmentation. The determining 
factor in its selection is the relatively small size of the initial data, with which U-Net shows satisfactory results in practice. 

The U-Net CNN architecture is based on the interaction of convolution layers + pooling, which first reduce the spatial 
resolution of the image (encoder), and then increase it, having previously combined it with the image data and passed it 
through other convolution layers (decoder) (Fig. 2). 
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Fig. 2. Architecture of U-Net SNN 


The convolutional blocks of the decoder and encoder are linked by end-to-end connections, or skip connections. This 
solves the problem of vanishing gradient, which is a challenge for computer vision [20]. In this study, we used the encoder 
from the ResNet—50 neural network, pre-trained on the ImageNet dataset. 

To select the hyperparameters of the U-Net CNN that affect the architecture and training process, the Optuna library 
was used. This made it possible to automate the model tuning to achieve better results. 
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Research Results. Table 1 shows the model parameters specified under training. 


Table 1 
Parameters for Training the U-Net Convolutional Neural Network 
No. Parameter Value 
1 | Number of images in training set 700 
2 | Number of images in validation set 200 
3 | Number of images in test sample 100 
4 | Batch size 10 
5 | Learning rate Ist-4 
6 | Overfitting detector Early stopping 
7 |Solver Adam 
The model was trained using optimization of the Dice loss function (2) based on the Dice coefficient (1). 
2X OY 
DSC = a (1) 
HeeiaeK4 2|x O 4 + Smooth (2) 


|X|+|Y|+ Smooth 

Here, X — a set of pixels defined during the mapping as a scope of a specific class; Y — a set of pixels assigned to a 
specific class according to the conclusions of the developed segmentation model. The Smooth coefficient is used to 
smooth the calculation result in the case when the values X and Y are close to zero. 

The Adam method for stochastic optimization was used to train the model. Early stopping was used as an overfitting 
detector. In machine learning, this is one of the most widely used regularization methods to prevent overfitting. The training 
process was performed on the basis of NVIDIA Tesla T4 computing accelerators, it was implemented in 100 epochs and 
took 55 minutes. 

Figure 3 shows the graph of CNN training on the training and validation sets. The OX axis shows the training epochs, 
and the OY axis shows the values of the loss function. Analyzing the graph, we can conclude that the quality of model 
training is good, since at the end of training on the training sample, small changes in the loss function are observed. 
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Fig. 3. Training U-Net CNN: — on training sample; — on validation sample 


When assessing the quality of segmentation models, the Dice coefficient and the metric of the degree of intersection 
between two bounding rectangles (Intersection over Union — IoU, Jaccard index), determined form the following 
formula, are used: 


loU =——_, (2) 


where X — a set of pixels defined under the mapping as a scope of a concrete class; Y — a set of pixels assigned to a 
concrete class according to the conclusions of the developed segmentation model. 
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Table 2 presents the values of per-pixel precision, recall, Fl-score, Dice coefficient, and IoU. To obtain the final IoU 
value, the weighted average is calculated for the values of this metric for each class. 


Table 2 
Results of Model Quality Assessment on the Test Sample 
Metric Precision Recall Fl Dice IoU 
Average value for test sample 0.89 0.88 0.87 0.87 0.79 


Figure 4 shows the results of the algorithm's work on segmenting regions of water resources, land and phytoplankton 
populations. The results obtained satisfy the tasks of water resource monitoring and have practical value. 


- 


d) 
D 


Fig. 4. Algorithm results for segmenting areas of water resources, land and phytoplankton populations: 
a, d — original image; b, e — manual mapping; c, f— model result 


e) 


The segmentation result in Figures 4 c and 4 fis visually close to manual mapping, which indicates the high quality 
of the model. The area of the segmented region of the phytoplankton population was calculated by estimating the area of 
one pixel. Each image provided has additional metadata indicating the image scale and its resolution. Based on this value, 
the area occupied by each pixel is calculated. In the case considered for Figure 4 a, the final value is 51202.5. This figure 
was obtained according to information on the number of pixels related to blue-green algae blooms from a set of segmented 
images of phytoplankton populations in coastal systems [21]. The calculation result for Figure 4 c is 51312. 
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Discussion and Conclusion. When assessing the state of water resources, computer vision and other machine learning 
algorithms allow specialists to become free from monotonous operations. They are performed by intelligent systems. In 
this case, monitoring can be carried out round-the-clock. The algorithm will adequately predict risks, model the 
development of situation, and support the adoption of operational decisions. Stored and replicated knowledge in the form 
of databases and registers can be used to create long-term sources of information that researchers can use to analyze the 
state of water bodies and build climate models. 

Processing ERS data in the form of semantic contours will provide verifying complex mathematical models through 
refining boundary and initial conditions, increasing the accuracy, speed and reliability of predictive modeling of 
hydrobiological processes. 
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Abstract 


Introduction. To improve the diagnostics of knee joint diseases, it is necessary to enhance the quality of processing 
radiographic images, i.e., to provide experts with more accurate information for pathology analysis. The objective of the 
study is to demonstrate the capabilities of fuzzy logic in improving the algorithm for determining reference lines and knee 
flexion angles. This requires a program that analyzes X-ray images. The methods known today, described in scientific 
and applied literature, are not sufficiently automated. In some cases, orthopedists and surgeons have to manually refine 
images and adjust lines. This gap is filled by the presented work. The algorithm developed by the author is described. It 
does not involve human participation and automatically identifies the lines and angles of knee flexion. Based on the result 
issued by the system, the doctor can, firstly, judge the presence of pathology. Secondly, the information provided by the 
program allows for more accurate planning, performing operations, and prescribing therapy. 

Materials and Methods. Images from two X-ray machines operating in Al-Basel Hospital (Latakia, Syria) were used. 
The Python language was used for the software implementation of the algorithm. The solution was tested on 500 patients 
at Al-Basel Hospital. The results generated by the new system and previous versions of X-ray image processing programs 
were compared. 

Results. An algorithm for constructing reference lines and angles for processing knee joint X-ray images is created, 
described, and implemented in practice. The capabilities of fuzzy logic in automating double threshold detection when 
identifying bone boundaries in images are shown. The operation of an improved Gaussian filter designed for processing 
X-ray images is described. The modified method of knee bone X-ray analysis includes the development of an algorithm 
for automatic detection of structures and anomalies in knee joints, determination and measurement of anatomical 
parameters, assessment of the degree of damage, etc. The method for determining the contour boundaries on radiographs 
combined the Canny detector, the watershed algorithm, and fuzzy logic. The program has been implemented in medical 
practice and shows 98% accuracy, spending less than 20 seconds to process the image. 

Discussion and Conclusion. The new system provides high accuracy, acceptable efficiency, and does not require manual 
correction of images. Experts are now able to identify subtle indicators of disorders. In addition, the new method makes 
it possible to understand complex cases when several factors are combined, indicating potential pathology. Widespread 
implementation of the method will improve the quality of medical services in orthopedics. Scientific research in this 
direction should be continued to expand the set of strategies for the treatment of diseases of the musculoskeletal system. 
It is necessary to create solutions with absolute accuracy, higher processing efficiency, as well as methods suitable for 
analyzing other joints. 


Keywords: reference lines in radiography, knee joint angles in radiography, Canny algorithm, improved Gaussian filter, 
watershed algorithm 
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AHHOTalna 

Beeoenue. Ina yay wena uarHocTuKH 3a0o0eBaHHi KOJICHHOFO CycTaBa HEOOXOJHMO MOBbICHTb KauecTBO OOpaboTKH 
peHTreHorpaduyeckux w300paxkeHuii, T.e. JaTb crielMamMctam Oonee TOUHYIO HHPOpMaLlHIo JIA aHau3a WaTOOrHH. 
Lem uccneqoBaHua — MOKa3aTb BO3MO2KHOCTH He4YeTKOM JIOTHKM B COBEPIICHCTBOBaHHN aJITOpHTMa olpeyeueHua 
ONOPHBIX JIMHHM MU yruioB cruOaHna Kosena. Ja 9TOTO HeOOxOAMMa IporpamMa, KOTOpasd aHallM3upyeT peHTreHOBCKHe 
cHuMKH. Vi3BecTHbIe Ha CeroqjHA MeTOJbI, ONMCaHHbIe B Hay4HOH UM UpukaqHol MMTeparype, HeocTaTOuHO 
aBTOMAaTH3HpoOBaHbl. B page ciyyaeB opTomeyqaM H XMpypraM MIpHxoANTcaA BpydHy!o AOpabaTEBaTb W300paxKeHHaA, 
KOPpeKTHpoBaTb JIMHUH. OTOT UpobeN BOCMOMHAeT IpeycTaBeHHas padota. OnvcaH CO3{aHHbIii aBTOPOM asIrOpHTM, 
KOTOpbIM He lIpeqMosaraeT yuacTHA YeOBeKa, ABTOMATHYeECKH HACHTHPUUMpyeT JIMHMM UM yribl crHOanua KoseHa. Ilo 
Pe3yJILTaTy, BbIJaHHOMYy CHCTeEMOM, Bpad MOXKET, BO-IepBbIX, CYHTb O HaIM4NH NaToOorMu. Bo-BTopbix, cBeqeHHA, 
IIpeOcTaBIAeMbIe IPOrpaMMOH, MO3BOJIAIOT TOUHEE MIaHHpOBaTh, IPOBOAHTb OllepallHu MW Ha3sHayaTb TepaliHio. 
Mamepuaaoi u memoodoi. Uctionp30Banucb CHAMKM JIBYX PCHTTCHOBCKHX alilapaToB, KOTOpble paOoraloT B OObHULe 
Asb-ba3enb (Jlatakus, Cupua). Ja mporpaMMHon peasm3alHv aIropHTMa 3aelcTBOBalIM A3bIK «IIuToH» (Python). 
Pemienue upotectupoBamm Ha 500 nauuentax OonbHuyst Asb-ba3serb. CpaBHHBasIMCb pe3yJIbTaTbl, KOTOpBIe 
creHepHpoBasia HOBad CHCTeMa HM Tipe WecTByIOWMe BepCHH IIporpaMM OOpaboTKH peHTreHOBCKHX CHHMKOB. 
Pesynomamoi uccnedoeanua. Co3aH, ONMCaH M peasiM30BaH Ha MpaKTHKe asIrOpHTM MOCTpOeHHA OMOPHbIX JIMHMM U 
YIIOB Jd OOpaOoTKH PeHTTeHOBCKHX CHHMKOB KOJIeHHoro cycTaBa. IloKa3aHbI BO3MOXKHOCTH HeYeTKOM JIOTHKH B 
aBTOMATH3allMH OOHapyxKeHHA FBOMHOTO Mopora MPH BbIABIIGHHM TpaHH KOCTH Ha H300paxKeHHAX. OrnmcaHa padota 
YCOBepIICHCTBOBaHHOTO rayccoBCKOrO (puJIbTpa, IIpeqHasHayeHHOro AIA OOpadoTKH peHTreHorpamM. 
MogudunnpopaHHblit MeTOy aHasIu3a peHTTeHOBCKHX CHHMKOB KOJICHHBIX KOCTel BKJIIOUAeT pa3spaoorTkKy aIropuTMa JIA 
aBTOMaTHYecKorO OOHapyxKeHHA CTpyKTyp HW aHOMaJIMii B KOJICHHBIX CYyCTaBax, OMpeyeeHia WM M3MepeHHa 
aHaTOMHUYeCKHXx MapaMeTpoBs, OLCHKY CTeMeHH NOBpexkeHHA H T.7. 

Merog onpeyeeHua rpaHl, KOHTypoB Ha peHTreHorpaMMax OObeMHUI TeTeKTop KoHHH, alropHTM Boyopa3yeua 
HeyeTKy1o Joruky. IIporpamma peasIM30BaHa B MeJMIMHCKOM MpakTHKe HM MOKa3bIBaeT TOUHOCTH 98 %, 3aTpayuBad Ha 
oOpadoTKy cHuMKa MeHee 20 ceKyHI. 

O6cyocdenue u 3akmio4uenue. Hopas cucTeMa JaeT BbICOKYIO TOUHOCTS, MIpHeMJIeMyIO OllepaTHBHOCTb HU He TpebyeT 
Py4dHo KOppeKTHpoBKH CHHMKOB. CrelMaiHcTbI MOYYHIM BO3MO2%KHOCTb BbIABUTb MaJIO3AMeTHbIC MHAMKATOpbl 
HapyuieHui. Kpome Toro, HOBbIi MeTOJ MO3BONAeT pas0OpaTbesA B CIIOMHBIX CyYaAX, KOTa COUCTAIOTCA HECKOJIbKO 
(pakTOpOB, yKa3bIBaIOWIMX Ha BO3MOXKHY!I0 MaTosoruio. Lupokoe BHeApeHHe MeTOsa MOBBICHT KaYeCTBO MeCAMUMHCKUX 
ycllyr B opToneguu. CreqyeT MpOMOMKHTh Hay4Hble W3bICKAHHA B JAHHOM HalipaBIeHHH JWIA paciuMpeHua HaOopa 
cTpaTerHi slevenua 3a00eBaHHit OMOpHO-yBuraTebHOro altmapata. IIpeqcTouT co3faTb pelieHua c a0coOTHOH 
TOUHOCTBIO, Oosee BLICOKOM OlepaTHBHOCThIO OOpaoorTKH, a TaKWKE METOABI, MOAXOTAMME [IA AHAIIM3a pyrux CyCTaBOB. 


KaroyeBble C10Ba: ONOPHbIe JIMHUM B PeHTreHorpadun, yIibl KOJICHHBIX CYCTaBOB B PeHTTeHorpaun, asIropHT™M 
K9HHH, yCOBepLIeCHCTBOBaHHBIM (puBTp aycca, almropuTM Bosopa3yeua 


Jia waTupopanua. Mannaa A.C. Anroputm oOpadoTKn peHTreHOBCKHX H300paKeHHH C HCTIOsIbB30BaHHeM HeYeTKOM JIOMHKH. 
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Introduction. Improving the quality of X-ray images is a critical task with direct implications for the quality of 
medical diagnosis and subsequent treatment, especially in cases like joint replacement surgeries. Advancements in 
medical imaging technologies, including X-ray imaging, are vital for providing more accurate and detailed information 
to healthcare professionals. There are some strategies and advancements that can help improve X-ray image quality and 
subsequently benefit joint replacement surgeries and other medical procedures. 

Analysis of X-ray images allows identifying and evaluate various pathological conditions [1]. Algorithms for image 
segmentation and bone delimitation have made a significant contribution to improving the clarity and accuracy of X-ray 
image interpretation. However, such visualization still needs improvement, as it is not always completely defined, leaving 
room for different interpretations [2]. 

Advances in intelligent computing make it possible to use fuzzy logic to remove imprecision and ambiguity in the 
medical image analysis. Fuzzy logic mimics human thinking, and it is a reliable basis for creating algorithms that can 
efficiently process and analyze X-ray images [3]. 
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In [4], reference lines and their relationships were studied for the purpose of correcting the shape of the legs. 
In [5], clinical results of treating multiaxial deformities of the lower extremities were compared to how this was reflected 
in X-ray images. The literature describes software products that can analyze X-ray images of bones. However, all these 
systems have a common drawback: sometimes they incorrectly recognize contours. If this happens, the specialist manually 
refines the image. This problem is solved in the presented scientific work. A new algorithm for analyzing X-ray images 
is proposed using fuzzy logic and improved traditional approaches. The solution eliminates the need for manual image 
processing. The program operates with high accuracy and ample efficiency. It allows for more accurate detection and 
characterization of anomalies in X-ray images, i.e., improves diagnostics. 

Materials and Methods. In carrying out this scientific work, the author proceeded from the fact that the assessment 
of limb deformations requires an analysis of the anatomical and mechanical axes of long bones, as well as the angles 
between them. Figure | shows the initial limits of angles for healthy people of different ages. If the indicator is outside 
these limits, then the patient has a deformed limb and needs surgical intervention. 


79-83° 
a) 85—90° 


85—90° 85—90° 


a) b) 


Fig. 1. Support angles between mechanical axes for healthy people: 
a — for young people; b — for elderly people [4] 


In orthopaedics, the concepts of mechanical and anatomical axis are crucial to understanding the arrangement and 
function of bones, especially in the lower extremities. The anatomical axis runs through the center of the bone shaft, 
through the middle of the diaphysis. The mechanical axis reflects the direction of the forces acting on it. In the lower 
extremities, it runs from the center of the hip joint to the center of the ankle joint. 

By examining the mechanical axis, it is possible to understand the load-bearing capacity and load distribution in the 
bone. This is particularly important for the diagnosis and treatment of osteoarthritis, fractures, and for determining the 
position of implants or prostheses. 

It is necessary to know the features of the anatomical axis to distinguish natural curvature from pathological 
curvature. This is also important under surgical interventions, since it helps surgeons understand whether the bone 
is positioned correctly. 

Two parameters are particularly important for assessing the anatomy of the femur and tibia: 

- ALDFA (axial length of the distal femoral axis); 

—- MPTA (medial proximal tibial angle). 

They are extremely informative when it comes to identifying knee deformities. ALDFA determines the orientation of 
the distal femur relative to its longitudinal axis, MPTA shows the inclination of the medial surface of the tibia relative to 
its mechanical axis (Fig. 2). 


Fig. 2. Angles between femur and tibia used to assess anatomy and detect pathology: 
1 — medial genu-diaphyseal angle (MGDA) = 130° (124°—136°); 2 — medial proximal femoral angle (MPFA) = 84° (80°—89°); 
3 — ALDFA = 81° (79°-83°); 4 — MPTA = 87° (85°—90°); 5 — lateral distal tibial angle (LDTA) = 89° (85°-92°) [5] 
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The angles shown in Figure 2, in combination, demonstrate the possibilities of knee joint alignment or deformation. 
Orthopedists use the norm (numbers in brackets) to draw conclusions about the type, severity of deformations, and 
develop a treatment plan [6]. 

The Canny edge detection algorithm is a robust method widely used to identify boundaries in images. It builds 
reference lines and angles from X-ray images as follows. The first step uses an improved Gaussian filter, which smoothes 
the image, reduces noise [7]. Next, gradients are calculated using the Sobel or Prewitt operators. The method is often used 
to detect the image edges. This allows both weak and strong edges to be taken into account. This way, the algorithm 
determines the edges more accurately and reduces the noise level more effectively. 

Despite its computational complexity, the Canny algorithm remains a fundamental part of computer vision due to its 
ability to cope with noise and clearly detect edges. 

To correctly select the standard deviation value, the absolute values of the point x(i, /) in the gray pixel are compared [8]. 
If kG | < eG) < beG, 7) — ul, then the standard deviation o = 1. In this case, the point is less affected by noise. Standard 
deviation o = 1.6 indicates significant noise. 

To automatically detect the edges (71, 72), Canny uses double thresholding, classifying pixels in the image by gray 
intensity. We calculate the average value of all pixels (71), after which 71 is used to divide the image into two 
areas — object and background. Pixel values that exceed the initial threshold 71 are taken as objects, and the rest are 
taken as background. Now there are two sets of values: one for the object (ob) and one for the background (bg). For each 
of the two regions (object, background), a threshold value is calculated. Let us calculate 72. We start with the fact that: 

T0 =(Tob+Tbg)/ 2. (1) 


New value 70 is taken as 71, and the process is repeated until 70 ~= 71. Final value 70 is taken as 72. 

The final step involves boundary binding to identify related boundary components in the VNSA (valgus knee shaft 
angle) image. This helps to combine disparate boundaries into larger structures, which improves image perception and 
analysis. It is advisable to integrate fuzzy logic principles and traditional software approaches to edge detection. In this 
case, membership functions are used that assign membership degrees to pixels based on gray intensity values, allowing 
for a detailed representation of edges. 

Unlike the binary approach, the fuzzy method facilitates a gradual transition in the range of membership values from 
0 (no membership at all) to 1 (full membership), capturing the inaccuracies of the boundaries. Fuzzification transforms 
crisp image data into fuzzy sets, taking into account the uncertainty. Defuzzification transforms fuzzy results into crisp 
ones for further analysis. 

A notable advantage of the method is its adaptability to real-world scripts, specifically to those with varying degrees 
of noise and ambiguity. However, careful parameter tuning may be required to achieve optimal performance; therefore, 
finding a balance between sensitivity and specificity in edge detection applications is critical [9]. 

A more detailed representation of edges can be achieved by disambiguating image boundaries. This requires assigning 
membership degrees to pixels. 

Each pixel is processed using fuzzy rules. If any of the fuzzy rules shown in Figure 3 are met, the pixel is processed 
as an edge, otherwise it is not taken into account. This provides highlighting the boundaries of objects that are different 
from noise and other image elements. The importance of this solution has been proven in relation to the analysis of the 
humerus from X-ray images [10]. 


I4=1 i. B=i,M= i, B=1, 
N&&I2GR3&EIS&&16&&ITKE&IB = 0 N1& SIZE KG KAT KEIB = 0 

I5=1 I1l=1,122=1,B=1, iene 
N&&I2G&I3& 814 & &16&&17K&I18 = 0 | V4& &IS&& 16K ITB = 0 

16=1 1I5=1,16=1,17=1, 
N&&I2GR3& E14 &&IS5&&1TK&IB = 0 a N&SI2KKI3 KKK EIB = 0 P| 
I8=1 Il =1,17=1,18=1, 
N&&I2ERB&&I4&&IS5&&16&&17 = 0 mI N&&I2Z&&B&&I4&&I8 = 0 


Fig. 3. Fuzzy logic rules 


The watershed segmentation algorithm is used in computer vision to recognize objects. It was originally developed 
for cartography. And now it perceives images as a topographic landscape, where the grayscale intensity corresponds to 
the terrain elevation. Watershed lines outline individual objects as boundaries between drainage basins. By filling the 
landscape from local minima, the algorithm separates objects by the intensity of gray color. 
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The watershed detection algorithm is effective for some applications, but it can be sensitive to noise and may cause 
over-segmentation. In practice, careful preprocessing and marker-driven watershed options are used to improve 
performance [11]. 

The watershed algorithm output is shown in Figure 4. 
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Fig. 4. Watershed algorithm output: 
a — original image; b — result of watersheds 


Research Results. After using the watershed algorithm, boundaries were drawn for the points under study [12]. To 
construct reference lines, the extreme points were calculated along the boundary of the mark (selected region) of the bone. 
Then, the center of mass was calculated, the image was converted to binary format, and its center was found. Having 
determined the extreme points on four sides, we took the middle of the top of the upper bone, the middle of the base of 
the lower bone and drew a line. 

We used the formula for calculating the angle between two lines and three coordinate points: aa, bb and cc. Two are 
the extreme points of both lines, and one is the common point of their intersection [13]. 


b=-axuit+y, (2) 
x1 =i —b+a, x2 =y2—b+a, (3) 
a =1.0x(y2-y1)+1.0x(x2—%1). (4) 


Figure 5 shows the result of constructing reference lines and angles on X-ray images. 
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Fig. 5. Results of using the improved watershed algorithm on X-ray images: 
a — construction of reference lines; b — calculation of angles 


The performance of the Canny edge detector [14] was compared to the results of the new algorithm proposed by the 
author. The higher performance of the latter solution is obvious. This is explained, firstly, by the improvement of the 
Canny operator. Secondly, the new method is integrated with an improved Gaussian filter. Therefore, there is no need for 
an additional filter, and the execution time of the algorithm is reduced. Its implementation has significantly increased the 
accuracy and yielded better results, specifically for images of radial lines with high noise levels. 
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The new solution was implemented at Al-Basel Hospital in Latakia, Syria. Images from two X-ray machines were 
used. The program was applied to X-ray images of 500 patients. In 490 cases, no manual processing was required, i.e., 
the accuracy of the results was 98%. Image processing took about 20 seconds. 

The proposed algorithm is implemented in Python with a graphic interface. The scheme for creating reference lines and 
angles for processing X-ray images of the knee joint is designed to improve the interpretation of X-ray images (Fig. 6). 


Input image 


Canny algorithm 
(automatic calculation 
of T1, T2) 


Morphological 
operations 


Calculation of the centers 
of mass of the knee joint 


Watershed 
algorithm 


Calculating points of each 
knee contour 


Caclculation and 
determination of reference 
lines and angles 


Fig. 6. Improved algorithm for creating reference lines and angles for processing knee joint X-ray images 


Discussion and Conclusion. A new approach to knee X-ray imaging is being used at Al-Basel Hospital (Latakia, Syria). 
The practice has shown the value of advanced automated image analysis in improving diagnostic accuracy. It should be noted 
that even an earlier version of the program effectively identified bone abnormalities that were missed by visual analysis [15]. 
With the help of the proposed solution, orthopedists can qualitatively work through subtle factors indicating violations, as well 
as combinations of indicators that can cause confusion and hinder the detection of pathology. It follows that the use of the 
approach described in this article allows us to reasonably expect an increase in the quality of medical services in the field of 
orthopedics and surgery. Specialists receive tools for more precise planning and execution of operations. In addition, the new 
solution opens up opportunities for better organization of conservative treatment taking into account the individual 
characteristics of the pathology. 
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Scientific research in this line should be continued to expand diagnostic capabilities and a range of treatment strategies 
for musculoskeletal diseases. Of particular interest is achieving absolute, 1.e., 100% accuracy in algorithmic determination 
of the lines. The time for image processing should also be reduced. Optimization of computational processes will provide 
obtaining results faster, which is particularly important in emergency diagnostics. In addition, it is advisable to refine the 
method for analyzing more complex clinical cases. Another promising direction is adapting the algorithm for analyzing 
other joints, such as hip or shoulder joints. 
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